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Summary. — D,D neutron polarization has been established again by 
measuring the azimuthal dependence of the differential cross-section 
of C2, Experiments allow assignment of a lower limit 20% for the pola- 
rization of 3,32 MeV neutrons, emitted at 45° from the deuteron direction 
for a deuteron energy of 610 keV. Measurements of the angular distri- 
bution of scattered neutrons at resonance energy E,„ = 2,08 MeV permit 
to identify .the level as a D state (J= 3/2 or 5/2). At the energy 
E, = 2,20 MeV, i. e. far from resonances, the angular distribution has 
the expected theoretical behaviour for the potential scattering. The 
effect of polarized neutrons on the angular distribution near an isolated 
level of the D,, type has been calculated. 


4. — Introduction. 


Recently it has been experimentally established that neutrons emitted 
in the D(d,n)He? reaction are polarized ("), (), as theoretically predicted 
and discussed by different authors (?)(‘). The polarization of protons emitted 
in the competing reaction D(d, p)T was measured by BisHop et al. (5). 

This work was undertaken with the purpose of studying polarization effects 
in the neighbourhood of an isolated resonance level. The azimuthal dependence 
of the differential cross-section in this case may give the magnitude of the 
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polarization of incident neutrons as one may use the Breit-Wigner scattering 
phases for the theoretical evaluations (°) (°). 

Because of the smallness of the effect, however, it is necessary to work 
with a neutron beam forming an angle « of about 45° or 135° in c.m.s. with 
the direction of the incident deuterons, since for these angles it is expected 
that the neutron-polarization P reaches its maximum value. 

The resonance level at EH, = 2,08 MeV of C" (Fig. 11) fulfills the above 
requirements, but unfortunately as its width is /’<10 keV no measurement 
of polarization effect is possible using D,D neutrons, whose spectral width is 
in every case several times 10 keV. At this energy we therefore limited 
ourselves to the measurement of the mean differential cross-section and 
established that the resonance corresponds to a D level. 


2. — Experimental method. 


We use a voltage of 610 keV (5) to accelerate a 15 uA current of deu- 
terons against a heavy ice target whose thickness corresponds to 0,25 of the 
saturation value. The neutrons are detected with anthracene scintillation 
counters previously described (°). 

The neutron source is completely surrounded by blocks of paraffin mixed 
with 15%, H,BO, wich forms a cube of about 75 cm a side. It is possible to 
rotate in a vertical plane about the neutron source a portion of the paraffin 
screening provided with a channel of 2,4 em diameter (1°). The scatterer 
consists of a graphite cylinder of the same diameter placed at 15cm from 
the end of the 60 cm long channel and can be in a very reproducible way 
placed or removed from the neutron beam that emerges from the channel. 

For each angle ß of detection, measurements with and without scatterer 
are repeated 10 times, About 10° counts are registered each time and this 
takes about a minute; one has therefore 104 counts for each fp. If 20% of the 
counts are due to neutrons scattered from graphite, this gives a statistical 
error of 2,2 %. Of course these are mean values as the counting rates 
depend upon the angles «x and f of Fig. 8. 

The most important fluctuations of measurements arise from the variation 
of the background. For « = 45° the background is 3% of the direct beam, 
but rises to 8% for « = 130°. 


) J. V. LEPORE: Phys. Rev., 79, 137 (1950). 
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(%) R. Rıcamo and W. Züntı: Helv. Phys. Acta, 24, 419 ( 951). 
1) R. MEIER and R. Ricamo: Helv. Phys. Acta, 26, 430 (1953). 
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3. — Theory. 


3-1. Unpolarized beam. — The differential cross-section o(ÿ) for an un- 
polarized neutron-beam may be written according to BLOCH (1!) 


il) = 
Di A 4), + LAS 12] Pılcos 9) + (LZ + 1er + Lye? s)21P,(cos 8) ir 
+S [Ava — Ai? ulPikeos 8) + Min — ne „1P\(cos 9) | , 
where we put: 
(2) aj. = exp [ide] sin O31) » 
(3) Mie = EXP [66717] sin Ortaja + EXP [256% 172] DR ein CNE 


E, is the neutron energy, E”,, the true resonance energy B,, T is the width 
of the level, k=1/k. The P, are the Legendre polinomials and P} the first 
associated Legendre functions defined as P} = dP, (cos 0)/dd. The ae 
Al 12 represent the contribution of the potential scattering alone and as 
the total eontribution from the potential plus the resonance scattering. 

If for simplieity we suppose that the potential scattering is not spin de- 
pendent, we have 


(4) Atta = aa A a = À, 
Putting 
D} 
5 iy = = rn —  — exp [ty] sin 
(5) ne ne 0 p [ty y 


and supposing that only one level is effective, we write (1) in the form 


(6) ko(9) =|a|?+ |B)? 5 


ee à = exp [2iö, u, P:(cos 9) , 


a = > (21 +1) A,Pı(cos d) bic == 
B = + iexp [2i0,]uP;(cos 9) , 


(1) F. Brocu: Phys. Rev., 58, 829 (1940). 
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ö means the phase due to potential scattering. The total phase 6 is 


(8) d=d+y, 
where 
E,— E, 
(9) y = are cotg - TR (see (5)). 


Equation (6) gives the mean value of the angular distribution averaged 
on all azimuthal angles, and, of course, this is independent of the polarization 
state of incoming neutrons. 

In Fig. 1 and 2 are given the curves calculated from (6) taking = — 90°, 
6,=— 140, 6, — 0° and y= 90°. It appears there that the angular distri- 
butions are very dependent upon the parity of the resonance level. 


un 


Fig. 1. — Potential scattering alone, re- Fig. 2. — D,,; Dsj. and F,, levels at 
sonance levels P,, and D,, at resonance. resonance. 

Fig. 1. and 2. — Differential cross-section o(@) in units of X? for scattering phases 
do = — 900, = — 149, §, = 5, = 0°, = y = are cotg (E,— E,)/1/2 = + 90° for 


an unpolarized neutron beam. 


3:2. Polarized beam. — If there are n, neutrons with the spin perpendicular 
to the scattering plane in one direction and n_ antiparallel to these, the relation 


(10) pass TER 
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gives the magnitude of the polarization vector perpendicular to the scattering 
plane of the D(d, n)He? reaction. 

The differential cross-section for a polarized neutron beam in function 
of Ÿ and @ is given (*) by: 


(11) k?o(d, p) = k?o(8) + P sinp[lo + Bl— (|x|? +1819]- 
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Fig. 3 and 4, — Differential cross-section in units of A? for a polarized neutron beam 
with 100% polarization; 5) = — 90°, 5, = — 14°, 5, = 0°; @ = 0°, 90°, 270°; L = 2; 
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In Fig. 3 we give the behaviour of (11) for the case of a D,, level, at 
neutron energy H, = E,— 1/2, and g = 0°, 90°, 270°. In Fig. 4 we give 
the analogous curves for #, = E, + I/2. 

It is interesting to note (see (11)) that 
there is no polarization effect if: 1)«=0; 
2) 6=0; 3) the phase difference between « 
and B is equal to (2n + 1)z/2 with n integer. 
Case 2) oceurs at Ÿ = 90° for isolated levels of 
\ even parity (as in the case chosen). In conse- 

-0,5 oF se quence it is not convenient to make measur- 
Fig. 5. — Representation of a ements of polarization effects at the fixed 
complex quantity of the type angle Ÿ — 90° only. We note also that in 


e'v sin y for graphical compu- any case the polarization effect is zero at 
tation of o(#) and o(#, y). Ÿ = 0°, 180°. 


We found it very convenient to make gra- 
phical calculations as all terms that enter in (11) are of the type e”* sin y 
(Eo) 
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4, — Experimental angular distributions. 


4:1. Potential scattering. (H, = 2,20 MeV, «=120°). — At this energy, 
we are nearly three times the width of the neutron spectrum above the 
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Fig. 6. — Experimenta] points for c() Fig. 7. — Experimental points for o(@) at 
at E,= 2,20 MeV, « = 120°, and theo- E, = E, = 2,08 MeV, « = 130°,5 and 
retical curve for potential scattering theoretical curve in arbitrary units, sup- 
with §) =— 90°, 6,=— 14%, 5, — 0° posing a D,, resonance level; 5,= — 90°, 
(arbitrary units). On = 14, 5, = 09, 65 =w= + 90°. 


level at 2,08 MeV. The experimental points are given in Fig. 6. The curve is 
theoretical and corresponds to (6) with 6, = — 90°, 6, = — 14°, 6, — 00. The 
theoretical cross-section 


(12) Ou = Anak? ¥ (21 + 1) sin? 6,, 
l 


amounts up to 1,61 barns; its experimental value is 1,65 barns. 
The measurements didn’t provide information on polarization because the 
background was very large (x = 1200). 


42. Resonance level at E, = 2,08 MeV (x = 1300,5). — The experimental 
point and the theoretical curves are given in Fig. 7. We take here 6, = — 90°, 
6, =— 14°, 6, = y = 900. 
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The resonance level belongs clearly to a D state. 


It is not possible to 


decide if the total angular momentum is J= 3/2 or 5/2 (?). 

The experimental anisotropy is less pronounced than the theoretical one 
because the level width [= 10 keV, while the spectrum width A = 30 keV. 
For this reason, we did not succeed in detecting a polarization effect here. 


43. Azimuthal dependence of the differential cross-section. (H, = 3,32 MeV, 
& = 45°). — Neutrons of 3,32 MeV are obtained in our experiments with 610 keV 


deuterons, at an angle x = 45°. 
polarization experiment. In addition 
at this angle the background is only 
3% of the incident beam and we are 
then favoured for the measurement of 
the small polarization effect. 


B 


gy =90° 8) 


Fig. 8. — Experimental arrangement for 
o(#, y) measurements. D is the detector that 
receives neutrons scattered from graphite 
G at an angle Bin lab. system (9 in cm:s.). 
Position I corresponds to the old experi- 
ments (1¢); position II to that of the new 
ones. Azimuthal angle p= 270° always cor- 
responds to neutron scattered inside the 
cone with 2 axis. 


This is the optimum value required for a 
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Fig. 9. — Azimuthal dependence of 


differential cross-section at E,— 3,32 
MeV, «= 45°. Dotted curve joins ex- 
perimental points for p = 90°, conti- 
nous curve that for y = 270°. Mea- 
surements are done in position II 
of Fig. 8. Polarisation effect is clear- 
ly outside of the experimental 
errors. 


As we didn’t succeed in detecting polarization effects in experiments 41 
and 42, we repeated the measurements at Æ, = 3,32 MeV wich were pre- 


viously described (1). 


In order to prove that no apparative asymmetry was 


responsible for the observed polarization effects, we measured the azimuthal 
dependence of differential cross-section in a d-n scattering plane rotated by 
180° around the direction 2 of incident deuterons (Fig. 8). 


(2) C. K. BOCKELMANN, D. W. Mrtter, R. K. Apatr and H. H. BARSHALL: Phys. 


Rev., 84, 69 (1951). 
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Experimental results appear in Fig. 9. We see that the intensity of scat- 
tered neutrons at 90°< 9 < 120° is larger at p= 270° than at 9= 900. The po- 
larization effect has here the same beha- 


10 = viour as in the first experiment (1). 

= 44, Resonance at H, = 3,66 MeV 

NS (z= 0°). - We checked our experimental 
| arrangement by making à measurement 
sl of the differential cross-section at x = 00. 


The experimental points fit very well with 
the old ones (Fig. 10) obtained by an 
entirely different geometry (1). 


a 


(e) 30 60 90 12051501807 


Fig. 10. — Differential cross-section 

at H, = 3,66 MeV, & = 0°. Experi- 

mental points: © present measure- 

ments; @ old measurements with the 
ring method (13). 


5. — Conclusions. 


With our new measurements | E, 
0 2 3 Mey 4 
we are able to confirm the pre- 
sence of polarization effects with Fig. 11. — Total neutron cross-section of C12- 


neutrons of the D, D reaction. from ref. (12), 


We can, however, still assign only 

a lover limit to the polarization degree, P > 20%, since the phases of the 
partially scattered waves are still unknown for an energy of 3,32 MeV 
where à polarization effect was successfully stated. The tentative phase 
analysis of HUBER et al. (14) cannot indeed be safely used at this energy. 
Another result of this work is the assignment of { — 2 to the level at 
2,08 MeV. 


We wish to thank Prof. P. SCHERRER for his kind and helpful interest and 
Prof. M. VERDE for his suggestions and encouragement of this work. 


(3) R. Ricamo: Nuovo Oimento, 8, 893 (1951). 
(4) P. HUBER, E. BALDINGER and R. Buppe: Helv. Phys. Acta, 25, 444 (1952). 
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RIASSUNTO 


Usando neutroni D, D viene misurata la sezione d’urto differenziale o(#) del CY 
per energia dei neutroni Æ, — 2,20 MeV, cioè lontano da ogni risonanza. L’aniso- 
tropia constatata viene interpretata come effetto del contributo di onde di momento 
orbitale 1 = 0 ed 1=1. Niente pud essere detto sull’effetto di polarizzazione per il 
«potential scattering ». Viene misurata la o(®) per HE, = 2,08 MeV corrispondente ad 
una risonanza della o,,, nel C1? (vedi fig. 11) con una larghezza del livello /=10 keV. 
I risultati delle misure confrontati con le curve calcolate teoricamente permettono 
di affermare che il livello € D,, oppure D,,. Poiché la larghezza dello spettro dei neu- 
troni impiegati era di 30 keV, & stato impossibile mettere in evidenza effetti di po 
arizzazione per questa risonanza. Sono state ripetute le misure degli effetti di po- 
larizzazione (14) ad Æ, = 3,32 MeV ma per un azimut dei neutroni emessi nella reazione 
D, D aumentato di 180°. Si confermano i risultati delle prime esperienze. La polariz- 
zazione dei neutroni D,D per E, = 610 keV e angolo di emissione nel labora- 
torio « — 45° risulta P >20%. Vengono calcolate le o(d, pP) secondo LEPORE (°) € 
VERDE (7) in funzione di 9, assumendo come parametro l’azimut g. Si mostra cosh 
l’'andamento della sezione d’urto differenziale ottenibile in presenza di effetti di pola- 
rizzazione. 
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Studio dell’azione di vari metodi di sviluppo nella formazione 
del fondo delle emulsioni nucleari. 


G. ALVIAL C. 


Departamento de Fisica, Facultad de Filosofia y Educacion - Universidad de Chile 
(ricevuto il 3 Agosto 1953) 
Riassunto. — Si investigano le leggi di formazione del fondo (fog) delle 


emulsioni nucleari di varia sensibilita in funzione della tecnica usata nel 
rispettivo processo di sviluppo. 


1. — Nel processo di formazione dell’immagine latente nelle emulsioni nu- 
cleari bisogna considerare, in primo luogo, la reazione fotolitica: 


(1) Br- + 4H > Br +e 
Agt—+ e— Ag. 


La reazione anteriore € valida qualunque sia il processo di formazione del- 
l’immagine latente nella rete cristallina degli alogenuri d’argento (12). EB oppor- 
tuno far risaltare che il AH della reazione fotolitica rappresenta un’energia 
di ionizzazione causata da diversi agenti: perdita di energia delle particelle 
cariche che attraversano l’emulsione, processi foto-ionizzanti, variazioni ter- 
miche, pressioni meccaniche, reazioni chimiche. 

Alla formazione del fondo delle emulsicni nucleari (fog) contribuiscono 
soprattutto i processi di manifattura, agenti termici, meccanici e chimici. 

L’autore di questo lavoro studia la formazione del «fog» in funzione dei 
fattori fisico-chimici corrispondenti ai diversi metodi di sviluppo delle emul- 


(1) R. M. Gurney e N. F. Morr: Proc. Roy. Soc., A 164, 151 (1938). 
(2) J. W. MITcHELL: Bull. Soe. Roy. Sci. Liège, 20, 300 (1951). 
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sioni: in altre parole, stabilisce condizioni sperimentali tali che il AH di (1) 
venga proporzionato preferibilmente dagli agenti fisico-chimici dei diversi 
metodi di sviluppo. 


2. - Le emulsioni nucleari usate furono di fabbricazione Eastman Kodak 
dei tipi NTA, NTB e NTB3 di 200 y di spessore. Soltanto in un esperimento 
si usarono emulsioni NTA di 150 u di spessore. 

Gruppi di lastre dei diversi tipi furono sottoposti a quattro metodi di svi- 
luppo: a) quello indicato da BLAU e DE FELICE (*); b) il metodo di « hot-plate » 
con sviluppatore D19b (metol-idrochinone) e con i tempi indicati da DAINTON, 
GATTIKER e Lock per emulsioni di 200 y di spessore (*); €) sviluppo in azol 
(paraminofenolo) identico al metodo indicato dagli autori precedentemente 
citati per emulsioni di 200 y: di spessore; e, d) un processo di sviluppo simile 
a quello che consiglia la fabbrica Eastman Kodak per lastre del medesimo 
spessore (5). Sostanzialmente ultimo processo è analogo a quello descritto da 
WILsoN e VANSELOW (*): D19 (elön-idrochinone) diluito di 1:1 a 5°C per 
30 minuti e D19 di 1:3 a 20°C per 30 minuti. Indicheremo ognuno di questi 
processi con i suoi nomi usuali: il primo come sviluppo in «due bagni >», il 
secondo, «hot-plate », il terzo, sviluppo in «due temperature » con azolo, € 
l’ultimo in «due temperature» con D19. 

Tutte le lastre furono sottoposte a un pre-lavaggio per 30 minuti; dopo, 
al corrispondente processo di sviluppo; poi, ad un bagno contenente acido ace- 
tico al 12%; il fissaggio si effettud in una soluzione di tiosolfato di sodio al 40% 
e, prima che le lastre si asciugassero, furono immerse 


per un’ora — in una 
soluzione di glicerina al 10%, come indica MIGNONE (?). Le lastre vennero 
asciugate senza corrente d’aria. 

Le emulsioni furono osservate per immersione in un mieroscopio Univer- 
sale Galileo con un ingrandimento di 3700 diametri, corrispondente ad una 
superficie osservabile di 1344 wy’. 

Le misure si realizzarono su un totale di 50 lastre. 


3. — In ogni caso & stata misurata la densita dei granuli del fondo dell’emul- 
sione (fog) in funzione della profondita. Tanto il primo strato quanto Pultimo 
misurato si presero a 5 u dalla superficie e dal vetro rispettivamente, allo scopo 
di evitare nelle misure il «fog» originato da cause estranee a quelle del pro- 
cesso di sviluppo. 


M. Buau e J. A. DE Ferice: Phys. Rev., 74, 1198 (1948). 

A. D. DAINTON, A. R. GATTIKER e W. O. Lock: Phyl. Mag., 42, 396 (1951). 
Bollettini della Eastman Kodak. 

M. J. Wuson e W. Vanserow: Phys. Rev. 75, 1144 (1949). 

G. MiGNoNE: Nuovo Cimento, 8, 896 (1951). 
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In ogni lastra si misurarono le densità corrispondenti per lo meno à quattro 
strati. I diversi grafici danno il numero di granuli per ogni 1344 u? di super- 
ficie osservabile in funzione della 

profondita. 
Si noti che solamente alcune 
emulsioni conservano lo spessore 
— ‘granulo primitivo. Nel caso delle NTB, 


2 
3 


‘ 
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3 ys curva b della fig. 3, la soluzione 
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Fig. 2. 


considerare che il tipo di emulsione & differente da quello usato da MIGNONE 
nel suo lavoro. 

Le curve a, b, e e d della fig. 1 corrispondono ad emulsioni del tipo NTA 
sviluppate in due temperature con azolo, in due temperature con D19, in due 
bagni a « hot-plate » rispettivamente. Ogni curva rappresenta le misure medie 
di tre lastre, eccezione fatta della curva d, che rappresenta soltanto le misure 
medie di due lastre. I valori particolari di ogni lastra sono segnalati nel gra- 
fico. Tutte le lastre osservate sono state scelte da un medesimo gruppo, sono 
di identica manifattura e sono state sviluppate in un periodo di 60 giorni, 
nell’ordine seguente: gruppo b, gruppo €, gruppo wu e gruppo d. Le condizioni 
di magazzinaggio in questo frattempo sono state normal. 

La fig. 2 rappresenta il numero di granuli per 1344 u? in funzione della 
profondità per due lastre tipo NTA delle quali una è stata sottoposta a 
«fading » in una soluzione al 3 %, di CrO, (5) e l’altra non & stata sottoposta 
a questo processo. Esse sono rappresentate dalle curve a e b rispettivamente. 


(5) N. A. PERFILOV: Compt. Rend. Acad. Sci. URSS, 42, 258 (1944). 
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Ambedue sono state sviluppate in due temperature con D19. Sostanzialmente 
appartengono alla famiglia corrispondente alla curva b della fig. 1. 
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La fig. 3 riassume i risultati degli esperimenti ottenuti per emulsioni del 
tipo NTB. La curva a é la curva media di misure realizzate su due lastre svi- 
luppate in due bagni; la curva b, di due lastre sviluppate in due temperature 
con D19; e la €, di due lastre sviluppate con il metodo di « hot-plate ». Tutte 


sono di una medesima serie, di 
uguale manifattura e sviluppate 
in un periodo di tre giorni. 


KH 2 granuli 


2 


r= 10.6 


—+ 1 granulo 


numero di granuli per 1344 u 


Fig. 4. 


Le curve a, b, ¢ della fig. 4 sono le curve medie di due lastre ognuna, del 
tipo NTB3, sviluppate rispettivamente con i metodi « hot-plate » in due bagni 
e in due temperature con D19. 

Infine, le curve a e b della fig. 5 rappresentano curve medie di misure 
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realizzate in tre lastre, ognuna dello spessore di 150 u, del tipo NTA. Il 
gruppo a venne sviluppato con il metodo di «hot-plate», ma negli ultimi 
cinque minuti dello sviluppo, la temperatura si elevd a 40 °C, di modo che 
la superficie delle emulsioni si asciugd completamente. Il gruppo b non passd 
al «plate », ma rimase un tempo equivalente nello sviluppo a 5 °C. 

In analoghi esperimenti realizzati con le rimanenti 19 lastre, é stato otte- 
nuto un accordo totale con le precedenti per le curve densita-profondita delle 
emulsioni del tipo NTA e NTB; e solo parzialmente per il tipo NTB3. Questo 
accordo non consiste tanto nei valori numerici (sebbene essi non si allonta- 
nino sensibilmente dai casi analoghi precedentemente descritti), quanto nella 
forma di variazione della densité del «fog» in funzione della profondita. 


4. — In ogni strato si determino la distribuzione della fluttuazione statistica 
della densità granulare, che risultO gaussiana. I pochi casi di fluttuazioni lon- 
tane dalle normali, corrisposero esattamente a misure di fenomeni spuri 0 
malamente realizzati (isolotti affetti da radioattività, osservazioni in cattive 
condizioni ottiche, ecc.). 

Nelle lastre rappresentate dalla curva a della fig. 1 si osserva un «fog » 
costituito da granuli d’argento di diametri normali (0,4 u) ed un altro fondo 
formato da innumerevoli e finissimi granuli d’argento. Nella detta curva inter- 
viene soltanto il primo tipo di «fog»: percid essa non riflette fedelmente il 
fenomeno in istudio. 

La densità granulare negli altri tipi di lastre sviluppate con azolo é notevol- 
mente elevata e tale da impedire il passaggio della luce per trasparenza all’os- 
servazione microscopica. Per questo motivo non appariscono nei diversi grafici 
le curve corrispondenti al trattamento con azolo. 

Nella fig. 2 si confrontano le curve densita-profondita di due lastre, di cui 
una — lastra b — & stata sottoposta al «fading» e l’altra no. Ambedue sono 
state sviluppate contemporaneamente in due temperature con D19 esatta- 
mente un mese dopo lo sviluppo della serie analoga rappresentata dalla curva b 
della fig. 1. Esse appartengono sostanzialmente alla medesima famiglia della 
curva b della fig. 1, nonostante la differenza di tempo nello sviluppo. Si 
osserva inoltre che la curva della lastra nella quale & stato effettuato il « fading » 
non si allontana notevolmente dall’altra. Il medesimo risultato & stato otte- 
nuto in esperimenti analoghi realizzati con emulsioni del tipo NTA e sotto- 
poste ad altri metodi di sviluppo. 

Un gruppo di lastre NTA di 200 » venne sottoposto a riscaldamento fra 
30 °C e 83 °C prima del rispettivo processo di sviluppo. Sviluppate in due tem- 
perature con D19, le curve densità-profondità non mostrarono ditferenze fon- 
damentali con la curva b della fie. 1. 

Questo trattamento, come quello del « fading » indica chiaramente che nella 
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formazione del fondo delle emulsioni nucleari intervengono direttamente i fat- 
tori fisico-chimici del metodo di sviluppo e quelli della stessa emulsione. 

Diversi autori hanno fatto riferimento a detti fattori in lavori di altro 
carattere. I fattori che ei interessano sono: 1) variazione del pH della solu- 
zione sviluppatrice nell’emulsione; 2) durata di penet razione dello sviluppa- 
tore nella emulsione; 3) conduttivita termica dell’emulsione; 4) natura chi- 
mica della soluzione sviluppatrice. 

Il fattore 3) rimane chiaramente esposto nell’esperimento rappresentato 
nella fig. 5. La curva b indica che la densité del «fog » delle emulsioni svilup- 
pate a 5 °C non è significativa, in cambio la curva @ mostra un massimo ben 
definito. Cid pud essere interpretato nel seguente modo: le lastre collocate 
nel « hot-plate » a 27 °C si svilupparono uniformemente, all’elevarsi della tempe- 
ratura a 40°C durante 5 minuti e all’evaporare del liquido sviluppatore dalla 
superficie delle emulsioni, l’azione sviluppatrice si intensificd negli strati inter- 
medi dell’emulsione dove ancora restava soluzione sviluppatrice, ma non negli 
strati vicini al vetro, che per ragioni di conduttivita termica non riuscirono a 
raggiungere la temperatura di 40 oC. 

Un’analisi della fig. 1 indica cid che segue: la retta quasi parallela all’asse 
delle profondita — retta d — mostra che lo sviluppatore penetro totalmente 
e che tutta l’emulsione si riscaldd ad una medesima temperatura. Pereiö, i 
fattori 2) e 3) non spiegano la diminuzione delle curve be ce. Posto che queste 
rappresentano i processi eseguiti con i metodi di Wilson-Vasenlow e Blau-De 
Felice, metodi che quasi non differiscono chimicamente, rimane soltanto il 
fattore 1). D’altro lato, le curve di densita costante per il primo di questi 
metodi in emulsioni del tipo NTB e NTB3 e le curve di densita sempre de- 
crescente corrispondenti ad emulsioni dei medesimi tipi sviluppate con il me- 
todo di Blau-De Felice inducono a far notare che, per il metodo in due tem- 
perature con D19, esiste una variazione del pH nelle emulsioni NTA.e che 
non esiste siffatta variazione negli altri tipi di emulsioni sviluppate con questo 
stesso metodo. D’altra parte, in. base a questa medesima analisi si pud accet- 
tare un’azione sviluppatrice della soluzione A nel metodo in due bagni per 
ogni tipo di emulsione. 

Le curve di densité costante corrispondenti ad emulsioni del tipo NTA e 
NTB sviluppate con il metodo di « hot-plate » e di densité decrescente per il 
tipo NTB3 sviluppate con il medesimo metodo, indicherebbero una variazione 
della conduttività termica per quest’ultimo tipo di emulsione. 

Tanto la variazione del pH, quanto la variazione della conduttivita termica 
si possono spiegare a causa della differenza nella granulazione delle emulsioni. 

Infine, la maggior densita granulare del «fog» per le emulsioni NTB3 è 
funzione della sensibilità delle suddette emulsioni. Detta sensibilità spieghe- 
rebbe anche le variazioni parziali delle curve densità-profondità in emulsioni 
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del tipo NTB3. La grande densità che si nota in tutte le emulsioni trattate 
con azolo & dovuta all'influenza del fattore 4). 


5. — Attualmente & in corso in questo laboratorio una ricerca per determi- 
nare la correlazione esistente tra i precedenti risultati e la densità granulare delle 
tracce nucleari. Per cause segnalate da altri autori (°), si lavora con prımari 
di grande energia della radiazione cosmica. 
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SUMMARY 


The Author studies formation laws of fog in nuclear emulsions of different sensi- 
bilities, as a function of the correspondent development method. 
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Determination of the Age of the Earth 
from the Isotopic Composition of Meteoritic Lead. 
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Summary — Under the assumption that a) the isotopic constitution of 
lead at the time of the formation of the lithosphere is represented by 
the figures found for lead from the troilite phase of the Cañon-Diablo 
meteorite [Pat 53] and b) the majority of tertiary lead ores analysed for 
isotopic constitution has been formed according to the simple model of 
primary origin (equ. 3), the age of the lithosphere w follows as w= 4,5 + 
+ 0.3:109 a. 


Several authors [Ger 42, Hol 46,47, Hou 46,47, Bu 49, Co 52] have cal- 
culated the time w since the formation of the lithosphere from the isotopic 
abundances of common lead, using the data given by A. O. Nıer [Ni 38] and © 
additional data, published by COLLINS, FARQUHAR, RUSSELL [Co 52}. All 
these calculations are based on the following simplified model for the for- 
mation of a lead ore or lead mineral. If we denominate the ratio of the 
numbers of Pb2°* and Pb?” atoms to the number of Pb? atoms in given 
* Pb-mineral as «, ß respectively 


(Pb?) (eps) = 
(1) (Pb2°1) en (Pb?°*) =p 


and if we assume the value of these ratios ‚at the time w to be universally 
ys Bw, and w is or would be the present geochemical ratio 


. 5 U238 
(2) Ht = pp 


105 — Il Nuovo Cimento. 
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in the particular rock or rock component from which the lead ore has been 
formed, then the amount of Pb?% and Pb?°7 present at the time p of the for- 
mation of the lead mineral as à function of p is given by 


(3a) C= Clap + u (ei Ber er») 
(35) p= Boe ME (erw A»), 
pcan! 


Equ. (3) gives in an «, # diagram the «lines of lead development » for constant 
values of « as a function of the time. By dividing (3b) by (3a) the equations 
of «isochrones » are obtained 


p— D ve il erw — er? 
Gao TER) CG | 


(4) 


These are a number of straight lines, intersecting at the point &,, P„ cor- 
responding to the isotopic constitution of «primeval lead» at the time w. 
Principally it is possible to extrapolate these isochrones graphically or to 
calculate the values of &,, B,, from the «, 6 values of 3 different lead samples 
of at least 2 different ages p. This has been done by the authors mentioned 
above, leading to the values given in Table I. The values for w obtained by 
all the authors coincide fairly closely, thus leading to the generally accepted 
value (*) of the earth’s age 


w = 3,3 + 0,3-10° a. 


TABLE I. 

| | | | 
| Author 4 | a | B y 
| (109 a) a ® Mf 
| SS Kr GEREING [Ger 22] Er a eee | 3.23 | — | os — 
IPASHOEMESNIAOAGIRE EEE 50) 12 502 SNS TRS 
| F. G. HourTEerRMANS [Hou 46, Hou 47] 2.9 11.52 14.03 SE 
| AS HOLMES [Hol AT INR 3.35 10.95 | 13.51 | — 
i. 0 BUBLARDI Ela | BURAIT Er meee INES OSs SOs _- 

Chis, Cones Oras (CS I © 20) EE | 11.33 13.55 31.10 


(*) G. Gamow: Dan. Mat. Fys. Medd., 27, Nr. 10 (1953); P. JORDAN: Schwerkraft 
und Weltall (Braunschweig, 1952); Address by Pope Pius XII to the Pontifical Academy 
of Sciences on November 22, 1951 (Tipografia Poliglotta Vaticana, Rome, 1951); C. F. 
VON WEISZÄCKER: Die Geschichte unserer Welt (1951). 
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It has been pointed out by the author [Hou 47, Br 47] that there are strong 
indications to find the Pb contained in iron meteorites to be very close to the 
isotopic constitution of primeval lead because the chemical ratio U/Pb in 
iron meteorites from then available data [Go 38] seemed to be 1.8:107?, i.e. con- 
siderably less than for instance in igneous rocks from which most Pb minerals 
are supposed to be derived. 

Through the courtesy of the authors a very important paper by C. PATTER- 
SON, G. TILTON, M. INGHRAM, H. BROWN [Pat 53] came to our knowledge before 
its publication which is apt to throw an entirely new light on the questions 
involved. Though these authors found the content of Pb in iron meteorites 
to be considerably less than given above, the experimental ratio U/Pb in 
the troilite phase of the Cafion-Diablo Meteorite was found to be as low 
as 1/2000. 


TABLE II. — Isotopic Abundance of Lead in Cation Diablo Meteorite. 
aa, geen). | ee 
| (Pb2%) | (Ph?%) ta (Pb204) 
| ia eee |e | : a 

role nl oss | 9.41 10.27 29.16 
Metal Phase ..... | 9.7 10.5 | 29.3 


The values found by PATTERSON et al. for « and f given in Table II indi- 
cate very strongly that the lead analysed has never been associated for any 
considerable time with a proportion of uranium or thorium of the same 
order of magnitude as in the lithosphere. Of course it is impossible to conclude 
from the present low value of # found by these authors that the value of u 
for the material from which the meteorite has been formed was as low as it 
is found at present since Paneth’s experiments on the He-age of meteorites 
[Pan 52] make it very probable that the age of the iron meteorites range from 
a few millions up to 10% years only. The maximum He-age for the Cañon 
Diablo Meteorite—even if all its Helium is ascribed to radioactive production 
has been given by PANETH et al. [Ar 42] as 5.5-10%a. Yet the low values of a, 6 
found for Pb in both the troilite and the iron phase of the meteorite in itself 
seem to justify at the present state of evidence the assumption that the actual 
values of &,, By are very close to the «, B values found for the troilite phase 
of the Cafion Diablo Meteorite. All experiences on the comparison of isotopic 
abundances of other medium and heavy elements have shown their identity 
of isotopic constitution in meteoritic and terrestric material. We therefore 
have good reasons to assume that the isotopic constitution of « primeval lead » 
at the time of the formation of the lithosphere was very similar to the values 
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found for meteoritic lead. We therefore assume, as an hypothesis. 
(5) i SOL ales 


from the mere fact that such lead does exist. We can now construct the iso- 
chrones for any given value of p, using the geological value for the parti- 
cular age of the galena under the assumption that most lead samples are of 
primary origin. We shall call a lead mineral to be of primary origin if its 
development has followed the evolution model given by equ. (3) from the 
formation of the lithosphere up to time p without change of the value « in 
the parent rock. Assuming specially that the majority of tertiary lead ores 
are of primary origin in the sense given above, we obtain for the inclination 
of the isochrone for all these lead ores 


(6) (ée—É 1 = 0,5793 + 0,004, 


X39 — Ly 


leading to a new value for w from (4) (p & 30-108 a) 


= 
= 


— (4.5 + 0.3)-10° a 


The list of tertiary lead minerals which have been used for the calculation 
of (6) and (7) is given in Table III together with their individual values for 
(P— P.)/(«—%). For calculating our value of w, we prefer to use only iso- 
topic abundances of young samples, since for these the resulting value w is 
exceedingly insensitive against an error in the determination of geological 
age, even if their ages are uncertain up vo 200 million years. 


TABLE III. 
a E | 2 5 
| Number (*) | Mineral and Source | Bb 
| | C7 
| ere 
N11 Wulfenite and Vanadinite, Tucson Mts., Arizona . 0.5851 
N16 Bournonite, Casapalea Mine, Peru Mages 
| 10.5528 
RES acer ae UT [0.5757 
N17 Galena, pe Peru 10.5669 
N21 Galena, sDurang om Mexicog a ES CRE QE En 0.5839 
N23 | Galena, Sonora Mine, Arizona Au à | 0.6014 
N24 | Galena, Freiberg, Saxony Pere gee AR | 0.5924 
C5 | Bolivia =>. 4 EN NY ET BE 0 
ere Vel |. Montenegro, Babin casero Ee Barat Neo 0.5911 
| G14 | Galena, sample 14, Wiesloch, Consens RENE oat | 0.5755 
G15 | Galena, sample 35, Wiesloch, Germany ..... 0.5755 
(*) Numbers given by the authors in original publications. Letters N, C, V, E, G mean 
(Ni38), (Co52), (Vi52), (Eh53), (Gei 53). 
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As has been pointed out by several authors [Hol 47, Je 48, Je 49, Bu 49, 
Gei 52], we find samples which cannot be considered as primary ores in the 
sense given above, since their x, B-values would lead to negative values for 
their «model age» p. These exceptions are the determinations given in the 
Table IV. It is of course arbitrary to a certain extent which samples to exclude 


TABLE IV. 


Number (*) Deposite 
| N7, N8, N14, V23, V24 | Joplin, U.S.A. 
N9 Metalline Falls, Wash. 
C20 Upper Canada Mine, Ontario 
| C22 | Box Mine, Goldfields, Canada 
| C25 Garson Mine, Sudbury, Ontario 
C26 | Worthington Mine, Sudbury, Ontario 
V27 Northern Ural 
V29 | Prijbalchashia, Central Kasakstan 


(*) See footnote from Table IIT. 


and it will be necessary to assemble more statistical data for a more reliable 
calculation of w. Yet the assumption that most lead minerals can indeed be 
considered as primary seems to be justified by the fact that only 8 deposits 
among 87 lead to «negative» model ages. This means that no appreciable 
change in the value of « in the parent rock material has taken place during 
their development from the time w till the time p of the formation of the lead 
minerals. The apparently negative values of p, can, however, be explained 
easily if these ores are considered as having secondary origin. This means 
that a change in the value of w has taken place in their parent material at 
a certain date. It has been shown by VINOGRADOV et al. [Vi 52] and in more 
detail by Geiss [Gei 53] that a change in the value of w can easily lead to 
values «, ß beyond the isochrone for p = 0, if the values of u had been in- 
creased. Indeed if in a parent rock with a given value «, a change of u to w 
has taken place at the time p’ the values «, ß reached at the time p are given by: 


| x= X% | ule er‘) | in (ere — e7”) ; 
1 


(8) | B — Bu + ae (eX el») + gins (er — @') . 


139 139 


It can, however, easily be seen, that even if we do not exclude the samples 
given above as of secondary origin because of their apparently negative « model 
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age» the value of w is not considerably altered. In fact, taking the values 
of «,, Pu given in (5) and assuming for the Joplin galena the average values 
of all its determinations 


& =-21.95 , B= 15.90% 
p = 200 Ma, 


the value of w obtained would be 


Wy opin — 405-109 a. 
This value may be considered as a kind of lower limit of w, resulting if the 
samples with extremely low percentages of Pb?! are to be considered as of 
primary origin. One sees therefore, that an average value of w taken from 
data on all lead samples will hardly be influenced by ascribing primary origin 
to the few samples which are most probably secondary. 

It would be of extreme importance for the verification of the assumed 
model of lead ore tormation, if it were possible to find purely geological evi- 
dence for the hypothesis that the lead ores given in Table IV are indeed of 
secondary origin. 

The new values &,, P, and w (5), (7) are in agreement with the x, 5 values 
not only for tertiary lead samples but also with those of the majority of most 
other dated samples of older age within the limits of error in «, 6 and their 
geological age. Table V shows a comparison of the geological ages with the 
«model ages » calculated from (4) for all dated lead samples, for which the 
isotopic abundance has been measured hitherto. 

In comparing geological ages with «model ages » it must be remembered 
that the geological ages of ore deposits do not necessarily mean the time of 
the separation of the lead mineral from its parent rock material, but in many 
cases may mean the time at which the deposit has reached its present position. 
A discrepancy between « geological age » and «model age » must therefore not 
necessarily mean that a lead deposit is to be considered of secondary origin. 

We have geologically to distinguish between 3 different classes of lead 
deposits: 


1) Those, which may be considered to be of primary origin. For these 
deposits geological age should coincide with the «model age» within their 
limits of error. 


2) Deposits, for which «model age » is inferior to geological age. Such 
deposits, especially those corresponding to a «model age» <0 must be con- 
sidered to be of secondary origin with an increase of # during the development 
between the times w and p. 
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Mineral and Source 


N4 
N5 | 
V21f 
N6 
N7 
N8 
N14 
v23 
v24 
N9 
N10 
N12] 
v22J 
N15 
N18 
N19 
N20 
C10 
C12 
C15 
C16 
C17 
v2 
v3 
v4 
V5 
V6 
V7 
V8 
V9 
V10 
Vil 
V12 
V13 
V14 
V15 
V16 
V17 
V18 


Galena, Great Bear gerne Canada . 
Galena, Broken Hill, N. S. Wales . 
Cerussite, Broken Hill, Re S. Wales 
Galena, Yancey, Col., N.C. 


Galena, Nassau, Germany . 


| Cerussite, Eifel, Germany . 
Galena, Joplin, Missouri 


Galena, Metalline Falls, Wash. 
Cerussite, Wallace, Idaho . 


Galena, Saxony, Germany 


Galena, Tetreault Mine, Canada . 
Galena, Clausthal, Harz, Germany . 
Galena, Przibram, Bohemia . 
Galena, Franklin, New Jersey . . 
Katanga, Belgian Congo 

Frontanac Co., Ontario . 


Rex, Stark Lake, Canada. . 
Sioux Lookout, Ontario 

Primaje DWK . 

Seravezza, Italy 

DWK . ; 

Sangesur, see a 

Northern Kaukasus 

Jakutia . : Sr 
Boleslav, Kan en: 
Northern Kaukasus 

Sadon, Northern Kaukasus 
Armenian Republie . 
Trans-Baikalia . 

Southern Kirgisia Sow EMI: 
Kuramenian Country (Central Asia) 
Northern Kirgisia 

Kirgisian Republie . re 
Kuramenian Country (Central Asia) 
Southern Saboi, Altai 


(*) See footnote from Table III. 
(**) Used values for calculation: a,=—9.41, By=10.27, 
À'=0,9815:107° a. 


Negus Mine, Great Slave ae ade | 


(106 a) 


1400 [Fa 
1500 [Fa 
1500 [Fa 
220 [Fa 


220 [Fa 
220 [Fa 


230 [Fa 


60 [Fa 
60 [Fa 


220 [Fa 


1030 [Fa 
220 [Fa 
220 [Fa 
400 [Fa 
630 [Fa 

1030 [Fa 


1850 [Fa : 
1850 [Fa 53] 


2480 [Fa 


50 [Vi ! 
50 [Vi 
60 [Vi é 
60 [Vi } 
60 [Vi 5 
100 [Vi 5 
130 [Vi 5 


130 [Vi 
130 [Vi 


130 [Vi 5 
180 [Vi 5 
200 [Vi 5 


200 [Vi 
200 [Vi 
200 [Vi 


220 [Vi : 
220 [Vi ! 


w=4.46 * 10° a, 


| Geological Age 


53] 
53] 
53] 
53] 


53] 
53] 


53] 


53] 
53] 
53] 
53] 
53] 
53] 
53] 
53] 


A 


(106 a) (**) 


1600 + 70 
1600 + 70 
1600 + 70 
180 + 70 
380 + 70 
180 + 70 

| 

| 
| negative 

| 
negative 
1400 + 70 
800 + 70 


1 300 + 70 


220 + 70 
500 + 70 
900 + 70 | 

| 1100 +100 | 

| 1100 + 100 

| 2500 + 100 

| 1 700 + 100 

2 700 -+ 150 

| 80 +250 

| 600 +250 
600 -+ 220 
400 + 250 
220 + 250 
160 + 250 
500 + 250 
500 + 220 
700 + 250 

| 400 + 250 

| 600 + 250 

| 1100 + 220 

1000 + 200 | 
500 + 250 
340 + 250 

| 600 + 250 

| 900 + 250 


=0,1535 10° a, 


Model Age p | 
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Nxr:2(%) Mineral and 


V19 West Tan Shan 
V20 | Karatau. 


Source 


V25 Nagolnij Krjash, Da h 


V26 Western Prijbaikalia 
V27 Northern Ural . 
V28 NW Kasakstan 


V29 Prijbalchashia, Central Kasakstan 


V30 Central Kasakstan 


V31 Karelo-Finnian Republic 
V32 Kinon-Kosiala, Finland 3 
Ela Galena Meggen, Gravel Deposite, N any 
Ib NES | 
G24 J Galena Meggen, Germany . 
E2 | Schapbach, Germany . 
G2 Galena Mühlenbach, Germany . . 
G3 
G4 i x re eg F ; $ 
Gs f Galena Christian Lewin Mine, Germany 
G6 | Galena Oberharz HG 687, Germany 
G7 | Galena Oberharz BW 292, Germany . 
G8 | | 2 
Re } Galena Schauinsland, Freiburg, Germany 
G10] ER we 
Gis Galena Liiderich, Germany 
GI2) > 
G13 f Galena Holzappel, Germany . 
+16 
at Galena Rammelsberg, Goslar, Germany | 
G18 | ; ; 
G19 f Galena Ramsbeck, Germany 
G20 | ; À 
Gals | Galena Maubach, Germany 
G22 | Galena Kônigstiel, Germany : 
G23 Galena St. Andreasberg, Harz, ans 
G25 Galena Altmannshausen near Würzburg, 
| Germany . ; 
G26 Galena Stiftsberg Houbronn ere : 
G27 Galena Hertmannsweiler near Stuttgart, 
| Germany . 
G28 | Galena Linsenberg, Ro ous Germanys | 
| | 


G29 | Galena Mitwaba, Belgian Congo 


(*) See footnote from Table III. 
(**) Used values for calculation: 
À'=0,9815-107° a 


Oy=9.41, 


By=10.27, 


Geological Age 
(10% a) 


220 
220 
260 
260 
320 
320 
320 
320 
600 
600 
300 


300 
200 


[= DE ern 


300 


w=4.46 10° a, 


[Vi 
[Vi 
[Vi 
[Vi 
[Vi 
[Vi 
[Vi 
[Vi 
[Vi 
[Vi 


[Gei 


[Grei 


[Gei 
[Gei 


[Gei 


[Gel 


[Gei 


[Gei 


[Gei ! 


[Gei 5 


[Gei ! 


[Gei : 


[Gei 5: 
[Grei £ 


[Gei 5 


52] 
52] 
52] 
52] 
52] 
52] 
52] 
52] 
52] 
52] 


[Eh 53a] 
[Eh 53a] 
[Eh 53b] 


53] 


53] 


53] 
53] 


53] 


Model Age p 
| (105 a) -(**) 


440 + 250 
600 +- 250 
700 + 250 
900 + 200 
negative 
500 + 300 
negative 
300 + 250 
2409 + 280 


| 2300 + 280 


480 +: 70 


380 + 80 
80 + 70 


360 + 50 


340 + 60 


180 -- 70 
120 + 60 


80 + 60 


380 + 70 


340 + 70 


360 + 60 
380 + 80 


200 + 70 
400 + 50 
220 + 50 


160 + 50 
200 + 60 


200 + 60 
200 + 50 


| 1100 +70 


A=0,1535°10~* a 


| 
| 


DETERMINATION OF THE AGE OF THE EARTH BIC. 1631 


3) Lead deposits, for which «model age » is superior to geological age. 
Samples of this type are V13 and N12. It seems that this type of deposits has 
been formed at a much earlier date than indicated by their geological age. 

A full discussion of the application of the formation model (3) and (8) of 
lead minerals and their geological applications shall be given elsewhere [Gei 54]. 
There have been many discussions on the actual meaning of the time w. 
As has been pointed out before [Hou 47], we may consider two limiting cases: 


1) If the cosmochemical value of u, i.e. the value of y for our planetary 
system as a whole is the same as for our lithosphere, the time means the 
«age of uranium » i.e. the time at which elementary processes by which ura- 
nium was produced ceased to exist in our planetary system. 


2) In case the geochemical value of y, which from the evidence of lead 
minerals analysed seems to vary between 8 and 11 in the rock material and 
which may be ~ 10 for the lithosphere as a whole from present geochemical 
evidence, did change considerably at the time of the formation of the litho- 
sphere, w means the time at which this formation has taken place. In all 
other cases, f.i. if there was a gradual increase of u during the time of the 
formation of the lithosphere, w is defined by a more complicated integral. 
The details of these calculations will be given elsewhere [Gei 54]. 


There are at present a number of arguments in favour of the second hypo- 
thesis. The low values of x, B found for iron meteorites as well as their low 
values of # and the well known arguments of the geothermical paradoxon 
seem to indicate very strongly that a considerable enrichment in U relative 
to Pb occurred at the time of the formation of the lithosphere, and we may 
therefore consider w te represent the time at which this increase in y did cecur. 
It seems very satisfactory that the value of w given here, which therefore 
may be considered as the «age of the earth», is not much lower than the 
upper limit of the age of the lithosphere given by ALPHER and HERMAN 
[AI 51] based on the argument that necessarily ©,> 0, which leads to an 


upper limit 
W = ).0 OP 727, 
and the argument of WEFELMEIER and the author [Hou 47] which allows to 


caleulate the «Age of Uranium» as S$ 6.0:10° a assuming that originally the 
ratio U 235/U 238 could not have been higher than 1. 


The author wishes to express his gratitude to Drs. PATTERSON, TILTON, 
INGHRAM, HARRISON BROWN who very kindly sent us the manuscript of their 
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very important paper [Pat 53] to be published in Physical Review. I thank 
Mr. J. ©. HoUTERMANS and Dr. J. GæIss for their helpful assistance in doing 
the calculations and him, Dr. EHRENBERG, Dr. Prcctorro and Dr. UREY for 
many valuable discussions. 


Note added in proof. 


A table of measurements of isotopic constitution of ordinary lead made available 
to us through the courtesy of Dr. HENRY FAUL of the Geological Survey of the United 
States contains more analysed samples which have not been included in the calculations 
in this paper. The figures given in (6) and (7) of this paper, however, are not altered 
within the limits of error by taking into account these new determinations. 

It may be also of interest, that NICOLAYSEN, ALDRICH and DoAk [Fa 53] have 
found a isotope-controlled Strontium/Rubidium-age for the Lepidolite of Letaba 
(Transvaal) of 3.85-10° a, which seems perfectly possible with the new value of w given 
in this paper and would be in contradiction with the older estimates of w given in 
table I. 
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RIASSUNTO (*) 


Nell’ipotesi che a) la costituzione isotopica del piombo all’epoca della formazione 
della litosfera sia rappresentata dalle cifre trovate per il piombo della troilite del me- 
teorite del Cañon-Diablo [Pat 53] e che b) la maggioranza dei minerali di piombo 
terziari di eui & stata analizzata la costituzione isotopica si sia formata secondo il 
semplice modello della creazione primaria (eq. 3), si deduce che Veta w della litosfera 
si ottiene dalla formula (4) w = 4,5 + 0,3:10° a. 


(*) Traduzione a cura della Redazione. 
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On Quantum Field Theory - I: Explicit Solution of Dyson’s 
Equation in Electrodynamics without Use of Feynman Graphs. 


E. R. OAIANIELLO (*) 


Istituto di Fisica dell’ Universita - Roma 
Istituto Nazionale di Fisica Nucleare - Sezione di Roma 


(ricevuto il 24 Settembre 1953) 


Summary. — A compact, explicit expression is given for the term of arbi- 
trary order in the power series expansion of Dyson’s S-matrix, for any 
arbitrary scattering process. The expression thus found for S has a very 
high degree of symmetry (a fact which the Feynman graph technique 
could not reveal; our procedure is also free from several shortcomings 
of this technique), which suggests the existence in the theory of deeper 
properties. These will be studied in the continuation to this work. 


1. — Introduetion. 


1'1.- In this work, the first of a series on questions concerning quantum field 
theory and in particular electrodynamics, we show that the customary method 
for computing the Dyson S-matrix (1), in which use is made of the so-called 
Feynman graphs, can be advantageously substituted with a more compact 
and simpler algorithm. Our method, although equivalent to the first one, 
demonstrates the existence of remarkable properties of symmetry and re- 
cursion in the various terms of the perturbation expansion of the S-matrix, 
which remain completely hidden by using the graph technique. This study, 
althongh essentially of a formal character, is a necessary, preliminary step 


(*) This work has been in part performed at the Institut for Theoretisk Fysik, 
Copenhagen, CERN (European Couneil for Nuclear Research). 
() F. J. Dyson: Phys. Rev., 75, 486, 1736 (1949). 
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toward the treatment of problems to be dealt with later, such as renorma- 
lization, investigation of the convergence of the perturbation expansion, search 
for a non-perturbative approach, etc. Our main purpose, here, is to intro- 
duce appropriate algebraic tools for dealing with questions of this sort; the 
fact that we consider only electrodynamics is not a limitation in this respect 
because, as we shall see, the extension to other types of interactions requires 
only trivial changes from the developments given here. This is borne out 
from the fact that our results depend, in their essential part, only on the 
commutation properties of the fields, so that the electrodynamical case, which 
requires the consideration of fields of both statistics, suffices to give all the 
information needed in more general instances. 

The Feynman graph technique, which has marked a tremendous impro- 
vement upon previous computational devices, suffers nevertheless from several, 
serious shortcomings. It can be used in practice only for the calculation of 
the very first terms of the power series expansion of an element of the S-matrix, 
and, even so, only if the number of particles created and destroyed is very 
small: this last feature forbids its application to problems involving, say, 
mesons and complex nuclei. More serious than this, from a conceptual stand- 
point, is the fact that the contribution of a single graph is not a gauge-invariant 
expression: this shows that the intuitive interpretation usually given to graphs 
is rather objectionable, and leads one furthermore to suspect that the standard 
renormalization procedures, in which each graph is treated apart and the 
different contributions are summed together after renormalization, may be 
unnecessarily complicated. Worst of all, from our point of view, is, finally, 
the fact that the graph technique does not lead to compact and formally 
simple expressions for the terms of the power series expansion; had one such 
expressions, one might hope to sum the series, at least in some simple cases, 
so to get a result which, although derived heuristically from a perturbation 
treatment, might well have a more general significance. 

The present considerations, which are intended to eliminate these draw- 
backs, are bound to remain of a merely formal interest until it will be possible 
to deal satisfactorily with the other questions mentioned before. Well aware 
as we are of the difficulties of the task still ahead, we feel, however, that this 
would be an impossible one, were not first clear the formal structure of the 
S-matrix. 


4:2. - Our investigation stems directly from the well known works : by 
Dyson (1) and Wick (?) and from a work (*) in which it was shown that the 


(2) G. C. Wick: Phys. Rev., 80, 268 (1950). 
. (8) E. R. CAIANIELLO and 8. FUBINI: Nuovo Cimento, 9, 1218 (1952) referred to 
in the text as A. ; 
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spur of the product of any number of Dirac matrices can be systematically 
reduced to the evaluation of a pfaffian. We show, by suitably generalizing 
the algebraic theorems of A, that the analysis of Wick’s chronologically or- 
dered products T can be improved so to yield a concise and symmetric ex- 
pression for any term of the power series expansion of an element of the Dyson 
S-matrix, irrespective of the number of particles and of the order of the term, 
gauge-invariance being self-evident: one, indeed, which can be written down 
at once, as soon as the process and the desired order are specified, by means 
of extremely simple mnemonic rules, without having at all to go through the 
rather elaborate algebra of Wick’s products. We shall have to use, for this, 
a new algorithm, which is defined in sect. 2°5. 

It may be said, at this point, that the algebraic results obtained thus far 
by us extend rather beyond those reported here. We have, for instance, 
proved Furry’s theorem, quite in gereral, by direct algebraic calculation, as 
a check on our formalism; also, we have reduced our final formula (34) to yield 
only the contributions of all inequivalent Feynman graphs, each only once 
and with the appropriate sign and weight. We forego, however, entering 
here into such details, due to their limited interest and to the quite cumber- 
some developments necessary for their derivation. 

Part 2 deals exclusively with such algebra as is strictly needed for our 
present purposes; the problem is stated and solved for the electromagnetic 
contribution in Part 3; the electron-positron field is studied in Part 4, and 
there we also give the final results. Finally, we extend in the Appendix to 
U(t, t) the results obtained in the text for S(co). 


2. — Algebraic Preliminaries. 


2:1. — It is most important, in a work of this sort, to pay particular attention 
to the notation; this should render writing as simple as possible, while exhi- 
biting the relevant features of the expressions to be considered. Thus, we 
shall encounter throughout expressions such as 14[S*(a, — de) |p Bs ete.: the 
only thing which matters, for us, is the behaviour of this quantity with respect 
to permutations of the indices h and k: we denote it, briefly, with (hk); and so 
on. Likewise, determinants, pfaffians and the other quantities introduced below 
are indicated, whenever convenient, with an abbreviated notation. By doing 
so we achieve also the advantage of giving our formulae a greater generality 
than would be surmised from the special problem to which we restrict here 
ourselves: the theory of the Dyson S-matrix, which we consider only in the 
absence of external fields, applies as well if such fields are present; one has 
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only to read (hk), say, as Yo[ S*(a,, 2)y"*]a,.,, to get the corresponding formula 
from our results. Also, as is shown in the Appendix, the formula for U (t, to) 
is again essentially the same. 

In the remaining sections of this Part we define the notation and the alge- 
braic quantities which occur in our discussion. 


2-9. Determinants. — We denote the determinant whose elements are the 
symbols (hk) with 


(hıkı) (like) «.. (Paka) | 
MO ehe (hoky) (ok) + (hzkn) | 

= fe Ta 2) pie A 
| (ink) (ik) … (Anka) 


This is Sylvester’s classical notation, for (hk) =a ,. No ambiguity may 
occur, although in general (hk) £ (kh), since the row index of an element, 
which comes from the first row of our symbol, must always precede the 
column index (from the second row). 

The antisymmetrized product of N, wave functions u™, w®,..., u» (omitting 
the indices which specify their components) is denoted with 


ud.) uM(x) ... UD (xx) 
5 wm un. an — | WC) ur). u) | 
h Ope ame eh Ni jr | 
| u (%,) Uo) (Lo) aes UV (Ly,) | 


(2) is a particular case of (1), when (uk) = u®(x,). Similar devices are used 
whenever convenient. 


2-3. Permanents. — A «permanent » (4) is quite similar to a determinant, 
the only difference being that each terms of its development is taken with 
the + sign. Only some of the theorems of determinant theory, of course, 
are valid for permanents; some either do not hold at all (such as that on the 
product) or require obvious modifications. We cannot enter here into such 
questions; the only permanent of interest to us is that defining the symme- 
trized product of P, wave functions; we denote it, in analogy with (2), with. 
the symbol 


| 2 (a) dm) … (Cr) | 
2D 22) … Po | 22(%,) 2 (a2) … 2 (x) 
3 — | | 
(9) NORME oP RE AT VRR Ce 
PR) Por) ... 2F(Hp,) | + 


(4) D. Lirrtewoop: The Theory of Group Characters (Oxford, 1950). 
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Here too, one may introduce the simbol [2%%] = 2%(x;) and obtain (3) as a 
particular case of the analogue of (1), which we forego writing. Square brackets 
are used, consistently, for elements of symbols, and the symbols themselves, 
which are expanded with à « + sign » rule. Round brackets (eventually with 
a prime) otherwise. 


2:4. Pfaffians — The pfaffian algorithm, well known since a long time, 
has been applied in A to questions of quantum physics. Its appearence, of 
course, is not accidental, but stems from its close connection with the sym- 
plectic group, which characterizes quantities obeying anticommutation rules. 
This remark, indeed, has been the starting point of our present considerations, 
as it provided the « missing link » in the determination of the symmetry pro- 
perties of the algebra of ordered products of operators. 

A pfaffian is defined by a symbol, or a triangular array, which is expanded 
according to the rule: 


(4) OA ZO) = 


= a NES ay et) een 
(2n _ 2 2n) | 


>’ means here, and throughout, in the following, summation over all the per- 
mutations À, t,..., %n Of 1, 2,..., 2n. which satisfy the limitations 


di de bgt ove ton < tony 


Vie (TA Ve PRE D ‘ 


A pfaffian (in analogy with a determinant) can be expanded by the elements 
of one if its lines: «line h» is that, in the triangular array, which contains 
all the elements carrying the index h, regardless of whether h is the first or 
the second index. The element (hk) belongs to the lines h and k; crossing 
these lines out, what is left is a new pfaffian, the «minor » of (hk). The rule 
for the expansion is then: « A pfaffian is equal to the sum of the products 
of each element (hk) of a prefixed line by its minor, each term in the sum 
being given the sign (—1)****1», Many other remarkable properties cha- 
racterize pfaffians, on which we do not insist here. 
It is well, to avoid possible misunderstandings, to remark explicity that 
the symbols used here to denote the elements of a pfaffian have nothing to 
do with the particular definition given them in A: for us, indeed, (hk)+ (kh), 
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while in A we had (hk) = (kh). However, the fundamental relation (A.12) 
(1 2 ... 2n)? = det. || nel 


remains true here too, when the elements of the (antisymmetric) determinant 
det.||a,,|| are defined by means of the relations 


(3) Op = 0 for shi =k: as (hk) for hk; au = — (kh): for h>% 
instead of (A.10). 


2:5. A new algorithm. — To complete the set of algebraic tools needed in this 
work we have to introduce an algorithm of which we were unable to find tra- 
ce in the mathematical literature. It is related to pfaffians as permanents 
are to determinants, and appears to be the «natural » computational tool in 
the algebra of operators obeying commutation laws. We use for it, in ana- 
logy with the symbolism introduced in sect.s 3 and 4, the notation: 


ey a ss ee 
23] … C 


(6) [1 2°... 2n] = a ; ee >” [ariel is da] -.- [énrien]- 


[on —1, 2n] + 


In particular, what was said on the expansion of pfaffians by the elements 
of a line holds also in this case, save for the sign which is always +. We 
need a name for it: we call it «hafnian », a choice which seems convenient 
for its assonance with the name of the similar pfaffian algorithm, and is sug- 
gested by the memory of our pleasant association with the place where its 
consideration first occurred to us. 

It would be quite interesting, and probably useful, to have a complete 
and systematic treatment of the relations among pfaffians, hafnians, deter- 
minants and permanents. This lies beyond our present scope; we confine 
ourselves, here, to stating without proof two theorems which emphasize both 
the analogies and the differences between pfaffians and hafnians: 

«A symmetric permanent of even order with vanishing elements in the 
leading diagonal (the analogue, that is, of an antisymmetric determinant) is 
the exact square of the hafnian formed by the elements lying above the leading 
diagonal ». 

«[1, 2,..., 2n] = P[1, 2,..., 2n], where P denotes any permutation of the 
indices 1, 2,..., 2n, provided [hk] = [kh] ». 


106 - Il Nuovo Cimento. 
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3. — Caleulation of an Element of the S-Matrix: Electromagnetie Field. 


3:1. — Since we are concerned exclusively with algebraic developments, we 
avoid any explicit consideration of such questions as the introduction of con- 
vergence factors, limiting processes and the like; this simplifies considerably 
our formulae, while no trouble ensues for the reader who wishes to modify 
the final results so to take these points into account. We work throughout 
in configuration space and use the interaction representation; we need not, 
and do not, restrict the initial and final particle states to plane waves. 

We are interested in the element M,, of the S-matrix which describes the 
process in which 


n electrons, m positrons, a photons 


are destroyed and 


p electrons, q positrons, 0b photons 


are created, in states specified by some prescribed wave functions; on these 
we make, for simplicity sake, the hypothesis that all states are orthogonal 
to one another. Clearly, n-+q=m-+p; we call N, this number, and set 
P, = a+b. The structure of our formulae depends essentially only upon N, 
and P,. If X, and Y, are the initial and final state vectors of the system, 
we have (4 = c= 1): 


où | 
= 


N 


en | NES | 
Mi Sail 2 | Al TT Fe avatar) B\ 
N Ö | U arr 


(the sign of ¢ is of course irrelevant). In (7), the integration and summation 
symbols operate upon all written variables and indices; 7 denotes Wick’s chrono- 
logically ordered product. We use Wick’s theorem and notations for 7 and 
S products. Since, in our case, the operators of the boson field commute 
freely with those of the fermion field and contractions between fields of dif- 
ferent statistics vanish identically, we can split our problem into two separate 
problems, one for each type of field (this justifies the remark made in the Intro- 
duction on the immediate generalizability of our results to any type of inter- 
action, provided use is made of the interaction representation). This requires 
an obvious decomposition of Y, and Y, into the Kronecker products of vectors 
which describe separately the fermion and the boson part of the field (denoted, 
respectively with, the labels f and b). 
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3:2. Electromagnetic field. — Let: 


| Pr = ag... afar |O>, 
(8) | Ye» = * * @* |O\ 
Ds Ay... Ay Ay |O>. 

(8) describes the initial and final situations, with respectively a and b photons, 
in states 1, 2)... avand 17, 27;:, 0, created by the operators a* acting on the 
vacuum state|0>. Writing will be simplified by using the same symbols 2 
to denote both destroyed and created photons, leaving to their labels the task 
of specifying the case of each particle. Thus a (Pe) will denote the wave 
function of a photon of state specified by the index h, which is destroyed 
at point æ,, if À <a; or the c.c. of the wave function of a photon created at +, 
ifa<h<P,. Call, for short, JJ the product A, (@)A,, (2) … Ay, (Ly) 5 Wick’s 
theorem states that 


(9) TP le SW Le a Sho HLS ae ea 


where JJ®” denotes any of the values /7 assumes when A pairs of operators 
are contracted, and > //?” denotes summation over all such values, for 
given h. Following WICK, the colons denote that all «free » operators (i.e. those 
which are not contracted) are to be split into their destruction and creation 
parts and then the products of these written in normal form (i.e. with all 
creation operators at the left). It is then immediately seen that, as is well 
known, the only term of the expansion (9) which contributes to the boson 
part of M is >: 11° :, having set N= P,+2Q: 


(10) Po | TD) PO) = > HY] TT? || PP». 


(M® 0 only if N and P, have the same parity). Let i denote the operator 

3 N 
A,@); consider the ig Q 
separating 2Q operators h,, Mry.., Nye (hy <h <...<h,,) from the remaining 
P, operators ji, jes) Jp: then, clearly : 


| combinations C,,, each of which is obtained by 


(11) (P| TUT) |P> = (D [ht] --- real | : fide ee Ir: | PP?» 
C2Q 


where [hk] is the value of the contraction between the operators A,,,(&n) and 
A uy (©) (chronological ordering): 


(12) [RE] = YoOpnugD? (tn — an) ; 


>’ has been defined in Part 2, any toy.) tz. are permutations of the indices 
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hy, h,.., h,, fulfilling the limitations implied in the definition of >’. Intro- 
ducing the hafnians defined in 2°5 we can write (12) as: 


(13) CPE | TD) | PP) = & [rhe … Nag PO |: fide --- fr, | PP) - 
C2Q S 


Clearly, the last expectation values reduce to the permanent 


DE ERO sq CU) 


yh Je + aes 


of the wave functions of the initial and final photons, which appear thus duly 
antisymmetrized: a well known consequence of the formalism. Also their 
indices are correctly specified, by defining the symbols [zk] as 


[eh] [Hee] = lm). 
Finally, we prove a theorem which permits to re-write (13) as a single 
hafnian, exhibiting thus at once all its symmetry properties as a consequence 


of the last theorem stated in Part 2: 


Theorem: «The expression (13) is the expansion of the hafnian of order 
N+P, (even!): 


(14) 12. Nee... 200] = 


[12][1 3]... [1 N] ER are. El zo 
[23]... [2 N] LES 29) 2 Se] a See Zee 
(MIEN Ne Lee oe BN Dal 
[N gh] [N ON or [N Po) | 
= 0 . . . . 0 | 
é à : 


with’ [22] 10), 


Proof. This statement can be proved in several ways, the simplest of 
which is perhaps to proceed in two steps: 


* a) to verify that each term of (14) coincides with a term of (13), no two 
terms beeing equal; this is apparent on inspection, having regard to the 
structure of (14) and to the rules given in sect. 2°5 for the expansion of a 
hafnian by the minors of one of its lines; 
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b) to verify that the expansion of (14) yields exactly as many terms 
as the expansion of (13), the last being, clearly, in number of 


N! 


(N — P;)!!° 


N 
P,! (2Q— 1) !! = 
(39) Pal @@—0 
Let H,,, denote the number of terms of a hafnian such as (14), with structure 
determined by the values of N, P,. On expanding it by the elements of its 
first line one easily finds, from inspection of the structures of their minors, 
the equation 


H yp, == (Ne 1)Hy_2,p, aie Poy 3,7, 5% ’ 


which is readily solved by taking V,, = N!/(N — P,)!!! q.e.d. 

The connection with the graph techique becomes apparent by observing 
that the terms of the expansion of (14) correspond to all, and only, the con- 
tributions one would obtain by drawing in all possible manners, for N fixed 
vertices, Q internal and P, external photon lines. 


4. — (Continued): Electron-Positron Field. 


4:1. — Denote the product of the operators of the fermion field with 
(15) P = 123... 2N= Wa,(1)Ve, (1) … Vay(@y) 5 


where h = 2(h)--1 stands for the operator Po )(®»), À = 2(h) for WBqy(Liany ) 
((h) =1, 2,..., N). We denote with (x) the maximum integer contained in 
(«+1)/2 (usual notation: [(«+1)/2]: («) = [a+ 1% — (—1))]. 

Wick’s theorem states that 


(16) TD) CPS ED ee 
where the symbols, have the same meaning as in Part 2. Let 


15 Woy = b#..b#at...ar|0), 
1 PO — b* ... DKA... a [0 , 
denote the state vectors which describe the initial and final situation of the 
fermion field; a* (b*) are the operators which create electrons (positrons) in 


the states specified by their labels. 
If one is interested in the calculation of M,, for a given process, its phase 
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is quite irrelevant; this is not the case if one requires the knowledge of all 
the matrix elements of S, say to evaluate some f(S), and it happens (as it 
does) that the phase of M,, is not simply a difference o(1) — o(F). We adopt, 
therefore, an ordering of all the Y such as typified by (17) and determine the 
phase of M,, with respect to it. 

Proceeding exactly as in Part 8, one sees immediately that the contri- 
bution of the fermion field is given by 


(ES) IL Vi)‘. (anatase) PP | ihe Ta: PP , 


where: N—N,+M; (hk)' denotes the contraction between the operators h 
and k (we use, below, primes also for pfaffian and determinant symbols: 
this only means that the elements of their arrays are primed symbols; we 
do so to simplify things in the final results and to avoid the introduction of 
a too eterogeneous notation): p is the parity of the permutation 


Vide ce. Laide ++ jew, 
with respect to the natural order 1, 2,..., 2N. One sees easily that p reduces 


to the parity p, of 7,7, ... 7 with respect to the natural order h, < hy < «.. < hey 
of them in 1, 2,... 2N. Indeed, 


p = Pi + P; 35 Di; , 
where p;, parity Of 12... Jey, With respect to their natural ordering, is = 0 


since this ordering is not changed, and p,,— number of inversions of hyhy...hoy 
with jie Jon, is given by: 


| au 2M +1 
Dis na = 1) 2) nom > | Le i. 


But (see below) the only terms 40 of (18) come from combinations having M 
even and M odd indices h: it follows that p,, is even: q.e.d. Formula (18) 
becomes, thus: 


(19) I (ple han) CPD |: jada eee fan, © ES 


Com 


This formula, which is the exact counterpart of (13), shows clearly the con- 
nections between hafnians and Bose statistics, pfaffians and Fermi statistics. 
We proceed now to the most tedious part of our task, the further simplifi- 
cation of (19); this requires that we study separately the pfaffians and the 
expectation value of normal products which appear in each term of (19). 
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4-2. Evaluation of (hyhy... hoy)’. — It is easily seen that the contraction 
symbols assume, for chronologically ordered products, the values: 


{ (h for h, k both even or both odd 
(20) (ol We= en xq) (EE — Li) for h odd, k even 


17 à > 5 
/e SF a (Du — Lu) for h even, k odd, 


with 8? = (OT m)Ar. 
The first of the relations (20) permits to reduce (hh, ... hes)’ to a much 
simpler form. The only terms #0 in its expansion are those in which each 


(hk)’ contains one odd and one even index. We can thus split hyh, ... han 
into two sets: 


kıky...k,, containing all the odd indices, 


hih,...hi,, containing all the even indices; 


we label them according to their natural sequence in Aıhz ... hex, So that 
bchi<..<h,, m<lm<..<h,. Remembering the observations at the 
end of sect. 2°4 (formula (5)) and formula (14) of A, we set 


(i = forsehr — 
Gy, = (bh) . for h>k, 
aux = — (hk) for h<k, 


so to have: 


(21) (hho ee ho) = 


1)?Q, ; Uy g, a+. Ui eh 


(2M) !! 2 ur into glg" gm-ıtam 


the summation being extended to all permutations of the indices (of parity p). 
(21) contains 2*(M!}? terms + 0. We prescribe now that dt, ... demi be 
given by ie bile ae Le in their natural order, so that 4, üay..., ty are à per- 
mutation of h,, here hy: This eliminates again the factor 1/(2M)!! and yields: 


(22) (aka de Ron) = D Lares, nhs, + 


2 Ryan ? 


the sum is BO over all the M! permutations of hih,... h, and p is the parity 


of kiki uw with respect to Ihe «Roy. Call 7’ the parity. with respect 
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u 


to Iyhy...hoy of kik,...k hihi... hy (that is, with all odd indices shifted to the 
left); if p’ is the parity of iss... ig, with respect to jh, … h,, p is given by 


, 


M~1 M ; 
p=r+ Sas +p'=r+(p) +2 ; 


S=1 
so that 
ee ee 
(23) he. hax)’ = (— ie} | a oe az | 
Gr n Un ay, nu 


The same result follows, of course, from formula (A.17), after a suitable gene- 
ralization. 


4:3. Evaluation of PP: jude jy: | Pr 2. - Clearly, there are N, odd and 
N, even indices j. To analyse this expression we write 


2e j intl 
(24) T= 2 jiJo + Ian, = Lu ” (Geigy) + yt CN) : 
where | 
UT = ANT = U) destroyes (creates) electrons, 
or == y (ù = 0) » » positrons. 


We are only interested in the products contained in expansion of J which 
destroy m positrons, destroy n electrons, create p electrons, create q positrons. 
To get the sign straight, we permute the factors of J so to arrive at the 
«standard ordering » 


LE AE WEE HWE Be RE RE nee ge 


p''a+1 "a+2* Ny ?+1 P+2° 


where 


N ss . 
he hy, are the even j in their natural sequence 
4 


al ta 


gic al ee OAA PD » » 

Thus, J = (—1)7.J', p,, being the parity of the permutation which leads 
from J to J’. The advantage of introducing J’ lies in the fact that the 
expansion of <W|J'|W> contains a term which, as is seen on inspection 
from (24) and the definition of J’, has certainly a positive sign, since no anti- 
commutations intervene in its evaluation. Denoting with w® (vo) the wave 
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function of the electron (positron) in state k, this term is: 


ta’) =) =D") ir) 
By! a) Paty’ no) ma) us) A Gp) 
Cm) @) (n) (1) 
mu) (20) N (543)) Peat) an) Pent aa (hia) 


There is no trouble then in seeing that 


We? oe UE MD vom! fo GED UM"! 
(f) (MN fs Der 4 .. 
cle cn Le 
CR NL RS el ele 


gd 


where we use the abbreviated notation defined in sect. 2:2, and (uA) = 
= up), etc, in an obvious manner. 

To evaluate p,,, call t” the parity with respect t0 jrjz + jay, OÙ Kee Kuh Mn, 
(that is, with all odd indices shifted to the left, in their natural sequence); 
this ordering becomes that of J’ after pq-+mN, transpositions, so that 
Por = T'+pq+mN. To get a more compact notation, we write 


pt Sey OED ye ee RO Nr LEE ORY a UP) = wo, 


so that, taking due care of the signs, (25) becomes: 


yy … 9 a0)\" UPy® … yy ro? ! 
26 pu) J wo) Se (tt ee u u 4 u 1 ! 
eo Er | | mn ( ) ie Toe a's ok | m Ve cet Ney )» 


where 0 = pq+mN,.+ (2)+ra+ (6 )+ mp; k" (h") are the odd (even) indices in 
the order in which they appear in J, and 


are 


27 1 " ! 
a | ara an) 


44, — Collecting our results we get, from (23) and (26): 
(28) PS |T(P)| Pr) = 
| Ayn! orale Oe Di yy | 
oF A “fi ay LE z 1“ ad =). Dao) ur) es uNo 
= (— 1) e Di i. |: 5 6 = ° Fe LS ki he ae hn 2 
; yh Uys | 


: 2 
The sum is over all the or combinations of 1, 2,.... 2N such that all h’s 
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are odd, all h’s even; t’+7" is the paritiy of the permutation 
(29) RE oR la rar Wh ah ehr 


(see discussion after formula (8): p,;; = 0). 
Consider now the determinant of order N+ Ny 


| day On we Ay ay ro) (1 vn)! 
| 
| @32 M34 se As on (STD) (3 yo)! 
| 
A = |Gev-1,2 Gew-1,4 …..  Ggy—s,2n (2N— 1,0%)" …  (2N— 1, vo)! 
| 
| (w2)’ (uma) … (u®, 2N) 0 7 0 
(UD) NEA) RE or DIN) 0 ape 0 


It is readily seen that, to within the sign, each term of A coincides with a 
term of (38) and that the total number of terms is the same in both cases. 
Expanding A by the minors of its last N, rows, and then the adjoint of each 
such minor by the minors of its last N, columns, one gets a sum of products 
of determinants like (28), save that now the sign of the running term is given 
by (—1)? with 


Tr oa 1 alt 1 
(ND) OSE) Remy eer 
4i=1 
os, 3 12%, 
3 2 (he nahe ee 


On the other hand, (29) is obtained from the natural ordering 1, 2,..., 2N 
through the following steps: 


@) chan TOME aes macaw GOm el mn sote IN I DNA, 
; (7) 
change = : 


b) change the latter into k,..k,k,..k,h.…k,h.…h,: parity change gi- 
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c) change into (29): parity change = MN). Therefore: 


T'+ a" Io) =, 


M N M 1 
=o+(5)+( \+5 + mt; 


and finally: 


(30) CO! T(P)|PP> = (— 1A, 


with 
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It will be noticed that this gives also the sign for the inverse process (matrix 


element M,,), a8 is seen after a trivial computation; also, 


that py, depends 


upon n, M, p, q in such a manner that is not amenable to the form (1) —œ(F). 


48. — Let us now give A its explicit form. 
the definition of a,., we find 


| 1)’ VS, (Le — 21) 
Ya) (22) 


ws” (1) wp, (@2) 


OF 
VY, Spyx,(@y— #1) ve (a) 


7 uf 
14 SByx, (Ex — 82) ve (a2) 


(NN) 
up} (Ly) 


(N, 
uby (Hy) 


ve (Ly) vee 


0 


It is at this point convenient to incorporate also the matrices y£g 
sums over the indices « into the expression of A. We take, for convenience 


of notation, the transpose of A instead of A. The result is, 


notation: 


Ay, bayer, AN 


Ss 12... N: VOVo-... 
60) ne re a 


TO). san DE ENDE oor 


From (20) and (27), remembering 


0 


and the 


in abbreviated 


}: 
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(hk) = 12 [87 (x, = Xx) ye |B, B, ’ 
(wh) = (hu®) = uP (a) 
(33) 7 (ie) 7(k) - (k) 
(VR) = (RV) => oF (en)yen 


i 
œB. 
= RPh 
nl 


(Vu) = (UV) ==). 


The Heisenberg prescription simply amounts, in this formalism, to setting 
(hh)'= (hh) = 0. If there are external fields, however, the symbols (hk) have 
to be re-defined in an obvious manner, and (hh)+0. Formula (32), and so 
formula (14), still hold. 

The connection of formula (32) with the graph technique is again evident: 
its expansion contains all, and only, the contribution one obtains by drawing 
in all possible manners, for N fixed vertices, M internal and 2N, external 
fermion lines. The antisymmetrization of the wave functions of the initial 
and final states is apparent. 


46. — Finally, from (7), (14) and (32): 


en PP ose NAO ST) 
84) M,,=(—1)" — | |: 125. Neo... 220) 
Fe ee ne DIS EN T0) | $ ] 


/ 


From the structure of (34) we see that the determinant and the hafnian cor- 
responding to the term of order N are simply the principal minors of orders, 
respectively, N+N, and N+P, of a determinant and a hafnian of infinite 
order. (The determinant is a well known one, in the Fredholm theory: this 
is more than a coincidence; we shall return on the question in the continuation 
to this work). The lowest orders may correspond to processes described by 
«non connected » Feynman graphs, which must, of course, be considered on 
the same footing as all other terms. 
In particular, the vacuum-vacuum transition gives: 


al LOTS NAN 
35 Mi ae DIN] 
2) = 22% = x) ] 


We expect that M,, and similar simple expressions, will play an important 
role in the investigation of deeper questions. Any M,,, in particular, should 
derive from (35) (written for the case of external sources) by means of suitable 
variational differentiations. 

Clearly, (34) contains, for each- order, many identical terms, as is with 
graphs. As we have said in the Introduction, we forego reporting the formulae 
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which accomplish its ultimate reduction, because much of its simplicity would 
be lost, and furthermore one would fall back into the old troubles. We just 
record here a trivial simplification of (34), which is symmetric enough to de- 
serve perhaps mention: 


ex 
36 M,, = (—1)?r ee ave TS 
{ ) FI ( ) eee (N Ce, P,) 1! 


TE PMU EN CT = 
aba Pern ite \21341 120 — 120) 2h (fat) ala: 


It is clear that formula (34) would yield, on inspection, the exact number 
of Feynman graphs pertaining to any wanted process (involving orthogonal 
states) and order — a number whose approximate evaluation has been the 
object of several investigations. We leave this to the reader. 


In conclusion, the author takes pleasure in expressing bis sincerest thanks 
to Prof. P. KRISTENSEN for many interesting and stimulating discussions. 


APPENDIX 


Formula (34) contains integrals ranging over all space and time values. 
It is just as easy, however, to obtain U(t, t,) instead of S = U(+00, — co): 
just take t,t, rather than Ho as limits in the integrations. It is instructive, 
however, to abandon the use of chronologically ordered products, and therefore 
define the contraction symbols as 


E.M. Field: 
(i) [hk | = du Pa (Ln 75 Lx) 2 
El. Pos. Field: 
— 18 — Tin) for h odd, k even 
(ii) hee | à (ny \% (k) BE ’ ’ 
| — Sr (La — Loo) for h even, k odd. 


so that (34) becomes: 
(iii) M7zı(t, to) = (1)Prr- 


t ty A 


ID NOT 2 
: ex | dx, | da... | dx | Ju NND, eo 
wollen | | | = > 11D pa, MONDE un : 


40 tn to 
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where the only changes with respect to the previous notation are given by (i) and 


ST (ar — 2r)y "x |p, 6, for h<k, 


oY ee apps for h>k. 


For t= co (iii) reduces, of course, to (34), as is seen from the well known 
relations: 


Ay RA) foray Or 
A(x) = — 2i47(x) for 4%< 0. 
RIASSUNTO 


Si da una espressione compatta per il termine generico dello sviluppo perturbativo 
della matrice S di Dyson, per il processo di scattering più generale. Tale sviluppo pos- 
siede un alto grado di simmetria, del tutto nascosto dalla tecnica usuale basata sull’uso 
dei diagrammi di Feynman, che presenta anche vari inconvenienti da cui il nostro 
metodo & esente. L’espressione trovata per S indica l’esistenza, nella teoria, di più ri- 
poste proprieta; queste saranno studiate nel seguito a questo lavoro. 
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Summary. — A comparison has been made between the theory of electron 
pair production by photons and experimental data, for energies up to 
400 MeV. The following points have been analysed: 1) Energy par- 
tition between the two members of the pair; 2) Statistical distribution 
of the momentum transfer to the pair; 3) Distribution of the angle. of 
divergence of the pair. 


1. — Introduction. 


In an investigation, with nuclear emulsions, of the soft component of cosmic 
rays at 2860 meters above sea level ('), enough data have been collected to 
enable us to compare the theory of pair production (?) with the experimental 
results at rather high energies. 

(4) G. BARONI, G. Cortini, A. MILONE, L. Scarsı and G. VANDERHAEGHE: Nuovo 


Cimento, 9, 867 (1952) Cited below as B. 
(2) W. HEITLER: The quantum theory of radiation (Oxford, 1944). 
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Verifications of the theory at relatively low energies (1-20 MeV) have been 
made by many authors and although some points need more discussion, one 
can say that the theoretical results are more or less in agreement with the 
experiments. Electron pairs by photons up to 70 MeV have been studied with 
nuclear emulsions by MCDIARMID (*); while DE WIRE and BEACH (*) investi- 
gated, by the pair spectrometer method, the energy distribution of positrons 
produced by 270 MeV photons. 

The comparison will be analysed in the following points: 


a) Energy partition between the two members of the pair. 
b) Statistical distribution of the momentum transfer to the pairs. 


c) Distribution of the angle of divergence of the pair. 


As in our case we have no indication of the direction of the generating 
photons, it is impossible to study the momentum transfer to the nucleus or 
the angle between either electron and the photon. It is possible however to 
make the comparison using a) and b) because the knowledge of the direction 
of the incident photons is not needed in this case. 


2. — Energy distribution between the two electrons of the pairs. 


We measured the energy distribution of the electrons from 196 pairs ob- 
served in our plates. As we are not able to recognise the sign of the charge 
of the electrons, we are indicating with H, the energy of the electron of lower 
energy and with #, the energy of the other electron. The statistical distri- 
bution of the pairs has thus been represented as a function of the ratio 
E,/(E, + E,), that ranges between 0 and 0.5. 

The pairs were divided into two groups: 


a) Pairs with total energy between 5 and 100 MeV. 


b) Pairs with total energy larger than 100 MeV. 


The distribution of the pairs as a function of the total energy is given in 
Table II, reference B; the pairs with an energy larger the 400 MeV are only 
three. In Fig. 2 we give the hystograms of the two groups a) and b) of pairs; 
the errors marked are only the statistical standard errors. Fig. 1 shows gra- 


() I. B. McDiarmip: Can. Journ. Phys., 30, 670 (1952). 
(*) J. W. De Wire and L. Bracu: Phys. Rev., 83, 476 (1951). 
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phically the theoretical curves of the cross-section for the production of po- 
sitrons of energy Æ. by photons of energy k= (EH, + E), evaluated with 
the formula: 


6 6 2,610 cm? 
6 / 
(1) o(E;) = o(E_) = a (HY. + EL )fily)— 5 in 2+ 


where o, = 2?r}/137 = 2.6-10-?5 cm? for Ilford G.5 
photographic emulsion; y = 100 me*/(E, E_ Z"*); 
fi(y) and f,(y) are functions known only numerically. 
Formula (1) was obtained by BETHE and HEIT- 
LER (5) in the Born approximation, taking into 
ae count; the screening au r ; Fig. 1. — Theoretical energy 

Each curve of Fig. 1 was obtained by taking distributions between the 
the average curves for the different elements two electrons of the pair. 
which constitute the photographic emulsion. 

We plotted the solid lines of Fig. 2 from (1) taking into account the known 
energy distribution of the pairs of each group. The hystograms and the theo- 


o OY 02 03 04 05 E,-mc° 


= 2 
hy -2mc 


30 
20 : 
Way 
Ss aS ee ı à 
ie) 0,1 0,2 0,3 0,4 0,5 ie} 0,7 0,2 0,3 0,4 0,5 
Es = mac” Ey = Moc? 
hv —2m oc? hy —2m,c* 
Fig. 2. — Energy distributions between the two electrons of the pairs. 


a) 5MeV <(E, + Ey) <100 MeV; 6) (E, + H,) > 100 MeV. 


retical curves of Fig. 2 are normalised at the same area. We may point out, 
finally, that the errors due to the use of Born approximation, are negligible 
in connection with the energies under study. It can be said that the agreement 
between the experimental distribution is rather satisfactory for the group b) 
of the pairs; for the group a) the deviations are clearly larger than the stati- 
stical errors, for the pairs with the value E,/k less than 0.3 were fewer than 
expected. 

A reasonable interpretation for this discrepancy shown by the group with 


() H. BerHe and W. HEITLER: Proc. Roy. Soc., A 146, 83 (1934). 
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total energy lower than 100 MeV, is that the pairs with large asimmetry are 
difficult to observe when one of the electrons has à very small energy. In 
fact in this case the effect of scattering reduces strongly the probability of 
finding out the pair: 1) because the scattering pushes rapidly one track away 
from the other, making ineffective the usual criterium of recognising a pair, 
based on the convergence of two near tracks; 2) because the slow electron 
has a great probability to leave the emulsion after a short path and conse- 
quently even if we can recognise the pair, we cannot utilise it as it is impos- 
sible to measure the energy of the shorter track with an acceptable precision. 
It is difficult, specially for the first of the two reasons of loss, to make a pre- 
cise correction. 

The drop below 30 MeV in the photonic spectrum (B, Fig. 4) can be explained 
in a similar way; the loss was evaluated as at least 20%. This reasoning be- 
comes plausible when comparing the number of the pairs detected by accurate 
search with that by a quick one (columns 2 and 3 of Table II in B). 

A similar drop was pointed out by MCDIARMID () in the study of pairs 
produced in plates by a synchrotron spectrum of 70 MeV. In this case too 
large deficiency was noticed of pairs with an energy below 30 MeV as well 
as of pairs with large asimmetry. This discrepancy still remains after the 
correction, made by MCDIARMID according to cause (2) of loss resulting from 
the scattering; however we cannot evaluate the importance of the effect (1). 

We mention the fact that the measures by KocH and CARTER (°), with a 
cloud chamber for energies below 20 MeV, have shown, in the case of asym- 
metric pairs, a discrepancy in the same direction, but of lower amount, for 
in a cloud chamber the above mentioned causes of loss for asymmetric pairs 
are practically ineffectual. 

The absence of disagreement for pairs with energy larger than 100 MeV 
substantiates the arguments given above (in fact the effect of the scattering 
on the tracks at these energies is much lower) and at the same time is in 
agreement with the results of DE WIRE and BEACH, obtained with photons 
of 270 MeV. 


3. — Distribution of the momentum transfer to the pairs. 


A point of interest in the theory of pair production is the momentum q 
transferred to the nucleus. The distribution of the modulus q and of the 
angle between q and the direction of the incident photon was studied in detail, 
with a cloud chamber, by GROSHEY (?) at the energy 2.6 MeV and by MODESITT 
and Koch (8) at energies up to 20 MeV. The experimental results of the latter 


(°) H. W. Kocx and R. E. CARTER: Phys. Rev., 77, 165 (1950). 
(?) L. V. GrosuEv: Compt. Rend. USSR, 26, 424 (1940). 
(8) G. E. Moprsirr and H. W. Kocu: Phys: Rev., 77, 175 (1950). 
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authors, compared with the theoretical ones obtained by JOST, LUTTINGER 
and SLOTNICK (°) in the Born approximation without screening, show a dis- 
agreement for large momentum transfers and for small angles. However these 
experimental results are probably affected by systematic causes of errors (*°), 
hence the significance of the disagreement must be ulteriorly clarified. 

As we mentioned in the introduction, in our case we do not know the 


Arbitrary units 


(0) 1 2 3 4 > 6 


Fig. 3. — Q distribution of the pairs a) with screening; b) without screening. 


direction of the incident photon and for this reason it is not possible from 
the measure of the total momentum of the pair P = p, + p, to obtain the 
momentum q = (k— p) transferred to the nucleus. Therefore we shall discuss 
directly the distribution of P, or even better, as it is more convenient, the 
distribution of the related quantity Q = (k?— P?)", which represents the 
energy of the pair in the center of mass system of the two electrons. This 
distribution has been calculated by one of us (11) with the Born approximation, 
taking also into account the screening effect at high energy. If we call the 
angle of divergence of the pair, the value of Q can be obtained from the 
experimental data with the formula: 


(2) Q=M:—P:—=(E, + E,)? — (Pı + Po)? = 
- E E m?ctk? 
~ 2m?c* + m°c4 É a cos © + 2H, E,(1 — cos w) = BE, + E,E,w° 


where mc? = 0.51 MeV. 


(2) R. Jost, J. M. LUTTINGER and M. SLOTNICK: Phys. Rev., 80, 189 (1950). 
(9) G. E. Mopesırr and H. W. Kocu: Private communication. 
(4) A. BORSELLINO: Phys. Rev., 89, 1023 (1953). 


the theoretical distributions. 
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100 


30 50 100 200 300 


Energy distribution of photons. + 
Q <1.5 (with screening); © 1 <Q < 1,5 (without 
screening); @ 1.5<Q <2.5 (with screening); 


<Q2.5 (without screening). 


In Fig. 3 the Q-distribution of the measured pairs is given togheter with 
These were obtained taking into account. the 


hy Mev 


200 


Fig. 4. - Theoretical cross-sections for pair production. a) integral; b,) @=1.25 (without 
screening); b,) Q = 1.25 (with sereening); €) Q = 2 (without screening); 
©) Q = 2 (with screening). 


actual energetic spectrum of the pairs and utilizing both the cross-section 
without screening and that with complete screening (form. 7, 8, 10, of ref. (11)). 


The distributions are 
normalised at the same 
total area. The exper- 
imental hystogram ap- 
pears a little narrower 
than both the theoret- 
ical curves, but the ge- 
neral shape is substan- 
tially in agreement. 
Another type of com- 


parison has been made 


by considering groups 
of pairs characterised by 
the same Q, and by trac- 
ing back for each group 
to the generating pho- 
tonic spectrum. This 
was done by using the 
cross-section correspond- 
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ing to each fixed value of Q. The obtained spectra ought to coincide bet- 
ween them and also with the spectrum obtained by means of the total 
cross-section from all the pairs. In Fig. 4 we give the theoretical cross-sections 
for the values of Q equal to 1.25 and 2; in Fig. 5 the spectra obtained from 
the groups of pairs with a Q value between 1 and 1.5 and with a Q value 
between 1.5 and 2.5. One can say that the spectra, obtained indipendently 
from the different groups coincide between them in the shape as well as in 
the absolute values. This test is satisfactory from the point of view of the 
consistency of the experimental data as well as of their theoretical treatment. 


4. — Distribution of the angle of divergence of the pairs. 


Many authors (%7-215) have already compared the previsions of the Bethe- 
Heitler theory in respect of the angles 9, and 0_ between the electrons and the 
photon or of the angle @ between the two electrons. 

The theoretical cross-section for the distribution of the angles 0, (or 6_, 
as in the Born approximation these angles occur simmetrically) is not com- 
pletely known. The following formula given by BETHE (15) was often used 
for the probability of an angle between 9. and 0, + d9,): 

63 + 0° 
(3) o(0,) 40, = C Pe (in M0 i 


where 9, = me/H,; B and C are factors not explicitly stated. Instead of (3) 
there has been also used the formula derived from it leaving out the loga- 
rithmic term. We remark that these formulae are valid for small angles, so 
that one can approximate sin 0 with 9. Now, while this is always the case 
at high energies, at low energies we can also have large angles and for this 
reason the use of (3) is questionnable; it seems that this has not been suffi- 
ciently noticed by the authors that worked at low energies. 

In the conditions in which the results we are discussing have been obtained, 
the only angle experimentally accessible is the angle w of divergence between 
the two members of the pair. If we call «, the angle kme?/H,E_ and x the 


12) [. Srwons and K. ZUBER: Proc. Roy. Soc., À 159, 383 (1937). 
) A. Z. HRYNKIEWICZ: Acta Phys. Pol., 10, 251 (1950). 

R. R. Roy: Proc. Phys. Soc., A 62, 499 (1949). 

I. B. McDrarmib: Can. Journ. Phys. 31, 337 (1953). 

H. B 


prune: Proc. Cam. Soc., 30, 524 (1934). 
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ratio w/w), we have for the cross-section as a function of x (1): 


x dx 
aay’ 
where a = E,/k and F is a function of x and a as given in ("). 

In Fig. 6 we represent the «-distribution as evaluated from (4), where 
a=0.5 and k=100 MeV which corresponds approximately to the average 
energy of the pairs (in any case the distribution is little sensitive to the 


(4) do = 16a,a(1 — a)F 


Arbitrary units 


ie) 0,5 1 1,5 2 25 3 3,5 4 4,5 5 SS 6 


Fig. 6. — Distribution of the angle of divergence » between the two members of the 
pairs. On the abscissa: ratio m/w), where m) = kme?/(E;E_). a) Distribution calculated 
with the formula x/(1 + x?)?; b) distribution calculated with the formula (4) of the text. 


variations of the total energy if this is sufficiently high). In the same figure 
we represent also the distribution calculated with the formula of Bethe 
æ|[(1+a?)?, which is obtained from (4) where F is regarded as a constant. 
One can see that, as an effect of the variation of F, the maximum of the 
distribution is not at the point 7 = 1/4/3 as we have in the case of the Bethe 
formula, but it is shifted to the right. The experimental points shown on 
the diagram were obtained by selecting the pairs with 0.25 < «a < 0.50; in 
fact only for these pairs the distribution is practically indipendent of a. In 
fact the position of the maximum varies less than 10% between the two above 
mentioned extreme values. The number of pairs here considered has been 
furthermore reduced to 107, because those with a too large error on w were 
disregarded. (See discussion in Appendix). The experimental and the theore- 
tical curves are, as before, normalizsed in the same area. 

We can remark that from this comparison one also finds an excess of pairs 
with small values of x, which can be referred to what has already been said 
for the Q-distribution. The limitation of the statistical data do not permit 
a more detailed analysis. Yet one can say perhaps that the Bethe formula 
gives à too narrow distribution for the small angles. On the other hand for- 
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mula (4) gives a distribution shifted in the opposite direction. It would be 
interesting to make a comparison utilising more precise experimental data 
and more completely elaborated theoretical expressions, i.e. a formula that 
gives the distribution as a function of the only «, eliminating the dependence 
upon the other parameters. 


We wish to express our thanks to Prof. M. ScHONBERG and Prof. G. 
OccHIALINI for helpful discussions and to Prof. M. Cosyns for some sug- 
gestion on the method discussed in the Appendix. 


APPENDIX 


The angles between the two members of a pair have been corrected for 
the effect of the scattering on the tracks. In fact the experimentally measured 
angle derives from the superposition of the real angle @ and the angle due 
to the scattering of the electrons. The effect of the scattering is to modify 
the angle of the tracks, defined as the angle between the cords from the point 
of attack A of the pair, and it is likely to increase it. For small angles the 
probability that the angle is increased is of the order (42 — Q)/4r, where © 
is the solid angle defined by the aperture of the pair. 

If 9, is the mean angle of scattering of the track on the length d, the angle 
between the cord and the initial direction of the track is 0,/./3 ("). 

If we call 1 and 2 the intersections of the tracks, in the case of absence 
of scattering, with a sphere of a center A and radius d, in consequence of the 
scattering each intersection shall be found around the above said points in 
circles of radii respectively 0,/43 and 0,13. In view of the indipendence 
of the two points inside the two circlés, we May assume ru = (0, — Ns); 
where 7? = (6,, + 0,,)/3. On the other hand when @ x, = Nr, WE took as a 
first approximation ea = Ns: 

In any case the correction on the values of Q was very small. 


(17) B. Rossı and K. GREISEN: Rev. Mod. Phys., 13, 262 (1941). 


RIASSUNTO (*) 


Ë stato fatto un confronto fra la teoria della produzione di coppie di elettroni da 
parte di fotoni e i dati sperimentali, per energie fino a 400 MeV. Sono state analizzate 
le seguenti questioni: 1) la ripartizione dell’energia fra i due elettroni nella coppia; 
2) la distribuzione statistica dell’impulso ceduto alla coppia; 3) la distribuzione del- 
Vangolo di divergenza della coppia. 


(*) Traduzione a cura della Redazione. 
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The Decay of the Space Charges 
in Intermittent Discharges in Neon and Argon. 


D. BRINI and P. VERONESI 


Istituto di Fisica dell’ Universita - Bologna 


(ricevuto il 14 Ottobre 1953) 


Summary. — The influence of the gas pressure and of the distance of the 
electrodes on the process of the decay of space charges in intermittent 
discharges in neon and argon is examined. The decay time appears in both 
gases to be of the order 103$, and are, for the distances and pressures 
examined, sensibly proportional to the pressure. The possible causes of 
the high values of such times are discussed, and the various parts of the 
static characteristic of a discharge tube are related to the functioning 
of counters and periodic discharges of the Ist and 2nd kind, which might 
there exist. 


1. — It is known that when a discharge tube containing a gas at low pres- 
sure inserted in a circuit of the type shown in Fig. 1 (*), intermittent discharges 
result. VALLE’s (1) studies have shown that three typically different situations 
can be realized by varying one of the three parameters of the external circuit 
1D Ve OS 


i) discharges in counter region; 
li) periodic discharges called «first kind » discharges; 


ili) periodic discharges called «second kind » discharges. 


Counter region is characterised by discharges which follow each other at 
irregular intervals. Such discharges occur when the gas is irradiated by means 
(*) Resistance R can conveniently be substituted by a diode in saturation, whose 


equivalent resistance can be varied by varying the heater current. 
() G. VALLE: Phys. Zeits., 27, 473 (1926); 48, 473 (1942). 
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of an external ionizing agent (the intensity of cosmic radiation at sea level 
is sufficient). The statistical frequency of these discharges increases with the 
increase of radiation intensity up to à maximum value: beyond this value 
the discharges become periodic, and any further increase in radiation intensity 
no longer produces any effect. With this intensity remaining constant more- 
over, the statistical frequency increases as the e.m.f. E increases. For a 
certain value of Æ, independent of the amount of irradia- 
tion, the discharges become periodic in character, instead 
of statistical, and the corresponding frequency is a cha- 
racteristie of the tube, called « block frequency » (?). Block 
frequency is also the lower limit of the frequencies of the 
periodic discharges of the 1-st kind, when the tube is only 
slightly irradiated. By varying the parameters of the 
external circuit, such frequencies can be continuously in- 


creased, until an unexpected change in the situation does arise, the frequency 
becomes abruptly higher, while the starting potential of each discharge di- 
minishes considerably and transition to 2-nd kind discharges takes place. 
In the field of 2-nd kind discharges also the frequency can be eradually in- 
creased, with a suitable variation in the parameters of the external circuit, 
until transition to continuous glow discharge is arrived at. Very often the 
transition between 1-st and 2-nd kind takes place through an intermediate 
form known as «group discharges ». 

The following interpretation has been suggested by VALLE [l. €. (*)] for 
the above mentioned phenomena: 


1) Transition to 2-nd kind discharges is due to the non complete disap- 
pearance, in the time interval between one and its following discharge, of the 
space charges which deform the field between the electrodes during the dis- 
charges: it follows, therefore, that every discharge starts when the field is 
still deformed and therefore at a potential less than the normal one. 


2) In the case of 1-st kind discharges instead, the time interval between 
two successive discharges is greater: what remains each time of the preceding 
discharge is then only a residue of ionization sufficent for eliminating the 
counter effect. 


3) When the tube functions as a counter, the time interval between two 
successive discharges is still greater, and in fact sufficient to determine the 
complete deionization of the interelectrode space, so that the discharges can 
start only with statistical frequency, dependent upon the tube irradiation. 
Such irradiation, however, must not be so strong as to eliminate, with the 


(2) G. Dascoza and G. VaLLE: Atti Ace. Scienze Torino, 69. 199 (1938-39). 
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ionization produced, the counter effect, as happens in 1-st kind discharges 
in virtue of the internal residual ionization. 

There are, therefore, two time intervals which are of interest to draw a 
schematical picture of these phenomena: 


a) the time t,, necessary for the disappearance of space charges de- 
forming the field, since the discharge has dead down; 


b) the time t,, which is obviously greater and necessary for the complete 
deionization of the gas. 


When the interval between two successive discharges is equal to t,, the 
transition from 1-st kind to 2-nd kind occurs; when it is equal to 7,, the 
transition from counter region to that of the 1-st kind takes place. If the 
duration of individual discharges is negligibly small in respect of these times, 
as is generally the case, 7, and 7, identify themselves with the reciprocal values 
of the corresponding transition frequencies. Therefore 7, is equal to the reci- 
procal of the block frequency, while 7, is equal to the reciprocal of the frequency 
of the end of 1-st kind, i.e. the reciprocal of the highest frequency of the 
1-st kind, as will be clarified later. VALLE has found that 7, is of the order 
of 10-3 s: VALLE and DASCOLA have found that 7, is of the order of 102$. 

Considering the interest presented by a greater precision of these times 
(or of the relative frequencies) in their dependency on the various factors 
capable of influencing the respective values, we have examined these pheno- 
mena by means of techniques available to-day. The present paper gives the 
results of a first series of measurements, regarding in particular the transition 
from the 1-st to 2-nd kind in neon and argon. 


2. — The frequency was measured by a scaling circuit, either, alone, or in 
conjunction with a counting rate meter (Fig. 2): use was also made of the 


m 


oscillographic method, with a gauged time marker (Fig. 3). With the aim, 
moreover, to observe in à single experiment 1-st and 2-nd kind discharges 
and the continuous transition, and to be able to make direct comparisons 
between them, a particular method was used (Fig. 4). This method consists 
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in obtaining frequency variations with a constant R and C, and an e.m.f. E 
increasing at a suitable rate. The required succession of phenomena was thus 
obtained and the periodic repetition of this succession with a rate low enough 
with respect to that of the discharges. The e.m.f. E is obtained from the 
relaxation oscillations of a thyratron: their rising part is practically linear 
and lasts from 1/50-th to 1/20-th of a second. With this method the oscillo- 
eram reveals in fact the time 
behaviour of the potential dif- 
ference at the electrodes of 
the discharge tube. The oscil- 
lations between maximum and 
minimum, corresponding re- 
spectively to the starting and 

: a ait ; Fig. 4. 
extinguishing of the discharges, 2 
follow each other with a fre- 


queney which increases slowly from beginning to end of the oscillogram (Fig. 5). 
But the variable which here in reality is of interest is not indeed the time, 
but rather the e.m.f. E, which is a funetion of it and derived separately 


Eıo,25r 


from the saw-toothed oscillations of the thyratron. A single oscillogram al- 
lows us to read out directly (and at a elance) the starting potential, the mi- 
nimum potential and the frequency of the intermittent discharges, both of 
the 1-st and 2-nd kind according to the values of C and R, as a function 
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and in a wide range of H-values. The method described here, does not how- 
ever permit to include in the same oscillogram the counter region of 
functioning and its transition to 1-st kind discharges, since the corresponding 
frequencies are of the same order of those of thyratron oscilations. 

In our case the repetetive rate of the discharges is gradually varying from 
left to right (# increasing); the measurement of such repetitive rate is how- 


Fig. 6. 


ever possible by using the time markers; this comes out to be particularly 
simple, if one.refers to the stationary points on the locus of the time marks. 
These points, on the other hand, can be easily moved along the same oscillo- 
gram, by modifying 

the fixed frequency of 

runter) the marker (Fig. 6). 
The measurements we- 

re carried out with a 

em tube (Fig. 7) with co- 

Rise axial iron filiform elect- 

rodes, of a diameter 

of 1.5 mm, one of which was mouvable for the purpose of varying the distance 
without opening the tube. The choice of filiform electrodes was due to the 
fact that, with such electrodes, intermittent discharges, in addition to starting 
with comparatively lower potentials, proved much more stable and regular 
than with other forms of electrodes. The results which were port were obtained 
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with the first of the methods quoted: the end frequencies of the 1-st kind 
were determined by varying the C capacity. 

We first of all measured the end frequencies of the 1-st kind using neon 
(99.5% purity) in a vial at pressure varying between ~ 40 torr and 1 torr. 
The e.m.f. was 500 V. The 
reciprocal of the frequency T 
was taken as 7,, while di- 5kx10°s | 
scharge times were neglect- 
ed, because they always 
resulted of the order of some 13 
percent of 7,. By plotting 7, 
versus distance of the elect- 
rodes, for different values 1 
of the gas pressure, we get 


19 roreC+) 


- 45 rorr Ce) 


| 


4torr (A) 


the curves of Fig. 8; the 
same results are also plotted 9 
in Fig. 9, as a T1 versus p A 
diagram. 


2,5torr Co) 


A second series of mea- 
surements was carried out 


in a similar way with ar- 
gon. The e.m.f. used was 
1000 V. We limited oursel- 
ves to exploring low pressu- 
res only, because, for higher 
pressures, the measurements 
did not prove to be easily 
reproducible. The curves in 
Fig. 10 and 11 show the ex- 
perimental results obtained. 

From an observation of 
these diagrams it can be 5 
said that the values of 7, are related to pand d by a formula of the type 


a 
PEA ie 
po fee 


5 10 15 


T1 = Kpd? 


were K seems to be characteristic of the gas, all other conditions of measure- 
ment (tube, electrodes, external circuit) being equal. : 


3. — Let us now consider the static characteristic of a discharge tube 
(Fig. 12). Our problem involves consideration of regions 1, 2 and 3 of the 
plan (i, V) (plan of the stable and transitory phases of discharge): field 2 
contains the static characteristic of Townsend’s phase, while in field 3 appears 
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that of the glow discharge. Experiment (relative in particular to G.M. counters) 
shows that in field 1, there can be no stable discharge conditions if the tube 
is not irradiated: in such a case, therefore, the part corresponding to the static 
characteristic is lacking. With sufficiently high radiation, the discharge current 


30 mm Ce) 
25mm Co) 


10mm (0) 


35) 40 45 p Crorn 


can be stabilized, however, even in field 1, and the static characteristic cor- 
responding there to comes out to be of the form of the dashed curve ABCD. 
In the stage of the curve AB the conditions are the same as in ionization 
chambers: in the stage BC, proportional multiplication appears, and in the 
CD stage Townsend’s phase is found stabilized by irradiation. The stages AB 
and BC of the curve depend on irradiation intensity, because the discharge 
current is not autonomous. The stage CD, however, does not depend on 
the amount of irradiation, provided that this is such as to ensure stabilization 
of the discharge, which in such conditions is autonomous. 

When radiation is very weak in field 1, transitory processes of discharge 
can only develop, pulses statistical in form and frequency: these depend gene- 
rally on the nature and intensity of the radiation which causes them, and on 
the way in which the tube is fed. 

In our experimental conditions at variance with the case of usual G.M. 
counters, the values of the discharge current extend not only to region 1, but 
also to region 3, so that pulses develop each time even as far as glow discharge. 
After a pulse is exhausted, the successive one is started while the potential 
at the tube electrodes is crossing, with a rapidity corresponding to the time 
constant of the supply circuit, the points of the static characteristic. It is 
now obvious that this can happen, as the case may be, either in field 1 at 
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the points of a dynamic charac- 
teristic d,, which are above the 
static characteristic CD, or in field 
2 at the points of a dynamic cha- 
racteristic d,, placed above the 
part DE of the static characte- 
ristic, or finally at the points of 
a dinamye characteristic d,, placed 
above the parte FGH of the static 
characteristic. If radiation is 
weak, the first case corresponds 
to the behaviour of the tube func- 
tioning as a (glow) counter, while, 
whatever may be the radiation 
amount, the 2-nd and 3-rd cases 
correspond respectively to the 
periodical discharges of the 1-st 
and 2-nd kind. 

In glow-counter functioning, 
every pulse starts on one of the 
possible d, characteristics at an 
overtension which is different in 
every single case and statistically 
distributed: stronger irradiation 
has the effect of causing every 
difference between the cases d, 


x4 torr (x) 
2 je ai ae JL 
ARGON | (= 
21 torr (+) 


(e) 


rr(o) 


and d, to disappear, thus eliminating counter effect. The succession of the 
three types of functioning of the intermittent discharges seems therefore, to be 


essentially determin- 
ed, both as regards 


30mm (0) the field of exist- 
ence and the start- 


ing potential of each 


kind, by the behav- 
iour of the static cha- 


nn racteristic of the dis- 


charge tube which is 
20mm (2) being used. The up- 


per prophile of the 


x— 15mm (x) 


10mm (+) 


4 p(torr) - Fig. 11. 
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oscillogram in Fig. 5 and 6 resembles in fact apparently, although only quali- 
tatively, the DEFG development of the static characteristic: but the ab- 
scissa is, in the two cases, very different, and the starting potentials are 


not, for the 1-st kind, those 
of the static characteristic, 
but higher. The appearance 
of supertensions depending 
on the e.m.f. and the fre- 
quency, even in the starting 
of the intermittent dischar- 
ges of the 1-st kind, has al- 
ready been proved by MAuz 
and SEELIGER (*) and interp- 
reted by these authors as 


due to a statistical delay in 
; the discharge. But actually 

this statistical character of 

the supertensions at the 
start of 1-st kind discharges is enterely lacking. For each particular fre- 
quency the starting potential of all partial discharges is in fact in this case 
the same. The fact that such a potential is higher than that on the static 
characteristic is merely due to the regular development of the dinamie pro- 
cess, which repeats itself identically in each pulse. 

Now the natural limits of the types of intermittent discharges under 
consideration appear well distinguished on the static characteristic simply by 
the mere limits of that part of it which refers to Townsend’s stable phase 
(points D and #). If the pulses start to the left of point D counter functioning 
w ll follow, with supertensions and intervals between the pulses statistically 
fluctuating. If the pulses start between D and Æ the normal case of 1-st kind 
discharges takes place: if these start to the right of point Æ the case of 2-nd kind 
discharges arises. Consequently, in order to determine the deionization time 
tT, and the decay time 7, of the space charges, one must respectively take into 
consideration the block frequency, or that of the beginning of the 1-st kind 
(starting point of the periodic pulses near D), and the end frequency of the 
1-st kind (starting point of the periodic pulses near Æ): this was in fact the 
procedure followed. 


Fig. 12. 


4. — Every attempt to interpret the disappearance of the positive space 
charges and to account quantitatively for the experimental results, must neces- 
sarily take into account the different processes which take place during this 


(5) E. Mauz and R. SEELIGER: Phys. Zeits., 26, 47 (1925). 
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phase of the intermittent discharge. These processes can be considered to be: 


1) migration of the charges under field action; 
2) migration of the charges owing to diffusion; 


3) recombination. 


In so far as regard point 1), if regeneration processes are not admitted, 
simple calculations conducted by considering the transit time of the positive 
ions in the average field of the cathodic space, give times smaller by à factor 
10-2 of the experimental ones. Moreover, with this hypothesis, the times 
referred to above would result to be independent from the distance of the 
electrodes. If therefore we want to maintain that the disappearance of the 
space charge in the cathodic space is only due to the field action, we must 
admit that regeneration processes play on important role, by partially com- 
pensating the charges that have disappeared, with a factor y(e*°— 1) always 
less than 1, but quite near to 1 (where J, is the cathodic space). 

Although that can explain the high experimental value of t, it would not, 
however, be able to account for the strong dependence of 7, on d. This strong 
dependence of 7, on d seems to suggest that the processes studied may also 
take place in the entire space between the electrodes. 

As for the other two points, a reasonable evaluation of their quantitative 
contribution seems even less simple. The data reported by RUMPF (*), rela- 
tive to recombination influence, refer to a well-defined particular case which 
cannot be reasonably applied in our processes. 

It seems, therefore, difficult at present to work out a reasonably simple 
theoretical scheme, which could be easily compared with our present exper- 
imental data. 


It is a pleasure to express here our thanks and deepest gratitude to 
Prof. VALLE, Director of the Institute, for his guidance in matters of research, 
as well as for his valuable suggestions and, moreover, for his counsel in advising 
us to offer the results of such studies to the Conference on Ionization at Oxford 
in July of the current year. At the same time we express our gratitude to 
Mr. L. Testont, chief of the technical works of the Institute, for the con- 
struction of the experimental equipment. 


(4) E. Rumer: Ann. d. Phys., 50, 66 (1921). 
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RIASSUNTO 


Viene studiata l’influenza della pressione del gas e della distanza degli elettrodi 
sul processo di decadimento delle cariche spaziali nelle scariche intermittenti in neon 
ein argon. I tempi di decadimento risultano in entrambi i gas dell’ordine di grandezza 
di 10-*s e sono, nell’intervallo esaminato, sensibilmente proporzionali alle pressioni. 
Si discutono le possibili cause dell’elevato valore di tali tempi e si mettono in rela- 
zione le varie parti della caratteristica statica di un tubo da scarica con i regimi di 
contatore e di scariche periodiche di prima e di seconda specie, che vi possono sus- 
sistere. 


IL NUOVO CIMENTO Vor. X, N. 12 1° Dicembre 1953 


On the Scattering Theory of Relativistic Equations. 


E. CORINALDESI (*) 


Division of Physics, National Research Council - Ottawa, Canada 


(ricevuto il 19 Ottobre 1953) 


Summary. — Parzen’s method of evaluating the phase shift at infinite 
energies is extended to the case when the potential has a pole at the 
origin. 

Introduction. 


An interesting investigation of the properties of the phase shift in the 
scattering theory of the Dirac equation was earried out some time ago by 
PARZEN, in the first paper of a series (1) concerned with the use of high energy 
electrons for the study of nuclei. PARZEN showed that, at relativistic energies, 
new features are encountered, quite unlike those of the scattering theory of 
the Schroedinger equation. One important result was that the phase shift 
for the Dirac and for the Klein-Gordon equation, instead of vanishing in the 


limit k— co, tends to the finite value | V(r) dr, if the potential V(r) has no 
0 


pole at 7 = 0, and decreases at large distances more rapidly than the Coulomb 
potential. In the second paper of the series, PARZEN added the extension of 
his formula to the case of potentials which behave like 1/r at infinity, but 
which have still no pole at the origin. 

It is the purpose of this note to consider precisely the opposite case, 
i.e. when V(r) does have a pole at the origin, but decreases at large distances 
more rapidly than 1/r. It will be shown that, on this assumption, the phase 
shift increases like «_,Ink for k — oo, the coefficient &_, being the residue 
of the potential at r= 0. 


(*) National Research Laboratories Postdoctorate Fellow. 
(4) G. Parzen: Phys. Rev., 80, 261, 355 (1950); 81, 808 (1951). 
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For the sake of generality we shall give formulae for the asymptotic value 


of the phase shift of the Dirac and of the Klein-Gordon equation, for a po- 
tential behaving like 


(1) V = + % + Or), 


for small 7, and like 


par uf 
(nr) ee = 
(IT) Per OSs 
for large r. 
1. — Phase shift of the Klein-Gordon equation. 
The equation 
(1.1) G + G 2EV + V?- ar Gr 0 


with H —4/k2-+ 1, is, in special units, the radial K.-G. equation for angular 
momentum /. If V(r) satisfies eq. (I), eq. (1.1) has a solution behaving like 


(1.2) G=r + O(n), 


for small 7, with 


(1.3) DO] de Aiea (OT ce) seed one 


while for large r 


| 
(1.4) G = A,(k) sin (i ne : In 2kr — 9.) 
with 
(1.5) b = EB... 


Eqs. (1), (2), (4) define the phase shift 6,(k) and the amplitude A,(k); re- 
strietions on &_, are easily derived from eq. (1.3) and the requirement that G,/r 
be regular at the origin. 

In order to apply a modified version of Parzen’s method, we write eq. (1.1) 
in the form 
(1.1’) GE G MA + 1) 0) TE 


yr? 
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where 

(1.6) iy re 

and 

(1.7) U(r) =~ +2EV—Vr+ Nee 


Note that U(r) has no singularity at the origin. 
_ We then consider the auxiliary equation 


ne 
(1.8) Goo" + G eset +) GO — 0, 
y? 


| 


which has a solution (2?) behaving like 
(1.9) Go = (1 + Ofr)) 
near the origin, and asymptotically 


(1.10) ao = A(k) sin | kr — 2 = „In 2kr — arg rs + À + 1)|. 


Following PARZEN, we now introduce the function F°, 


et 


(1.11) Re Go! — Go, 


We define also the functions S® and 7\° by the formulae 


1 
Go = 5 (So eile T) à 


(1.12) , 
| F0 = =: CS = LEQ) ; 


\ 


Note that 7 is the complex conjugate of 8S, so that in the following it 
will be sufficient to give relations referring to the latter. Thus, for small 7, 


(1.13) so = ee pi. (1 ae: =) pH +, 
and asymptotically 


(1.14) 8 — — i4f°(k) exp i (i - a — „In 2kr — arg rs +A+ 1))f- 


(2) H. BETHE: Handbuch der Physik, 24/1, p. 289. 
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The functions S and TYP satisfy the equation 


(1.15)  2k8t0 =) Qk? + = | MA + 2S Kl + =| so 


| Be 1) ‚2a. A4 +1)— +1) 


j 
r Y ne 


and that obtained by taking the complex conjugate. 
Similarly, from the solution @, of the original equation (1.1), we define 
the functions S, and 7, by means of the formulae 


L +] 
—kF, = G— = G,, 
1 Y IT 
(1.16) G= 5 (8, + Ti) 
ji, = : N fe 
I een 1) 


Such functions satisfy the equation 


L2EV re) T,, 


ik 
(1.17) Di! = (— 2k? + 2EV. pag, + (ARC à) 


and the complex conjugate. 
Let now 
j S, = exp [io,]8® , 


ak ; 
vie) | T, = exp [—ir,]T, 


so that o, and r, are clearly complex conjugate, and satisfy the equations 


tainty I+1 To 
10, = dk U(r) + ey yo (exp [ UG, | T)| 1) | 
1 To 20, WERL) ASE) 


(1.19) | 


de sw [BY V*) exp [— i(o, + 1) —— + 
7 ı 


r y2 ? 


and the complex conjugate. 
From eqs. (1.4), (1.14), (1.18) it is seen that, for large r, 


(1.20) or) ~ — 6,(k) + (A- 5 + = 2 In 2kr + 


+ arg: rs x? - 2) En 
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is obtained by changing the sign of 


The corresponding expression for 7,(7) 
1.2), (1.13), (1.16) it is clear that o, and 


the last term. Moreover, from eqs. ( 
tT, vanish for r= 0. 

Subtracting the second equation of (1.19) from the first, and integrating 
with respect to r from 0 to co, we obtain the formula 


(1.21) lim (6,(k) + (B-1— a1) In 2k) = (À ED) ee ave LU ler) 


k—> 2 
: Pr 
+tim| | (5 un Ja (ce, — B_1) In B)., 


by letting k tend to infinity. The terms which have dropped out contain 
integrals which are convergent at 0 and at co by virtue of the properties of 
the functions involved, and which vanish for k— oo. 

In order to check the above result, we consider the case of the potential 


(1.22) Vir) = 


i.e. a Coulomb potential which breaks off at a distance 7, assumed to be very 
large. This example can be treated exactly by adapting to the Klein-Gordon 
equation the method given in a well-known paper (?). We get 


(1.23) lim (6,(k) — Z In 2b — Ape n= + arg D(A, +1—iZ) + Z Inn, 


k>o 


where A, = A(@_, = Z), in agreement with eq. (1.21). 

Although it would obviously be nonsensical to attempt to use the Born 
approximation for the calculation of the phases at high energies, as they are 
divergent in the general case which we are considering, it is easy to satisfy 
oneself that such a method would yield a result compatible with eq. (1.21). 
In fact, one would obtain an expression consisting of the first two terms of 
the expansion of (1.21) in powers of x_,. This is also true in the case of the 
Dirae equation. 


(3) W. GorDon: Zeils. f. Phys., 48, 180 (1928). 
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2. — Phase shift of the Dirac equation. 


The system of equations 


I 
| Gr en: EV LEN, 
(2.1) 
11 
| FAN een Neem 


with V satisfying eq. (I), has a solution which behaves like 


(2.2) = OA 0m) F,=r*1(e + O(r)), 


(2.3) A+1=V (+1) — 2, 
and 
(2.4) le 

%-ı 


The constant À occurring in this section is not identical with that of the 
K.-G. case defined in eq. (1.3). 
If V also satisfies eq. (II), the asymptotic behaviour of @, and F, will be 


En ls E 
| G,~ A,(k)V E +1 sin (ir — > Den = In 2kr — A) : 


(2.5) 
= l E | 
| Rn — Aılk)VB— 1 cos (er 5 ah 5:0) 


Proceeding as for the K.-G. case, we introduce the functions 8, and T, 
defined by 


VE+1 
Gi == er 
(2.6) = 
NE I 
er 
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and obeying the equations 


PR: : (LS ig lee ee 
| = m ars +( re, 
2. 
Sei h a ts pie O~ 
Peak RP) ES | V8 
k r k 


Use will also be made of two auxiliary functions S\? and T°, which are 
solutions of the equations obtained from (2.1) on replacement of V by œ:/r. 
Such functions will behave according to (2.2) near the origin, and asympto- 
tically 

— l 
GO ~ A(k) VE + 1sin (ar — 5 on 
23) == Ix E 
| FO ~ — AM(k) VE — 1 cos (er — 5 Hr In 2kr — so) D 
with (4) 
(2.9) exp [— 2369 |= 
dan —(L+ 1k TA +2 + ÜBlk)e-,) 
— Bas + (À + 1)k P(A + 2— iElk)e_,) 


exp [— mi(A—1+1)]- 


The functions S® and 7, expressed in terms of G® and Fi” by relations 
identical in form with (2.6), obey the equations obtained from (2.7) by substi- 
tuting «,/r for V. The remaining part of the procedure, which consists in 
defining o, and 7, as in eq. (1.18), and in integrating with respect to.r from 0 
to co the equations for these new functions, runs precisely as in section 1, 
and it seems superfluous to reproduce it here. The final result is 


(2.10) lim (814) + (B-1— %1) In 24) = (A— N= + arg P(A + 1—iaa)+ 


be Sin; He )+ im | | (r— =) dr + (a, — 64) In R}- 
Za L. 1 


3. — Discussion of results. 
We can distinguish several special cases: 

1) &_,—=ß-ı. According to eqs. (I) and (II), this is not necessarily the 
case of a purely 1/r potential. At high energy 6,(k) tends to a finite value. 
The integrals in (1.21) and (2.10), which depend essentially on the potential 
at finite distances, vanish identically if V= a,_/r everywhere. 


2) a,=—0, B #0. For this case, already considered by PARZEN, the 
(4) N. F. Morr and H. S. W. Masspy: The Theory of Atomic Collisions (Oxford. 
1949), p. 79. 


1680 E. CORINALDESI 


phase shift diverges like —- ß_,Ink for k > co. This logarithmic divergence 
is of a rather spurious character, and depends on the special form adopted 
in eqs. (1.4) and (2.5) for the asymptotic behaviour of @,. It arises from the 
fact that V has not a polar singularity at the origin to match the 1/r decrease 
at infinity, so that, at very high energies, the phase may differ considerably 
from the finite Coulomb-like phase in case 1). 


3) a, #0, B_, —0. This case was not considered by PARZEN, probably 
on account of its irrelevance to the problem of scattering of electrons by nuclei, 
the latter being usually represented as continuous distributions of charge pro- 
ducing potentials which have no singularity at the origin. It might be of 
interest, however, in problems of atomic physics, such as the scattering of 
high energy charged particles by hydrogen atoms in the one-body approxi- 
mation, with (5) 

“ 1 1 

V(r) =— & E + N exp [— 2r/a| - 

nn 
The phase shift ö,(k) diverges like «_,In%k for k — co, so that the slope of a 
graph of 6,(k) versus In k would immediately yield the physical quantity «_,. 
It is not perhaps amiss to consider briefly an important feature of the 
scattering by potentiäls which have a pole at the origin. For those which 
do not (x, = 0), all the phases are equal at high energies, as lim Ô,(k) is in- 
dependent of /. Therefore, the scattered amplitude for k — oo is propor- 


tional to the quantity > (21 + 1)P,(cos 0), which, though divergent, is summ- 
1=0 

able for 6 40, the sum being zero. This implies the existence of an extre- 

mely narrow maximum at 6 = 0 for the scattered intensity. On the other 

hand, when &_, #0, the above property does not hold, and the distribution 

may be expected to flatten out. A more careful study of this point is in 

progress. 


My thanks are due to Dr. G. HERZBERG for enabling me to spend the 
summer months at the National Research Council of Canada, and to Dr. Ta- 
You Wu and my colleagues for their interest in my work. 


(5) Reference (1), p. 184. 


RIASSUNTO (*) 


Si estende il metodo di Parzen per il calcolo dello spostamento di fase alle energie 
infinite al caso in cui il potenziale ha un polo all’origine. 


(*) Traduzione a cura della Redazione. 
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Remarks on the Validity of the Tamm-Dancofi Method (*). 


G. MorpurGo and B. F. TOUSCHEK 


Istituto di Fisica dell’ Universita - Roma 
Istituto Nazionale di Fisica Nucleare - Sezione di Roma 


(ricevuto il 19 Ottobre 1953) 


Summary. — An example of a field which has already been investigated 
by WENTZEL and BLATT is discussed in detail. The exact solution of 
the problem corresponding to the seattering of a meson by a nucleon 
is compared with that obtained by the first approximation of the Tamm- 
Dancoff method. It is found that this method leads to results, which, 
though better than those obtained in perturbation theory are essentially 
qualitative. The disagreement is particularly strong, if the coupling 
constant is nearly sufficient to give rise to bound states. 


Introduction. 


The Tamm-Dancoff method (!) has been recently applied to various pro- 
blems of meson theory; DANCOFF, LEVY (?) and KLEIN (?) have treated the 
deuteron problem, CHEW (‘), DYSON (5) and FUBINI (°) the scattering of mesons 
by nucleons and CINI (?) has given a covariant formulation of the method, 
which admits the treatment of the renormalisation question. The validity 
of the method, however, has not yet been investigated; in the case of the 


(*) Part of this work was communicated at the Cagliari meeting of the Italian 
Physical Society, September 23-27, 1953. 

(1) I. Tamm: Journ. Phys. USSR, 9, 449 (1945); and 8. M. DANCOFF: Phys. Rev., 
78, 382 (1950). 

(2) M. Levy: Phys. Rev., 88, 72, 725 (1952). 


(3) A. Krein: Phys. Rev., 90, 1101 (1953). 

(4) G. F. CHew: Phys. Rev., 89, 591 (1953). 

6) F: Dyson et al.: Proceedings of the Rochester Conference, December 18-20, (1952). 
(6) F. Fuginr: Nuovo Cimento, 10, 564 (1953). 

(7) M. Crnt: Nuovo Cimento, 10, 526, 614 (1953). 
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scattering of mesons the results of the first order approximation have merely 
been compared with the experimental data. 

It therefore appears useful to study in more detail what degree of accuracy 
can be expected from this method. For this purpose we choose an example 
which on the one hand is sufficiently simple to allow an accurate solution and 
which on the other hand should be sufficiently complex to give significant 
results. We shall restrict ourselves to a comparison between the accurate 
solutions of the example, which describes the scattering of a «meson» by a 
«nucleon », and the first approximation of the Tamm-Dancoff method. 

The example which we chose is one which has been used by WENTZEL (8) 
in various publications and by BLATT (°) in an investigation of the validity 
of the Heitler method. It is a field which can be described in terms of the: 
Hamiltonian 


RE SIEH 


yes 


(a*(a) + (Vep)? + wp?) de + 5 9 R Inte joe) a | 


where 


We observe that this Hamiltonian describes the interaction of a «meson- 
field » with a fixed «nucleon-source ». It gives rise to the virtual creation and 
annihilation of pairs of mesons and contains first order matrix elements for the 
scattering of a meson. This last fact simplifies greatly the first order cal- 
culation with the T.D. method. The limitation to the first order approxi- 
mation is in keeping with the treatments hitherto given to the meson scattering 
problem, in which all those states are neglected the introduction of which 
would lead to coupled integral equations. 

Although BLATT has already obtained an expression for the cross-section 
of meson scattering, it appears to us that some qualifications have to be 
made to this derivation, regarding the possibility of bound states as well as. 
the interpretation of the results of his canonical transformation. We shall 
therefore give a derivation of the cross-section in the first four paragraphs. 
In particular we shall treat the diagonalisation of the Hamiltonian in the first, 
and the question of renormalisation in the second. In the third we shall 
discuss the physical significance of the transformed states and variables and 
in the fourth we shall give the actual derivation of the cross-section. The 
5th paragraph contains the solution of the T.D. equation and a comparison 
of the approximate results with the accurate ones. 

(8) G. WENTzEL: Helv Phys. Acta, 25, 569 (1952) and preceeding works. 

(?) J. M. BLarr: Phys. Rev., 72, 466 aye 
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1. — Diagonalisation of the Hamiltonian. 
Preliminary Considerations. - We denote by 


Pr, 1m Ir, 1,m 9 


the operators associated with the oscillator of the component of the meson 
field with energy 


on = VIE +, 


angular moment / and magnetic quantum-number m. Since the interaction 
in the Hamiltonian (1) only acts on the s-component of the meson we put 
Proo = Pry Too = % and therefore 


il ; sin kr 1 | sin kr 
(1) pr = —— |a(2) —— dx = ——— | pie) dé, 
oer ty re, Pa A HAT 
where R is the radius of the sphere of normalisation. It thus follows from (1) 
that the part H, of the Hamiltonian corresponding to the s-component has 
the form: 


1 1 
(2) H,= 5 > (pi: Oni) 59 > Are VrVr! » 
% CRT 


a 


were we have put 


For the following it is convenient to introduce the usual operators of creation 
and annihilation, for which one has, owing to the reality of the q, and px 


il Fr à Ox, 
zn, ara) (VA — (ar — &) - 
20% 


The transformation 
(3) Ck = > Exe Qxr y 
= 
in which &,. is real and satisfies the equation 


(4) er Oy )Erx = 9x > VE ur 
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(where W, is the eigenvalue of the equation) diagonalises (1) on account of 
the orthogonality relation 


> Eure nrs = Ors . 
PR 


For the possibility of the expansion (3) we have also to require that the &x- 
form a complete set: 


> Err& sk = Or; = 
k 


The Solution of (4). — It is easily seen that equation (4) admits a set of 
solutions of the form 


F VgV 1 
2 fun = den + 9P a 798,’ 
where 
> V4 


To obtain (5) one has to put W; = u?+ k? and to assume that k is not a 
solution of 


(7) AMIS 0 


The symbol P in (5) and (6) denotes the principal value, defined by the following 


f-plane 
path of integration in the complex t-plane 1 eS 
| 2 es 
# 


Other solutions of (4) are of course obtained by à different treatment of the sin- 
gularity. BLATT has chosen a solution which represents outgoing waves only 


k 
and which is defined by a path of integration of the form un 


His solution is of course completely equivalent to ours but leads to a complex 
£x, which we have excluded by our choice (3) of the diagonalising transfor- 
mation. 


If k = k* is a solution of (7) it is easy to see that 


Sane == lim Ex 
k—>k* 
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and it therefore follows that the solutions (5) completely exhaust the con- 
tinuous spectrum. 

There are, however, other solutions (4) corresponding to bound states of 
the system. Putting k = ik, in equation (4) it is easy to see that if 


(8) L+9 game 


has roots |k|< y (and no roots with |k|> u, in which case the system, due 
to the appearance of run away solutions, probably cannot represent a phy- 
sical field) there exist solutions of the form 


Vp! 
(9) Erin, — EL RZ Cx, 5) 
where Ce, is a constant of normalisation. The Fourier transform of wn, va- 
nishes exponentially for r—co. The solutions (9) can therefore be interpreted 
as belonging to bound states of energy nW, (Wy = APR WT, 25-2) OE 
the meson field. 


Orthogonality and Completeness. — The proof of orthogonality can be ob- 
tained immediately by inserting from equation (5) and (8) for the discrete 
states. 

The demonstration of completeness is more complex and we have only 
succeeded in showing the validity of the completeness relation for the case 


(10) 198% |< i, 
for every real k. The demonstration will be given in appendix I. 


The Diagonal form of the Hamiltonian is obtained by putting 


(11) = 


1 .1/ Wi 
Pr oi) P, = il: (er Ay) - 


The transformed Hamiltonian then has the form 


(12) Hi= % 3 Wlaat + ai) = >dafa.W. + % > Wr. 
k k k 


It is useful to note the expressions for the relation between the a, a; and 
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the «,,«,;. These are obtained by using (3) and (11): 


ay = > IN (23 KA = ag Ae) 
= 


(13) , \ 
a > Error A® & + pale 12) 
= 
where 
1 (1/W; Ve = WE eee 
Krk — TER | kik — __ kk’ — = pe À 
iP JR! 5 | read hare and A Ale W. 


The inverse formulae can be easily obtained by using the completeness relations 
for the &,,.. 


2. — Mass Renormalisation. 


The energy of the lowest state recovered from (12) is in general different 
from the mass of the nucleon. If one puts (in the continuous spectrum) 


= Ex 
Wem oy + 20 à 


© 
REES 
Ir 
sel 


the energy of this lowest state becomes 


and the last term gives a contribution even for R— co. If therefore one 
adds to H* the constant 


the energy of the lowest state is just equal to M. The same M’ could have 
been added at the start to the Hamiltonian H. It therefore follows that the 
prescription for renormalisation in our case is just the substitution of H’ by 
H'= H'+ M'(g. We can immediately see that in first order in g this pre- 
scription gives the right result. 

The independence of our field from the nuclear variables implies that in 
our case the renormalisation can be effected simply by the addition of a 
constant to the Hamiltonian. 
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We therefore have to determine «,(g). If 
Il rs 
sinken 
m 
is a stationary free meson wave in the sphere of normalisation and 

Ehe x 

— sin (k,r -+ 0,) 

= 


is the corresponding asymptotical behaviour of the solution in the presence 
of an interaction, one has 


Ön 
ek 2 — |) 
(in lea) Se R 
and therefore 
Ex kn i kn Ön 
(WW. SE > 
PO a well Wi, ae Wr, R 


Using the expressions of the following section 4 for the scattering phases one 
obtains 


(14) GX VE 


If the normalisation procedure has to be carried out by means of a perturb- 
ation treatment, one has to expand (14) in powers of g. The first term el? 
of this expansion is 
2 
a — 9% 
TN. 


and this term is exactly the term needed for the normalization of the first 
approximation in the T.D. method. We shall come back to this point in 


section 5. 


3. — Interpretation of the Theory. 
In this section we want to discuss in some detail the relation between the 
variables of the original Hamiltonian (1), which have à simple physical signi- 


ficance and the variables of the transformed Hamiltonian which have to be 


109 — Il Nuovo Cimento. 
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interpreted in terms of the physical variables. This necessity of interpretation 
occurs whenever a canonical transformation is carried out and it therefore 
appears to us that the following treatment, though cumbersome, may be of 
a methodological interest. To give a greater generality to the argument, we 
shall ignore certain properties of the Hamiltonian, which could be discovered 
by an application of perturbation theory — as e.g. the fact that the matrix- 
elements leading to the production of three mesons vanish on the energy shell 
(at least in first and second order). We shall denote by 7(z,, %i,5--+) (ax 4x Y = 4) 
the eigenstates of the Hamiltonian H, and by X(N,, N;,,...) the stationary 
states of the system (a'a,X=N,.X). We put y(0, 0,...) =v) and X(0, 0,...) = Xo. 

Passing to the case of the scattering of a meson of energy Z,, we shall 
have to consider. a combination of states X, which represent asymptotically 
in-going and outgoing mesonwaves of energy H,. If E,<3u (it has to be 
remembered that in our example mesons could only be created in pairs) only 
the state satisfying a, a,X— X will have this property. This state we shall 
denote by X(1,). Expanding X(1,) in terms of the complete system of the states 
of the Hamiltonian H, one finds 


(15) X (1p) = À UK| ka) xAn,) + I VO) + …, 


ky keykokg 


where U'®(k|k,k,k,) is Symmetrical in k,k,k,. No components with an even 
number of mesons need be considered. The 1-particle component of X(1,) 
is given by U(k|k,), or its Fourier transform 


sin kyr 
(16) PO Ir) æ % —— 


ky 


U%(k| ky), 


which may be interpreted as the wave function of the meson. Since £,< 3u 
only U® can be different from zero asymptotically. It is also clear that U™ 
must show a singularity for k +k, and it is just the behaviour in the vici- 
nity of this singularity that will lead to the expressions for the scattering 
cross-section. 


4. — Derivation of the Scattering Cross-Section. 


Determination of U®. — For U® one may write 


(16) UO (k | ky) = (yo ea Xo) = (X e421 Yo + D UP(O| kk jai xf x0) 
k'k" 


where. U® is defined by the expansion similar to (15) of the ground state of 
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the system : 


(17) X= yo + D UP O|WE\yAede) +. 
k'k" 


U®—symmetrical in k', k'— is therefore the two particle component of the 
groundstate. It should also be noted that the only terms of the expansion (17) 
which contribute to U® in (16) are those actually written down. Using equa- 
tion (13) one finds the following integral equation relating U and ue 


(18) Uu(k 


ky) = Nr k a 3 ConU?(0| ky k’), 
m 


where 
Nik = Ed Mk ; Crk = Er »Ah* : 


An integral equation for U® can be found in the following manner: X, can 
be defined by 
eX 0% 


Multiplying on the left by Xe de, one obtains immediately by using equation (13) 


(19) Ci 2 > ex UP (0 | k,k') — 0. 
er 


Equations (18) and (19) represent a system of integral equations from which 
in principle it should be possible to recover U® and U®. We have not found 
a general solution of this system of integral equations. For the determination 
of the eross-section, however, it is sufficient to study the behaviour of U® 
for k—k,, and this can easily be determined from the integral equations. 
It is perhaps worth noting, that, as our example shows, the complete know- 
ledge of the Heisenberg operators of the system, does not immediately lead. 
to the quantities necessary for the interpretation of the development of the 
system. The problem of interpretation may in general be much more com- 
plicated, leading, as in our case, to a system of integral equations. 


The Singularity of U®. - We now proceed to discuss the singularity of U® 
by means of equations (18) and (19). First it is clear that in the following 
calculations we need not distinguish between W, and w, in the expres- 
sions (13). Taking the sum and difference of the equations (18) and (19) 
one can then easily obtain an integral equation for U™, namely: 


‘Wi, U(k | ky) 
(20) UD(k | A ~ a Ene > RES ! - > Or Erne Fur! + 
rk’ 


x” 
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We try to find a solution of this equation of the form 


f (ky | k) 
k? — ki 


(21) UM Kk] ky) = du, + P 


in which f(k|k,) is a regular function. Remembering that Ex is of the form 


k,|k 
Fink = On,k In piel? ’ 


with g(k,|k) defined by equation (5), one finds that 
(22) fik|k,) = g(k,|k) + X(k, |k), 


where X(k I is a regular function of k,, k and X(k|k) — 0. If now one 
defines 7(k) (|7|?= 1) from the asymptotical behaviour of the wavefunction 
Wo(k|r) by means of 


P(k |r) > exp [— ikr]— o(k) exp [ikr] asymptotically, 
one has 
2iRg(k|\ kb) 
2k + iRg(k|k) ' 


dr 0) — 


For the phase used in section 2 one obtains: 


Rg(k|k) 


exp [2404] = mo(k) ; te 6, = = 


In the paper quoted above, BLATT has determined the eross-seetion simply 
by considering the singularity of the transformation matrix &,x. Equation (22) 
is essentially Blatt’s result. It has here been derived in accordance with the 
more general viewpoint adopted in this section and after the introduction of 
the two conditions 1) E<3u and 2) the coupling constant is chosen such that 
the system has no bound states. The result is probably more general than 
these restrictions. 


5. — The First T. D. Approximation; Comparison of Results. 


In the first T.D. approximation we only consider the one-meson component 
of the generic state vector. Calling w, the amplitude of the 1-meson state 
we obtain the following integral equation 


Ve 


(23) : AO Does cS 
Wy K 


a tn 
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where a term gu, > (v?/4w,) does not appear because it .is exactly cancelled 
LA 


by our prescription for mass-renormalisation. If one would want to continue 
_ this method one would have to introduce suecessively the amplitudes u(k), 
u(k,, ky, ky)... corresponding to the amplitude of 1-, 3, -,... meson-states re- 
spectively. This would lead to coupled integral equations. 

The singular part of the solutions of (23) is formally identical with that 
of the solutions (5) of the equation (4) for the transformation function &r,. 
The expression S, however has a different meaning, namely 


F ir: 
23 SED) — P en u 
(23) i 2 20 (Wx — Wx’) 


For the transition amplitudes one can therefore write: 


a x __ —2iRgo 
1 N Ok — gSo)— iRgr?’ 
cm HR 


1— ne) = oe Wee i = ee 
| le ~ Bk — ge) — Re 

where the indices (e) and (T.D.) respectively indicate the exact and the T.D. 
solution. It is immediately seen that the dissipative (HEITLER) term is the 
same in both expressions, the difference rests with the reactive term. For 
completeness we add the first Born approximation, which is 


; Riera’ 
(24') 1 hee ae ae 
To discuss these expressions it is convenient to consider a particular form 
of v;. 
We take 
ark 1 
(25) VE T—— ; 

vV2nR k? + a? 
which corresponds to a function U(r) = (a?[4x)(e-“ fr). 

This choice is only made for reasons of simplicity. Since our example 
can be considered as representing the scattering of mesons by a nucleon the 
recoil of which is neglected it is reasonable to choose a = M, which is in 
fact equivalent to a cut-off of frequencies — M. In this way one obtains 
for 8; 


3 k— M! 
= Jus (ke + MR: 


(26) se 
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In the nonrelativistic limit k< u < M one has 


we, 


the last line serving to define 6. In the T.D. case one finds 


Mt kak’ 1 
STD) — Ko 4 a ne et a a 
27 ) (k° + M?) 2œoy(œox + Ox’) 

0 
Neglecting terms of the order k?/u? and w?/M? one obtains +f/2 for the 
integral and therefore 
(27) STD) > — 2 
There is however an inequality which / must satisfy in order that our deri- 
vations should be valid. This inequality is connected with the fact that at 
the moment we have only demonstrated the completeness when (10) is satis- 
fied; in terms of ß and with our choice of v,, (10) can be rewritten 


(10') —1<f<l. 


It is however possible that the above restrictions (10’) are too severe and we 
have reasons to believe that the upper limit could be replaced by at least +8. 

’ The lower limit in equation (10') corresponds to the fact that for 6<—1 
bound states with energy nv/p2— k? (n = 1, 2,...) are possible; in fact with 
our choice of v, equation (8) can be rewritten 


IQ | ko | DS 
(28) (3 +1) ZB, 


which if 6<—1 has one (and only one) solution for |k,|. 
Using the expressions (26) and (27) and inserting into equations (23) 
and (24) we now have in the limit k —0 


— 2igk 
1—n® — 3 Er 
7e Anl — B)— igh’ 
— 2igk 
HDI TE ee 
"0 An(1 — B/2) — igk ? 
— igk 
1! = re = tas eos , 
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and therefore in the same limit 


og TD) (al = pb)? o'B) es à 
con A Bj2)? ’ ar (1 + B) 


In the case of positive B (repulsive interaction) the ratio o® /o® compatible with 
condition (10’) is at most 4—the Born-approximation giving an overestimate 
for the cross-section owing to the neglect of the reactive effeet—the corres- 
ponding ratio for the T.D. method is 16/9 and this is due to the fact that the 
reactive term of the T.D. approximation is only 1/2 of the correct reactive term. 

The situation is even worse in the case of an attractive interaction. In 
the neighbourhood of the region in which bound states are possible the ratio 
oP) /6® may be made arbitrarily small. 

We therefore conclude that the results given by the first T.D. approxim- 
ation, although better than the ones of perturbation theory have a qualitative 
character; the results are slightly better in the case of repulsion than in the 
case of attraction. 


APPENDIX 


Substituting from equation (5) into the completeness relation one has 


Urs i oy al À VOTE 
= (6), fo —- a ——: 
2 End SE Bet ee s? À, Ian s? A, 72 — §? 
Dew el Vs Deel 
23 2 2 j g° 2 2 > 2. 2 Rye 
-pP—t A, rs? + 92? A 


For the completeness it is herefore necessary that 


3 v 1 1 i x 
(a) DI um = const =—1. 


r 


The value of the constant must be — 1 because of 1/4, >1 for roo. For 
equation (a) one may write 


1 1 On 


IS à Er A? Sen 


(b) 20 


Expanding 1/4, in a power series with respect to gS, (Supposing that 1gS+1< 1 
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in the interval in question) the proof of (b) is reduced to showing that | 


2 
n v 
(n+1)> 8S} ae, = tt, 


t 
The latter follows from the general relation 


1 1 L 
me Die (a, — 2)... (&n-ı — 


En) — En) (Bm) — a) ... (© — u) 


7 


where P denotes the permutation of the n objects x, ... x, among themselves. 


RIASSUNTO 


Il modello di campo gia considerato da WENTZEL e BLATT, & discusso confrontan- 
done la soluzione esatta con la soluzione ottenuta col metodo di Tamm e Dancoff in 
prima approssimazione nel caso dello scattering mesone-nucleone. Risulta che il metodo 
di Tamm e Dancoff conduce sempre a risultati che, sebbene siano migliori di quelli 
della teoria delle perturbazioni, hanno carattere qualitativo; in particolare, se la costante 
di accoppiamento della teoria è tale da essere in vicinanza di uno stato legato i risultati 
del metodo di Tamm-Dancoff sono particolarmente cattivi. 
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Summary. — A method for the calculation of the adiabatic potential 
between two fixed nucleons is proposed giving a sequence of approxi- 
mations in which at each stage only a limited number of mesons and 
pairs is taken into account. By means of a variational principle the 
adiabatic ponent is approximated with a sequence of fractions of poly- 
nomials in g?, whose coefficients contain only the potentials calculated 
by Livy and KLEIN with perturbation methods. The convergence of 
the sequence does not depend on the magnitude of the coupling constant 
and, when the limit exists, it is an analytic continuation of the series in 
powers of g? outside its radius of convergence. 


1. — Introduction. 


Recent attempts have been made to deduce a potential between two 
nucleons with the use of the Tamm-Dancoff method (1*). However, as is 
well known, the kernel of the integral equation for the two nucleon amplitude 
is obtained by means of an iteration procedure. KLEIN has shown that the 
lowest order potentials obtained with this method give results in strong dis- 
agreement with experiment, and furthermore that the series in powers of g? 
has a very bad convergence. Both these circumstances suggest that a better 
solution of the problem might be obtained by using a method whose con- 
vergence does not depend on the magnitude of the coupling constant. 


(1) M. Lévy: Phys. Rev., 88, 72, 725 (1952). 
(2) A. KuEin: Phys. Rev., 90, Ga (1953); 91, 740 (1953). 
(3) M. RupERMANN: Phys. Rev., 90, 183 (1953). 
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It may be worthwhile to notice that the Lévy-Klein approximation differs 
essentially from a method in which the number of mesons and pairs is limited, 
without making any assumption on the relative magnitude of the different 
amplitudes (%:). A method of this kind has been recently used for the pion 
nucleon scattering problem (47), taking into account only amplitudes with 
no more than two mesons: the fact that the results show a substantial im- 
provement in comparison with the Born approximation suggests that the 
amplitudes with two or more mesons cannot be simply eliminated with an 
iteration procedure. 

In the present paper we propose a method yielding a sequence of approx- 
imations to the nuclear potential in which in each stage only a limited number 
of mesons and pairs is taken into account. The convergence of this sequence 
of potentials does not depend on the magnitude of the coupling constant, at 
least as long as one considers adiabatic potentials. The method can be 
extended in a natural way to the evaluation of the reaction operator Æ (5) 
for scattering problems. 


2. — A compact expression for the power series expansion of the adiabatic nuclear 
potential. 


In a previous paper (?) the bound state integral equation for two nucleons 
interacting through meson field has been deduced in a covariant way by means 
of a tedious and lengthy iteration procedure. In this section a new formulation 
will be given by means of which a compact form of the same integral equation 
can be obtained, and the expression for the adiabatic potential between two 
nucleons readily deduced with the introduction of suitable approximations. 

Instead of starting from eq. (1.1) for the interaction representation state 
vector Y(t) we start from 


+ © 


(1) P(t) = P(— 00) + P(+oo) — i farine ~—t,) dt, , 


(2) 


(is) F. J. Dyson: Phys. Rev., 91, 1543 (1953). 

) G. Cuew: Phys. Rev., 89, 591 (1953). 

) S. Fugini: Nuovo Cimento, 10, 564, 851 (1953). 

(°) F. J. Dyson et al.: Proc. of the Rochester Conference, December 1952. 

) R. Karpzus, M. Kıverson and P. C. MARTIN: Phys. Rev., 90, 1072 (1953); 

K. A. BRUCKNER, M. GELL-MAnN and M. GOLDBERGER: Phys. Rev., 90, 476 (1953). 
(8) J. SCHWINGER: Phys. Rev., 74, 1439 (1948). 
(?) M. Cini: Nuovo Oimento, 10, 526 (1953) denoted in the following as I. All 

symbols not explicitly defined will have the same meaning as in I. 


re 
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where 


— 
bo) — 


fester 


pate 


SE 


In the present problem of a bound stationary state of two nucleons both 
the contributions from Y(+oo) and Y(— co) vanish (7°). 

Let us introduce a projection operator n with eigenvalue zero for the 
states with only two nucleons and 1 for all other states with any number of 
mesons and pairs. Then 


(3) A(t) =—1 fo n)H(t)P(t)e(t — t,) dé, . 


Obviously @(t) is the component of Y(t) describing two nucleons only. 
To derive an expansion in powers of g? of the non adiabatic kernel for the 
two body problem one iterates eq. (2) an odd number of times, and intro- 
duces the value of Y(t) thus obtained into eq. (3): this gives a chain of factors 


(1 — n)HnHn ... nH , 


in which all the unwanted kernels containing intermediate states with two 
nucleons only are automatically ruled out. 

Two further assumptions are necessary to derive the adiabatic potentials 
of Lévy-Klein. Let us define V(t) by means of 


(4) Y(t) = V(t) D(t) . 


(1%) Both the contributions from W(— co) and W(+ oo) oscillate, at t = — co and 
t= + co respectively, for the bound state and a suitable limiting procedure (see M. 
GELL-MANN and M. L. GOLDBERGER: Phys. Rev., 91, 398 (1953)) makes these terms 
vanish. Eq. (1) is exactly equal to eq. (1.1) for the bound state because in the former 
the integration on time gives terms of the form P/(W— Hp — Ep) (P means principal 
value), while in the latter one gets terms 1/(W— Ep — Ep: + in). Since W < Ep + Ep: 
there is no difference between the two. We shall see later the reason of our choice. 
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Inserting this expression in eq. (2) and integrating by parts one gets 


+ © 
(5) V(t)D(t) = 1 =o nH (t')V (t')e(t na) P(t) + 
Mt ! i (/4 It 7 " d i 1/4 MW 
+i any fae MOINE (t! = [art fom nH(t") V(t") m ® 


The first assumption referred to above consists in neglecting the time deri- 
vative of ®(t) in eq. (5) (1). Thus the equation for V(t) is 


— 


(6) VQ | nH(t)V(t)e(t—t)dt,, 


— © 


u | 


=p) V = 17: 


We introduce the representatives of states with two free nucleons, denoting 


them shortly as <. Then eq. (3), with the help of eq. (6) becomes 
(8) son) [au IKB) 4) dt. 


This is the compact form of the two body non adiabatic integral equation. 
Let us consider 


(9) W= | H(t)V(t) dt. 


This operator connects states with equal energy. 
In configuration space its matrix elements, when the nucleon recoils are 
neglected, are of the form 


(10) Wieser) == | r\H()V(G)|r du = 6(r—r')w(r) 


('') The time variation of @(t) is much slower than that of any component of Y(t). 
This corresponds in the usual time independent theory to neglect (W— E, — Ep’) in 
comparison with all other energy denominators. 
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w(r) is the required Lévy-Klein adiabatic potential when V(t) is obtained 
by the iteration of eq. (6). The reason of our choice for eq. (1) is now ap- 
parent: we need an expression for the adiabatic potential W in a form expli- 
citly hermitian. 

The potential (9) contains also all sorts of divergent contributions cor- 
responding to radiative processes. Following Lévy and Klein we will dis- 
regard for the time being these contributions: nevertheless we notice that 
the covariant form of all the operators introduced is appropriate for an in- 
vestigation of the divergences. 


3. — The variational principle. 


We will produce now a variational principle for W, showing that this quan- 
tity is stationary for small arbitrary variations of V(t). The procedure follows 
closely the formalism of LIPPMANN and SCHWINGER (©) for the K reaction 
operator. 

The analogous of eq. (1.42) of reference 2) is 


+ © + 
(11) W'— [mono + H(t)V(t)] dt — | V*(t)H(t) V(t) dt + 
TEE CS 
4 ja | dt’ V*(t)H(tmH(t}V(t'}e(t-—t'). 


It is easy to verify that in eq. (11) W' is stationary and equal to W for 
V(t) satisfying eq. (6). 
We shall introduce the following operator as «trial function », 


(12) V{t) = 1 + V ()4 + V(t) A +... + Valt)An ; 
where 
ae + © 
Bert (= 4)" | dm [ ae, nH (tin --- nH (t,je(t— t eth ts elle sla) 


The A; are time independent operators which will be determined in order 
to make W’ stationary. The first term in expansion (12) is kept equal to 1 
in order to satisfy eq. (7). The introduction of the «trial functions» in 


(2) B. A. Lıppmann and J. SCHWINGER: Phys. Rev., 79, 469 (1950); see also 
M. Cinı and L. A. RADICATI: Nuovo Cimento, 7, 905 (1950). 
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eq. (11) gives for W’ the following quadratic form: 


a) Wi Wet Sf AP Wa + Weeds ZA Wer Wunsch, 
1=1 


= 1 
where 
ap +a 
(15) W,= (ii) | dé, .. [andre … NH (t;)E(ts — tn) ... Eli — ti). 


The variation of eq. (13) yields the following equations 


n 


(16) Wie = > (Wan Wera), 1<i<n 
t=1 

(16’) Win = > A (Wis — Wire). 
14 


Since the operators W, are hermitian, eq. (16’) is simply the complex con- 
jugate of eq. (16). The stationary value of W is then 


ul n 1 n ; 
(17) W= Wt 5 > AW +5 D Wind. 
Si=1 


1=1 


We notice that at all stages of approximation W remains hermitian. As 
a consequence of eq. (10) all the operators W, are simultaneously diagonal (1°) 
in the configuration space of the two nucleons, and obviously all the W, 
with odd à vanish. The system (16) and eq. (17) become 


(18) Wald = > |r| Wier — Warua|r") <r" | Air ar", 
Et, 


LUE 
(19) 7 |W ltr Walt > | rl Arr | War) dr" -8.c.. 
2 i=1, 


i 


With the help of eq. (10) we obtain that 


<r" | Ai|r> = Alr')d(r'— r") 


(13) Strictly speaking the w,(r) are still operators in the spin space. Reasons of 
invariance restrict greatly the possible spin dependence (L. EISENBUD and E. WIGNER: 
Proc. Nat. Acad. Wash., 27, 281 (1941)) and it is easy to see that a representation 
may be found in which all the w,(r) are simultaneously diagonal. 
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and eq. (18) reduces to the following system of linear numerical equations: 


(18)' Walt) = D Loisir) — Wirualr) ar) . 1<i<n 


l=1 


This shows that the A,(r) are real and eq. (19) becomes 
(19) w(r) = w,(r) + > Wita(r)Ai(r) « 
i=1 


Since the odd w, vanish the expression (19') for w contains only even hie 
Let us denote with 2m the largest even number contained in n. An inspection 
of the system (18’) shows that by adding the equation with 1/=2j to the 
equation with | = 2j—1 one obtains the following m equations in the m 
even À, which are the only interesting ones (14): 


Ja(W_e— Ws) + Ag(wy — We) + + + Aam(Womt1 — Wamte) = — Ws 
(20) Aw, — We) m Aw — Ws) + + Am Want — Wom-+a) = lin 

Ao(Wom+1— Wamt2) + A Wamt+s — Womtr) + see + AomWamte — Wim) — Wome - 
Furthermore 
(21) w—=W-+ Aw; + u. + AymWom+2 - 


The order of approximation for w is evidently given by the number m of 
even À, introduced in the «trial functions >». 

At this point it is useful to introduce explieitly the dependence of all quan- 
tities on the coupling constant g. Let us define 


(22) WE Use 


From equations (20) (21) (22) we obtain 


m 
(m) 21 
NP ser 


en 


m 

(m) 21 
da: 41-19 
1=0 


(23) wg 


(1) We have for simplicity dropped the explicit dependence of our quantities on 7. 
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where the suffix (m) explicitly introduced stresses the order of approximation 
for w. The quantities P°” and Q/” are defined in the following way. Let us 
introduce the rectangular matrix Q(® with m columns and m-+1 rows with 
elements 


(24) Vij = Usitoi+4 1<i<cm, 1<j<m +1, 


We define Q{”’ as the determinant obtained by canceling in Q” the j-th row, 
taken with the sign (— 1)’*!. Furthermore we define 


m+s+2 


(25) PO = Son QO |: 
k 


=1 
The first approximation given by eq. (23) is 


: Us 

26 WEGE et re 5 

(26) w Put 9 DR 

It is clear that for large g? the difference from the Levy-Klein potential 
g’u, + gu, may be quite significant. 


4, — Discussion. 


Our method gives an approximation to w(r) by means of a sequence of 
fractions of polynomials in g?. It is worthwhile to notice that the approx- 
imation is not a ratio of two series, because all the coefficients Q’” and Pi” 
change at each order. It is well known that in this way it is possible to 
represent a much wider class of functions of g. As an example it is easy to 
show that the sequence converges for arbitrary values of g if w is supposed 
to be a meromorphic function of g (15). 

In any case the convergence of the sequence does not depend on the magni- 
tude of g, but only on the sequence 


Us 5 Ug 5 Uony + 


determined only by the nature of the interacting fields. 
For values of g for which the power series 


(27) Pu + gut... +9? an + oy 


(5) It is trivial that if the function is rational the exact value is already obtained 
by the approximation of an appropriate order. 


EL a ae 5) Wry L 
À À ; 
.. > 
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converges, the limit of our sequence coincides with the limit of the series (26). 
The fact that the present approximation, deduced from a variational prin- 
ciple, yields an analytic continuation of the power series (26) for values of g? 
outside of its radius of convergence, may be considered as an argument in 
favour of the reliability of the method. 

It may be noticed that our choice of the 4;, determined in order to make W' 
stationary is certain much better than the choice A, = Ay =... = dom = 1, 
leading to the Lévy-Klein potential. 

The application of the present method to scattering problems, including 
the contributions of radiative corrections is being studied. Preliminary results 
show that for the pion-nucleon scattering our first approximation gives results 
analogous to those obtained by solving the Tamm-Dancoff integral equation 
by means of a variational method. 


The authors wish to thank Prof. G. WATAGHIN for his kind interest in 


the work. 


RIASSUNTO 


Si fornisce un metodo di approssimazioni successive per il calcolo del potenziale 
adiabatico tra due nucleoni fissi, nel quale ad ogni stadio si tiene conto di un numero 
limitato di mesoni e di coppie. Il metodo, facendo uso di un principio variazionale nel 
quale si rende stazionario il potenziale adiabatico, fornisce un’approssimazione di questo 
come successione di frazioni di polinomi in g2, i cui coefficienti contengono solo i poten- 
ziali ealeolati da Lévy e Kirin con metodo perturbativo. La convergenza di questa 
suecessione non dipende dalla grandezza della costante di accoppiamento e il suo limite, 
se esiste, rappresenta il prolungamento analitico della serie di potenze in g? al difuori 
del suo cerchio di convergenza. 


110 — Il Nuovo Cimento. 


See 
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Eccitazione dei nuclei per assorbimento dei mesoni =. 


II. Risultati e conclusioni. 


V. DE SABBATA, E. MANARESI e G. PUPPI 


Istituto di Fisica dell Universita - Bologna 
Istituto Nazionale di Fisica Nucleare - Sezione di Padova 


(ricevuto il 21 Ottobre 1953) 


Riassunto. — Nel presente lavoro si completa una ricerca sull’assorbi- 
mento dei mesoni x, legati nell’orbita A da parte dei nuclei medio- 
pesanti (4 =100). Il problema viene trattato col metodo di Montecarlo, 
ed 1 risultati confrontati con l’esperienza. Distribuzione delle stelle in 
funzione del numero dei rami visibili, percentuale dei protoni veloci 
(E > 30 MeV), rapporto di emissione particelle «/protoni, percentuale 
delle catture senza rami visibili, sono in accordo con i dati sperimentali 
se si suppone che l’atto primario di assorbimento, in seno alla materia 
nucleare, sia compiuto da una coppia protone-neutrone. Ricordando i 
risultati della cattura dei mesoni x in volo da parte dei nuclei si potrebbe 
tentare di concludere che la cattura dei 7 in seno alla materia nucleare 
a tutte le energie (almeno fino a 100 MeV) avviene essenzialmente da parte 
di una coppia protone-neutrone. 


Introduzione. 


In un precedente lavoro (!) abbiamo discusso, nelle sue linee generali, il 


problema della eccitazione di un nucleo pesante delle lastre, a seguito della 
cattura di un mesone z~ legato nell’orbita piu profonda. Avevamo fatto la 
ipotesi semplice che il mesone fosse catturato da una coppia di nucleoni, cosicche 
il risultato netto dell’assorbimento si poteva riassumere nel fatto che da un 
punto interno del nucleo partivano due nucleoni veloci, all’incirca nella stessa 


direzione ma in versi opposti. 


Avevamo distinto due fasi essenziali del processo (2), delle quali una fase 


1) G. Puppı, V. DE SaBBarA e E. MANARESI: Nuovo Cimento. 9, 726 (1952). 
2) R. SERBER: Phys. Rev., 75, 1508 (1949). 
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istantanea in cui i nucleoni veloci, originati dall’assorbimento del mesone, 
attraversano la materia nucleare .generando per urto dei secondari veloci e 
costituendo cosi, ciascun nucleone primario coi suoi secondari una piccola 
valanghetta di nucleoni. Avevamo seguito il destino di questi nucleoni fino a 
che essi raggiungono la superficie del nucleo dove in parte vengono riflessi, 
e in parte fuoriescono. I nucleoni emessi in questa fase sono generalmente 
detti nucleoni di knock-on, e trasportano con sé circa 2/3 della energia globale 
di eccitazione. L’energia delle particelle che non escono immediatamente, si 
pud immaginare che venga distribuita rapidamente a tutti i componenti del 
nucleo ed emessa nella seconda fase del processo, più propriamente detta di 
evaporazione. 

Avevamo discusso le ragioni teoriche che giustificano la schematizzazione 
del processo globale in due fasi. Un esame teorico (?) della fase transitoria tra 
le due ha sostanzialmente confermato l’assunto. 

Nel presente lavoro discuteremo i dettagli del calcolo delle due fasi del 
processo, per entrambe condotto col metodo di Montecarlo, e dal confronto 
con i dati sperimentali disponibili vedremo quali conclusioni e riprove si pos- 
sono trarre sul modello impiegato. 


PARTE I. 


FASE ISTANTANEA DEL PROCESSO 


1. — Punto di eattura del mesone rt”. 


Abbiamo supposto, per questo processo, che la cattura avvenga mentre il 
mesone si trova in uno stato legato, e precisamente nello stato più profondo, 
e che le forze responsabili del legame siano solo le forze coulombiane. Questo 
significa, in altre parole, che consideriamo la perturbazione dovuta alle forze 
nucleari, responsabili dell’assorbimento del mesone, come piccola, come si pud 
vedere in realtà. Abbiamo pertanto considerato il mesone, legato da forze 
coulombiane come una particella e lo abbiamo descritto quindi con una equa- 
zione di Schrédinger ordinaria. 

Come potenziale si è preso, come di consueto, il potenziale di Thomson 


per l’interno 
y Ze?|\3 1 Ne 
i R |2 2 Re 


dove con R si é indicato il raggio del nucleo. 


(3) A. Kinp: Nuovo Cimento, 10, 176 (1953). 
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Impiegare il modello di Thomson, che alla superficie del nucleo si raccorda 
col potenziale coulombiano puro con continuità, significa naturalmente consi- 
derare che la densità di carica positiva nel nucleo sia costante, e per di più 
che improvvisamente tale densità alla superficie del nucleo cada al valore 
zero. In altre parole la densità della carica positiva, e con essa la densità dei 
protoni, e con essa, ancora più in generale, la densita della materia nucleare 
ha, in questo modello, una discontinuita alla superficie del nucleo per cui passa 
da un valore zero ad un valore costante. BE evidente che la unica ragione che 
consiglia l’uso di tale modello & la sua semplicita, e un poco forse il fatto che 
per cambiarlo si desidererebbe poterlo fare in base ad una concreta conoscenza 
dell’andamento della densita dei protoni e dei neutroni nel nucleo, su cui si 
hanno invece soltanto delle indicazioni teoriche (?). 

Il problema ha un notevole interesse soprattutto per la definizione della 
barriera coulombiana presentata dal nucleo quando una particella carica 
tenta di entrarvi o di uscirvi. Quello che si puö dire, dal confronto dei risul- 
tati ottenuti con distribuzioni continue di densità di carica positiva, è che 
il modello di Thomson è un modello estremo e che le barriere calcolate con 
distribuzioni continue sono tutte inferiori a quelle calcolate col Thomson, a 
parita di altri parametri. Per lo scopo del presente lavoro non abbiamo rite- 
nuto necessario allontanarci dal modello convenzionale, pur sapendo di trat- 
tare un caso ideale, le cui previsioni, perd, non si discostano tanto da quelle 
degli altri modelli da alterare le conclusioni. D’altra parte, l’uso del modello 
di Thomson era richiesto anche da ragioni di coerenza, avendo poi descritto 
Vinsieme dei nucleoni, nei riguardi delle forze nucleari, come un gas di Fermi 
racchiuso in una buca di potenziale di raggio pari al raggio nucleare, ed essendo 
sempre opportuno l’uso di un modello definito agli effetti della critica dei 
risultati. Con le ipotesi fatte dunque, la autofunzione del mesone legato nello 
stato fondamentale si pud scrivere come 


'uZe: 


1 exp [—  ar°] 
7.30; Hrn RR: 


Hola) — 2 a = | 


a R (u = massa del mesone), 


ed & praticamente contenuta entro il nucleo. Infatti, l’autofunzione del mesone 
si puö ricondurre a quella di un oscillatore nello stato fondamentale, la cui 
escursione, classica, à circa 1/3 del raggio nucleare. Poichè nel nostro problema 
quello che interessa è trovare le probabilitä relative di cattura delle varie zone 
del nucleo, la autofunzione globale deve essere normalizzata sul volume del 
nucleo e cioè 


R 
An [ : 
a) exp or dr. 
0 


(*) L. ROSENFELD: Nuclear Forces (Amsterdam, 1948), p. 242. 
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Data la simmetria del problema è conveniente suddividere il nucleo in tante 
intercapedini individuate dai raggi delle superficie sferiche concentriche che le 
racchiudono. Se diciamo N il numero delle intercapedini, e vogliamo che a 
tutte le intercapedini corrisponda la stessa probabilità di presenza (e quindi 
la stessa probabilità di cattura), dovremo porre la condizione, indicando con ”, 
il raggio della n-ma superficie sferica, 


Tnt 


sa | Rire dr = — 


N° n= 0, 1, 2, 35...) 9s 


Tn 


Noi abbiamo preso N— 10 ed abbiamo scelto a rappresentare i punti di 
cattura per ogni intercapedine un valore di r intermedio tra i raggi delle super- 
ficie sferiche che la racchiudono. Nella tabella I sono riportati i 10 valori che 
individuano le distanze del punto di cattura dal centro del nucleo. 


TABELLA I. 


an 0 ; 1 2 | 3. 


a a CARE IT NES 


0,746 | 0,804 | 0,861 0,917 


| | 


(7/R) | 0,303 | 0,439 | 0,536 


| 
0,615 | 0,684 


2. — Modalità della eattura. 


Poichè il modello nucleare da noi tentato è, come gia accennato, quello 
del gas di Fermi, e poichè il minimo numero di nucleoni che possono inter- 
venire nel processo di cattura del mesone, per la conservazione dell’energia 
e dell’impulso è due, abbiamo trattato due casi distinti in cui il mesone & assor- 
bito da una coppia protone-neutrone (P-N) oppure da una coppia protone- 
protone (P-P); naturalmente, nel modello adottato la probabilità che due 
nucleoni si trovino in un certo elemento di volume, entro il raggio di azione 
delle forze nucleari, & retta dalla legge del caso, e cosi si pud parimenti dire 
a proposito dell’impulso degli stessi nucleoni. Cosi, la probabilita che uno dei 
due nucleoni che assorbono abbia un certo impulso dipende unicamente dalla 
distribuzione degli impulsi in un gas di Fermi degenerato. Con questo si viene 
a dire che la probabilita che due nucleoni, correlati spazialmente, assorbano 
il mesone non dipende dallo stato di moto dei due nucleoni. Comunque, per 
trovare le caratteristiche della distribuzione degli impulsi bisogna calcolare 
Vimpulso massimo per i neutroni e per i protoni dati dalle relazioni ben note (°) 

(920)2/*h (2.N\1/* (920) "8h [2Z\13 
en, na OOM 


27 


A 


(5) L. ROSENFELD: Nuclear Forces (Amsterdam, 1948), p. 192. 
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A conti fatti, per questi impulsi massimi si ottengono rispettivamente delle 
energie massime di 26 MeV e 22 MeV circa. Assumendo come energia di legame 
dell’ultimo nucleone il consueto 8 MeV, risulta che le buche di potenziale in 
cui debbono essere rinchiusi questi due gas degeneri, rispettivamente di neutroni 
e di protoni, hanno profondità di 34 MeV e di 30 MeV. Questo modello nucleare 
a due buche di differente profondità porta come conseguenza che i neutroni 
sono leggermente più caldi dei protoni, ma la differenza è molto modesta. La 
probabilita che un nucleone abbia un impulso di modulo compreso tra p e 
p+dp e direzione e verso compresi nell’angolo solido infinitesimo dQ rispetto 
ad un riferimento arbitrario € pertanto 


L’espressione scritta è identica a quella dell’elemento di volume infinite- 
simo di una sfera unitaria; si pud cosi pensare di rappresentare l’insieme degli 
impulsi in una sfera in modo tale che il rapporto tra un volumetto elementare 
e il volume totale misuri anche la probabilita che Vestremo di un vettore 
impulso spiccato dal centro cada nel volumetto in parola. Noi abbiamo diviso 
la sfera di Fermi in tante cellette (10000) a mezzo delle seguenti condizioni: 


Pitt Orta Pr+1 


t es, 1 
a | Pape 5 | sin bad = = | a= 


Pr 9% Pr 


Con N, = 10, No = 20, N, = 40 Nea SI, 102000: 


Il fatto che l’estremo di un vettore impulso cada in una o in un’altra cel- 
letta @ regolato puramente dal caso e le cellette sono dello stesso volume e 
pertanto rappresentano casi equiprobabili. 

Supponiamo ora di conoscere gli impulsi dei due nucieoni che assorbono 
e diciamoli rispettivamente p, e p,. Dopo aver assorbito il mesone, ed esser- 
sene divisa la energia totale, che, considerando la energia cinetica trascurabile, 
si pud confondere con la sua energia di riposo, i due nucleoni ripartono con 
due nuovi impulsi p, e p, che in approssimazione non relativistica si possono 
calcolare a partire dalle ovvie relazioni di conservazione seguenti (M— massa 
del nucleone): 


! 


P, +P:= Pi tps, » EP —=p LP +2Me. 


E opportuno introdurre un vettore ausiliario K legato agli impulsi dalle 
relazioni: 


Pp=mEK, pp -K: 
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Se diciamo #, e #, gli angoli formati dai vettori p, e p, con la direzione 
orientata di K, avremo ancora: 


pe = pi + K? + 2p, K cos 0, 3 
Dy = pe + K*— 2p.K cos d,, 


e quindi in definitiva: 
Ke + 


+ K(p, cos 3, — D: 608 d,)= Muc?. 


Come si vede nella scelta di K 
e’& della arbitrarietä e si pud farne 
Voggetto di una estrazione a caso, 
come diremo più avanti. Noti P, 
e p, possiamo considerare con- 
eluso il processo di cattura, e 


E (MeV) 
40 60 80 100 120 140 160 


dobbiamo oceuparei di quello che 


accade mentre questi due nu- Fig. 1. - Distribuzione energetica dei nucleoni 
eleoni veloci, originati dalla cat- che hanno catturato il nr. 

tura, attraversano la materia nu- 

cleare, dal punto di cattura fino alla periferia del nucleo. La fig. 1 mostra 
la distribuzione energetica in seno alla materia nucleare dei nucleoni originati 
dall’assorbimento, indipendentemente dalla loro carica; per confronto è ripor- 
tata anche la distribuzione energetica della materia nucleare ordinaria, cioe 
prima dell’assorbimento. 


3. — Eccitazione del nucleo. 


Come gia spiegato, dal punto di cattura partono due neutroni, oppure un 
neutrone e un protone, a seconda del tipo di coppia che ha assorbito, di dire- 
zione, verso ed energia nota. Attraversando la materia nucleare, essi interagi- 
ranno con questa scambiando energia, € precisamente cedendo energia, dato 
che la materia nucleare incontrata à fredda, cioè i nucleoni incontrati occu- 
pano i pit bassi stati di traslazione del gas di Fermi. Coerentemente col nostro 
modello, e come gia ampiamente discusso in (1), tratteremo le interazioni con 
la materia nucleare come una successione di urti tra due particelle le cui moda- 
lita sono regolate dalle ben note sezioni d’urto protone-neutrone e protone- 
protone, ed ammetteremo che la sezione d’urto neutrone-neutrone sia identica 
a quella protone-protone. In fig. 2 sono riportate le sezioni d’urto integrali 
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che sono state impiegate per questo lavoro (°). La sezione d’urto differenziale 
nel centro di massa protone-protone e neutrone-neutrone & stata considerata 
isotropa a tutte le energie che interessano, e la sezione d’urto protone-neu- 
trone (?), invece, di andamento 
mbarn . : . 
| parabolico in cos d, e precisa- 
300 
mente del tipo 


280 
al od =a + b cos? À ; 


| a e b risultano funzioni della 
| energia e si possono facilmente 
200 ed 
| determinare dal valore delle 
180 es s 5 B 
sezioni d’urto sperimentali a 
90° e a 1800. Fissati questi 


749 valori a tre energie, come ri- 
120 N-P portato in tabella II, per ener- 
00 À gie intermedie si è interpolato. 
sol Fissiamo ora l’attenzione su 
cel es uno dei nucleoni veloci che 
N partono dalla cattura; la prima 

cosa da fare sarà quella di de- 
is E(MeV) finirne il cammino libero medio 


etl D at eae Oe ee ee 
GTS" CORNE EE ee in materia nucleare, e questo 


si pud fare in base alla cono- 
scenza delle sezioni d’urto e 
della densité della materia nucleare. Vi è ancora da dire, perd, che le se- 
zioni d’urto riportate si riferiscono ad urto di particelle veloci contro parti- 
celle ferme, mentre nel caso presente le particelle incontrate sono in moto. 


Fig. 2. — Sezioni d’urto totali. 


TABELLA II. 


E - (MeV) 40 90 260 | 
a (mb/sterad) 11 2,8 0,5 
b  (mby/sterad) 8 9,5 7,5 


Poichè, perd, nell’intervallo energetico che ci interessa le sezioni d’urto inte- 
grali hanno approssimativamente una dipendenza dalla energia del tipo 1/E, 


(5) J. DE JUREN e B. J. Moyer: Phys. Rev., 81, 919 (1951); R. W. BIRGE: Phys. 


Rev., 80, 490 (1950); R. W. BIRGE, U. E. Kruse e N. F. Ramsey: Phys. Rev., 83, 
274 (1951), 


(7) E. Kerry, C. Lerru, E. SeGr e C. WIRGAND: Phys. Rev., 79, 96 (1950). 
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le sezioni d’urto mediate sulla sfera di Fermi forniscono, in approssimazione 
non relativistica, lo stesso valore che si ottiene non tenendo conto del moto 
dei nucleoni urtati. Dette pertanto o,, € o,, le dette sezioni d’urto integrali, 
possiamo definire il cammino libero medio in materia nucleare nel modo con- 
venzionale : 


1 1 1 1 il 


= = 0x0xx + OpOnr : — — +-— = QpOpp + OnOpn - 
Ay Axe Ayn As Ae Roe 


Agli effetti del Montecarlo è opportuno suddividere il percorso del nucleone 
in tanti intervalli di lunghezze tali per cui in ciascuno di essi si abbia la stessa 
probabilità a priori di una collisione. Si vede immediatamente che, detto N 
il numero totale degli intervalli, ed essendo exp [— x/A] la probabilita che 
la collisione avvenga dopo un percorso %, la detta condizione di equiproba- 
bilità si pud serivere nella forma: 


exp {— %_1/4] — exp [— À] = 1/N - 


Di qui si trae che la suddivisione & definita da certi punti di ascisse Xp, 
a partire dal punto di cattura, calcolabili dalla relazione 


“u 
JA 


din m: 

Ad individuare ogni intervallo, abbiamo al solito assunto un punto inter- 
medio tra gli estremi. Questo serve a calcolare il punto in cui avviene la col- 
lisione; per sapere poi se la collisione del neutrone avviene contro un protone 
o contro un neutrone basta ricorrere ai rapporti Ales Ant Aww» ed analoga- 
mente per un protone urtante. 

Supponiamo ora, di fissare che in un certo punto della sua traiettoria il 
nucleone proveniente dalla cattura, ad esempio un neutrone, ha urtato contro 
un protone. Bisogna sapere ora quale & l’impulso del nucleone urtato, passare 
al sistema di riferimento del centro di massa, fare una rotazione a caso dei 
vettori rappresentativi degli impulsi nel centro di massa, e ritornare nel si- 
stema di riferimento iniziale. A questo scopo bisogna, come al solito, suddivi- 
dere la sezione d’urto differenziale nel centro di massa in intervalii equipro- 
babili, che poi andranno sorteggiati, e questo si pud realizzare con la condizione 


Pets 
Qn | 1 
= o(9) sind db = —. 
o N 


Dy 


Quando si ritorna, dopo avvenuta la diffusione nel centro di massa, al 
sistema di riferimento iniziale, pud accadere che uno degli impulsi finali cada 
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dentro la sfera di Fermi, cioè si riferisca ad uno stato di traslazione gia occu- 
pato; poiché questo & proibito dal principio di esclusione, in questo caso do- 
vremo considerare che l’urto non sia avvenuto, perchè proibito, e far conti- 
nuare al nucleone veloce la sua corsa. Ancora una osservazione a proposito 
della scelta dell’impulso del nucleone urtato: Vapprossimazione da noi usata 
è stata quella di sceglierlo tra uno dei possibili valori all’interno della sfera di 
Fermi, mentre in realta, dato che la sezione d’urto contiene il fattore 
(Mine = Port) ’ 

sarebbe pit giusto sceglierlo in una distribuzione deformata, essendo relativa- 
mente più probabile l’urto contro un nucleone di verso contrario che contro 
uno di verso concorde. Sta di fatto, perd, che la dipendenza della sezione d’urto 
dalla energia è tale da compensare largamente questo fattore, cosi che la nostra 
approssimazione, fino a che è valida approssimativamente la legge 1/H per 
le sezioni d’urto integrali, pud ritenersi soddisfacente. 

Se l’urto € permesso, da esso hanno origine due nucleoni più veloci dei 
nucleoni che formano la materia nucleare ordinaria e di cui bisogna seguire 
Vulteriore destino attraverso i successivi urti, che si trattano con le modalita 
specificate dianzi: da ogni nucleone proveniente dalla cattura ha quindi ori- 
gine in materia nucleare una valanghetta di nucleoni che bisogna seguire fino 
alla superficie del nucleo. Qui giunti, i nucleoni che compongono la valanghetta 
subiscono varia sorte. I neutroni che hanno energia superiore si alla energia 
massima della sfera di Fermi, ma per una quantita inferiore ad 8 MeV ven- 
gono senz’altro riflessi, i protoni, a causa della barriera coulombiana, vengono 
praticamente sempre riflessi se l’energia non supera il valore massimo della 
sfera di Fermi di 13 MeV. 

Nello svolgimento del Montecarlo, quando un nucleone cadeva in questo 
intervallo energetico veniva abbandonato e considerato catturato definitiva- 
mente, Senza seguirne il destino fin proprio alla superficie del nucleo. I nucleoni 
invece che giungono alla superficie del nucleo con energia cinetica superiore, 
rispettivamente a 34 MeV per i neutroni e a 35 MeV per i protoni hanno una 
certa probabilita di essere riflessi ed una certa probabilita di fuoriuscire con 
energia ridotta, pari alla somma algebrica della energia cinetica e della energia 
potenziale. Per i neutroni abbiamo calcolato la probabilita di trasmissione per 


il cambio di lunghezza d’onda tra l’interno e l’esterno del nucleo, con la for- 
mola ben nota 


(1) p = Sh 


Ee Cee Eee 


ove k e K sono rispettivamente i numeri d’onda all’interno e all’esterno del 


a: are Gt. r FES À 5 = \ N 
u. ene re EN MED EE \ 
= 3 - = - ; 
‘ ; 
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nucleo e legati alla energia dalle relazioni: 
i 
b=> V2ME, K= 


(dove — V & la profondita della buca). 
Per i protoni bisogna tener conto della barriera coulombiana, e percid 
abbiamo impiegato per il coefficiente di trasmissione la formola composta (°) 
4kK . 22e? 


2 T = =: exp | — 
(2) = pm | Fe, 


(2W — sin 2W) 


5 


ove il primo fattore rappresenta il coefficiente di trasmissione dovuto al cambio 
di lunghezza d’onda, come per i neutroni, mentre il secondo rappresenta la 
probabilità di passaggio della barriera, sal- 
colato col metodo di W.K.B.. Al solito i 
simboli hanno il significato convenzionalee — 002 Fy = 


—R = 


1 
' 
! 
1 


W rappresenta il rapporto R;/R, tra il rag- 
gio interno e il raggio esterno della bar- 


—— 30 MeV 


riera, come segnato in fig. 3). Tutti i pa- 

rametri sono funzioni della energia e i ri- TRE an 

sultati sono stati tabulati e messi in grafico. Fie. 3. — Buche di potenziale. 
Abbiamo ora tutti gli elementi per | 

descrivere come si & in effetti proceduto col metodo di Montecarlo. 


4. — Organizzazione del Montecarlo. 


L’uso del metodo di Montecarlo, nei più svariati problemi, si € ‘andato 
diffondendo negli ultimi anni, come lo testimonia la gid rieca letteratura in 
proposito. Esso consta sostanzialmente nella costruzione di un certo numero 
di possibili esemplari del fenomeno che interessa, quando questo fenomeno 
sia il risultato di una successione, non troppo numerosa, di processi elementari, 
perfettamente noti ed indipendenti l’uno dall’altro. Il numero degli esemplari 
deve essere sufficientemente elevato da permettere di rilevarne i risultati stati- 
stici, che debbono essere confrontati con i risultati statistici di un adeguato 
numero di casi effettivi. Il fenomeno in esame, della eccitazione dei nuclei per 
assorbimento di un mesone # si presta a questo trattamento, purche non 
vi sia luogo a dubitare che il modello impiegato cada in difetto. Il modello, 
come abbiamo già accennato, riduce le interazioni in materia nucleare ad una 


(8) I. M. Buarr e V. F. WEISSKOPF: Theoretical Nuclear Physies (1952); G. BER- 
NARDINI, E. T. BootH e $. J. LINDENBAUM: Phys. Rev., 83, 669 (1951). 
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successione di urti nucleone-nucleone, trattati come gli urti tra due nucleoni 
isolati (°), salvo le conseguenze del principio di Pauli. Questo € solo un modello 
approssimato, perché, certamente, quando si considerano nucleoni di energia 


Fig. 4. 


relativamente bassa, le forze a piu corpi intervengono sensibilmente, ma à 
purtuttavia il migliore che si possa attualmente usare ed il più facile a trat- 
tare; vediamo dunque come lo si è organizzato. Il primo quesito che si pone 
è il punto di cattura del mesone x; a tale scopo, come gia detto, il nucleo 
è stato diviso in 10 intercapedini sferiche tali che in ciascuna di esse la proba- 
an: bilità a priori di cattura del mesone sia 1/10. Si associa quindi a ciascuna inter- 

:apedine un numero dal 0 al 9, i 10 numeri si mettono in un’urna e se ne 


(®) M. L. GOLDBERG: Phys. Rev., 74, 1269 (1948). 


i. eo 
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estrae uno a caso. In realtà sono state impiegate tavole di numeri estratti a 
caso (1%); a partire da un punto arbitrario di queste tabelle, fissiamo che la 
prima cifra dica in quale intercapedine è stato eatturato il mesone. Estratto 
questo e supponendo di trattare l’assorbimento da parte di coppia P-N, il 
secondo quesito che si pone & il sapere gli impulsi dei due nucleoni che assor- 


+ AN 
AE À... AN 
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2 ' C2) pea) 
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Fig. 5 a). — Cattura da parte di P-N. Fig. 5 b). — Cattura da parte di P-P. 
—— Spettro dei neutroni e dei protoni di knock-on.  ----- Spettro comprensivo 


dell’evaporazione. 


biranno il mesone. Per questo bisogna estrarre à sorte la celletta della sfera 
di Fermi in cui cade l’estremo del vettore impulso del protone, ed analoga- 
mente del neutrone. Poichè avevamo suddiviso il modulo dell’impulso in 10 inter- 
valli, la seconda cifra a caso dard il modulo dell’impulso del protone; per 
Vorientamento, dato che avevamo suddiviso l’angolo 9 in 25 parti e l’angolo y 
in 40, occorre estrarre due cifre per 9 e due cifre per y. La identica operazione 
si fa per il neutrone; gli orientamenti si riferiscono ad un asse polare princi- 
pale che coincide con un diametro del nucleo, e ad un piano passante per 


(9) M. G. KenpaLL: Tables of random sampling numbers (Cambridge, 1951). 


ES 
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questo asse; ma per questa prima fase abbiamo trovato più conveniente rife- 
rire p, e p, al vettore ausiliario K, ed estrarre la orientazione di K rispetto 
all’asse polare principale e costruire le risultanti, dopo l’assorbimento, p, e p,- 
Questo ci permetteva di separare la prima fase del processo dal resto e di 
partire addirittura con la condizione dopo l’assorbimento. Per lo svolgimento 
ci siamo valsi di due modelli meccanici: uno, fig. 4, rappresentante il nucleo, 
consisteva di un’asta verticale di lunghezza di 2 metri corrispondente al detto 
diametro del nucleo. La meta superiore di questa asta era suddivisa in 10 tratti 
rappresentanti i punti intermedi dei 10 intervalli equiprobabili di cattura. 
Uno scorrevole veniva fissato nel punto di cattura e su esso erano imper- 
niate due aste orientabili a piacere che rappresentavano le direzioni dei 
due nucleoni veloci dopo la cattura. La modalita della cattura era invece risolta 
con un altro modellino meccanico a parte, che non descriveremo per brevita. 
Una volta sistemate nel modello le direzioni di moto dei due nucleoni si comin- 
ciava a seguirne uno e si estraeva a caso con la tecnica ormai nota il punto 
di collisione con un nucleone del nucleo e le direzioni di moto del nucleone 
urtante e di quello urtato dopo Vurto. Le traiettorie di queste asticelle veni- 
vano visualizzate in scala, essendo le aste, a cannocchiale, di lunghezza varia- 
bile, e provviste alle estremita di cerniere e goniometri. Era cosi possibile passo 
per passo visualizzare la traiettoria di tutti i nucleoni che entravano a far 
parte della valanghetta, fino alla periferia del nucleo. Quando accadeva che 
N un nucleone passasse, senza urtare ul- 


AE teriormente, oltre la superficie del nu- 
1,6 


cleo, allora si estraeva ancora à sorte, 
12 in base alla conoscenza del coeffi- 
ciente di trasmissione T e di quello 
i. di riflessione (1— T) per quella data 
14 energia € per quel certo tipo di par- 
rime ticella, se effettivamente la particella 

2 40 60 6 100 120 140 


0 usciva o veniva riflessa. Un esemplare 
Fig. 6. — Energia residua. è completo quando di tutte le parti- 
celle costituenti le due valanghette se 
ne € seguita la sorte fino alla superficie, ei dati che si hanno a disposizione 
dopo questa prima fase del processo sono: tipo ed energia delle particelle che 
escono dal nucleo in questa fase istantanea, tipo del nucleo residuo e stato di 
eccitazione del nucleo. Questi risultati sono raccolti, per 100 catture PN e 
100 catture PP nelle figure 5a, 5b, 6). Anche la distribuzione angolare delle 
particelle emesse istantaneamente si pud ottenere da questi risultati, ma in 
questo lavoro non vi siamo particolarmente interessati. I risultati qui ri- 
portati sono sostanzialmente identiei a quelli riferiti in I e basati su un nu- 
mero di esemplari all’incirea meta; questa stabilità ci conforta a pensare che 
la statistica che presentiamo sia sufficientemente indicativa. 
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PARTE 1. 


FASE EVAPORATIVA DEL PROCESSO 


1. — Nucleo residuo. 


Come conclusione della prima fase del processo, abbiamo, come gia detto, 
un gruppo di particelle uscite istantaneamente ed un gruppo di particelle riflesse 
alla superficie del nucleo, che rinunciamo a seguire nella loro ulteriore evolu- 
zione entro il nucleo per le ragioni discusse in I, e che in breve tempo tendono 
a realizzare approssimativamente entro il nucleo una condizione di quasi equi- 
librio. II nucleo cosi eccitato sara trattato con i metodi della teoria evapora- 
tiva, ed i risultati globali confrontati coi risultati sperimentali. Nella prima fase 
del processo abbiamo trattato un nucleo ipotetico con A= 100 e Z— 43, che 
dovrebbe avere proprietà intermedie tra quelle dei nuclei pesanti (Br, Ag, I) 
delle lastre fotografiche. Avremmo dovuto, in realta trattare separatamente 
i due casi del Br e del Ag (I entra solo in piccola percentuale) dato che questi 
nuclei hanno raggi sensibilmente differenti e precisamente nel rapporto 0,9 
ed anche barriere coulombiane sensibilmente differenti e precisamente nel 
rapporto 0,82. Tenuto conto, perd, che nella prima fase del processo intervengono 
solo proprieta che sono funzioni continue e monotone del numero di massa, 
la scelta di un nucleo intermedio non pud alterare i risultati e si risolve in 
una notevole semplificazione del metodo, di per sé gia abbastanza gravoso e 
lungo. Nella seconda fase del processo, invece, abbiamo pensato che non fosse piu 
lecito ricorrere alla detta approssimazione ed abbiamo separato i processi iniziati 
nell’Ag da quelli iniziati nel Br. Abbiamo ripetuto per ogni nucleo residuo il 
processo di evaporazione, considerando il nucleo iniziale una volta Ag e una 
volta Br e stabilendo con una estrazione quale degli isotopi si dovesse per 
ogni caso, scegliere, conoscendo l’abbondanza isotopica. Caso per caso, dalla 
conoscenza del numero e tipo delle particelle emesse nella fase istaneanea e 
del tipo di nucleo che ha assorbito, & possibile immediatamente stabilire quale 
sia il particolare nucleo eccitato che inizia la evaporazione. Abbiamo ritenuto 
opportuno in questa seconda fase impiegare questi criteri per il fatto che le 
probabilità relative di evaporazione di un tipo piuttosto che un altro di parti- 
cella (neutrone, protone, alfa) dipendono in maniera sensibile dal particolare 
tipo di nucleo che evapora, dal suo numero di massa, dal suo numero isotopico, 
e non sono senz’altro funzioni monotone del numero di massa. Non ci nascon- 
diamo che l’aver trascurato la differenza tra Ag e Br nella prima fase, ed aver 
usato sempre uno stesso nucleo intermedio, ha alterato leggermente i risultati, 
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nel senso che se Ag e Br fossero stati tenuti distinti, all’Ag sarebbe in media 
rimasta un poco più di energia residua e al Br un poco meno che non all’ipo- 
tetico nucleo con A=100. Ma poiché alla fine i risultati vengono conglobati, 
& nostra impressione che gli eventuali scarti non siano tali da non essere a priori 
conglobabili nella stessa indefinizione dei parametri caratteristici dei nuclei, 
come il raggio e la barriera coulombiana, e la funzione di densita della materia 


nucleare, ecc.. 


2. — Processo evaporativo. 


Per trattare il processo evaporativo si € fatto uso del modello termodina- 
mico di Weisskopf, supponendo cioé che il nucleo eccitato emetta particelle 
il cui spettro energetico sia quello che deriva da una distribuzione energetica 
di un gas di Fermi a temperatura diversa da zero. Questa distribuzione energe- 
tica & data dalla 


(3) I(e) de = cost Ageo.(e)w(E — e— B) de, 


dove E & l’energia di eccitazione, e quella della particella emessa, B l’energia 
di legame della particella evaporata; E— e— B rappresenta quindi l’energia 
residua d’eccitazione del nucleo d’arrivo; o.(e) è la sezione d’urto per il pro- 
cesso a ritroso e g è una costante che da il numero degli stati di spin molti- 
plicato per la massa della particella evaporata ; œ@(E — e— B) rappresenta la: 
densità dei livelli energetici del nucleo di arrivo per intervallo unitario di 


energia: 
(4) o(E— e— B) = cost exp[2?Va(E— e— B)], 


a è una costante determinata empiricamente in maniera che la (4) riproduca 
le nostre conoscenze sulla densitä dei livelli energetici. 

Noi abbiamo limitato le nostre considerazioni alla evaporazione di tre tipi 
di particelle: neutroni, protoni e particelle « perchè per altri raggruppamenti 
le probabilitä di emissione sono trascurabili. Se fosse valida l’approssimazione 
geometrica classica, si dovrebbe prendere per i neutroni 


oe) che: 


Nel caso dei protoni e delle particelle « si dovrebbe poi tener conto del 
potenziale coulombiano e i calcoli, sempre in approssimazione classica, darebbero 


[ xR(1— V(R)le) per e> V(R), 
note | © per e< V(R), 
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dove V(R) rappresenta Valtezza della barriera. In tal caso la (3) diventerebbe 
(5) I(e) de = cost Agle— V(R)] exp [2Va(H — e— B)] de. 


Nella (5) perd la barriera di potenziale figura come un ostacolo classico; 
in realtà bisognera tener conto del fenomeno della trasparenza della barriera. 
Aleuni autori (1) hanno tenuto conto di questo fatto riducendo opportuna- 
mente il valore dell’altezza della barriera V(R) a un valore V'(R) < V(R). 
Ad esempio, per un nucleo standard A=100 e Z=43 per cui V,(R)— 9 MeV 


e V,(R)~18 MeV, LE COUTEUR prende Vi(R) = 5 MeV e V,(Rk) = 10 MeV. 


Noi, anziché fare uso della (5), abbiamo tenuto conto della barriera moltipli- 
cando la (3) per il coefficiente di trasparenza della barriera; nel caso poi dei 
neutroni, la (3) si & moltiplicata per il coefficiente di trasmissione. La (1) percid 
diviene 
(6) I,(e)de = AT g,e exp [2Va(B— e— B)]de (i = n, p, à), 
con T,, T, e T, date dalla (1) (2). 

In questa seconda fase del processo, non ei siamo accontentati di partire 
da un nucleo intermedio fra quelli delle lastre, ma abbiamo fatto il calcolo 


esplicitamente per ogni nucleo residuo dal processo di knock-on e quindi per 
ogni nuovo nucleo che via via si fosse presentato in seguito alla evaporazione ; 
e questo per tener conto dell’energia dello stato fondamentale del nucleo residuo, 


poiché evidentemente le probabilita relative di evaporazione dei neutroni, 
protoni e particelle « dipendono dalla energia dello stato fondamentale dei 
diversi nuclei di arrivo. Per questo, secondo il LE COUTEUR si é usata una 
formula semi-empirica che ci da le energie degli stati fondamentali: 


Ed) = — aA + BA — 8.) ; 


9, & il valore di @ per il nucleo più stabile N,, Z, di massa A ed è molto bene 
rappresentato dalla formula empirica 


A2l3 
Yee Tate 181)" 
mentre Ÿ = (N — Z)/A. 
x e B sono delle costanti determinate empiricamente. 
Si ha 
0H, (N, Z) 
ON 


OH, (N, Z 
AE = a + 26,(0— 8.) 


= — a + 28, (0—*#-), 


(2) K. J. Le CoUTEUR: Proc. Phys. Soc. London, A 68, 259 (1950). 
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con 
CM EEE 
pu = Bi aA nD ): 
dd, -9 9 
= pl1 44344 - ): 


Se allora un nucleo a con una energia d’eceitazione E, emette un fram- 
mento (n,2) lasciando il nucleo residuo b con una energia di eceitazione E,, 
la conservazione dell’energia ei dice che dovrà essere 


E,+E°=E#H,+E —w+e, 


dove E® e EH sono le energie degli stati fondamentali rispettivamente del 
nucleo a e b; e & l’energia cinetica del frammento emesso e w la sua energia 
di legame (da non confondere con la energia di legame del frammento nel 
nucleo). 
Cioé: 
EH 0e", 


avendo posto 


Q = EY — EO — w = a(n +2) 2(9 — d,)(nßrn — 2B») — W ; 


Q rappresenta cioè l’energia minima necessaria per l’emissione dal nucleo di 
un frammento (n, 2) con energia cinetica nulla. Percid la (6) andra ancora 
scritta alla stessa maniera solo che al posto di B si dovrà mettere Q che è 
comprensivo delle differenze delle energie degli stati fondamentali del nucleo 
di partenza e di arrivo. 


3. — Organizzazione del Montecarlo. 


Il metodo di Montecarlo é stato applicato nella solita maniera. La condizione 
CE 
A | >: 9:7; & exp [2Va(E en Q)] de — 1, 


0 


mn 
1 
— 


da le probabilità relative di evaporazione dei neutroni, protoni e particelle « 
da un dato nucleo. Si & quindi suddiviso l’intervallo energetico (0, H— Q:) 
in N parti equiprobabili, tali cioè da aversi 


En+1 


fe it 
ECTS 


En 
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questo & stato fatto per vari valori delle energie di eceitazione Æ del nucleo 
e per i tre tipi di particelle dette, riportando i punti di divisione in funzione 
di #; naturalmente tale preparazione € stata fatta per i diversi nuclei che si 
sono presentati nell’evaporazione. Infine si ¢ proceduto alle estrazioni: il primo 
numero estratto, o letto sulle tabelle, da il tipo di particella che evapora 
(n, p, ~) mentre il secondo dä l’energia cinetica di questa. Non si è tenuto 
conto dell’evaporazione di altre particelle perché la probabilità risulta molto 
piccola. Questo si pud vedere considerando il valore di Q. 
Avendo preso 


p = 23 MeV x = 8,6 MeV 


Bn = 0,88 By = 1,26 


questi ultimi due d’accordo con il valore di dd,/dA = 0,0009 corrispondente 
a un nucleo con A= 80, si hanno subito, per il termine a(n+z)—w che 
compare in Q, i seguenti valori 


n | p | = | ne | 


(8) 


| a(n+z)—w | 8,6 MeV | 8,6 MeV | 6,2 MeV | 15,1 MeV 


Come si vede dalla (8), per i deutoni il termine 15 MeV che compare nello 
exp della (7), rende del tutto trascurabile le probabilita relative di evapora- 
zione di tale particella. L’altro termine che compare in Q e cioè 29 — #,)- 
-(nB,—2B,) da nei vari casi: 


(per n) 29 — 0.) Bn = OD — 0) 
(per p) — 2(Ÿ — 20.) = — 2,4B(Ÿ — 3.) 
2 Ÿ — Ÿ, ME DES 
: ee (per «) 2(0 — B.)(2Bn— 2B») — 1,6$(® — À) 
(per d) 2(9 — 8.)\(Br— Bp) =— 9,86(0— 2%.) . 


Tale termine, cioè, favorirà l’emissione dei neutroni o dei protoni a seconda 
che (9 —9,) & positivo o negativo. Esso controlla l’emissione in modo che 
N—N, non diventi troppo grande. 

Esplicitamente nei tre casi considerati, la (4) del paragrafo precedente si 
scrive: 


— per i neutroni 


I,(e)de oc 7, ¢ exp [2Va(B— 8,6 + 1,68(9— 8.) — e)] de ; 


— per i protoni 


I,(e)de & Pre exp [2Va(B — 8,6 — 2,48(9— 0.) — e)]de ; 
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— per le « 


L(e)de o 27, ¢ exp [2Va(E — 6,2 — 1,66(9 — 9.) — e)]de. 


Il fattore 2 per le particelle « deriva dal fattore g che € comprensivo degli 
stati di spin e della massa del frammento evaporato. Per i nuclei intorno 
all’argento & stato preso, seguendo il WEISSKOPF, a = 6 e per quelli attorno 
al bromo a = 3. Questi valori di @ sono in accordo anche con i valori teorici 
dati da FERRARI e VILLI (2), i quali considerano a variabile con la temperatura. 
Se si pensa il nucleo come un gas ideale completamente degenere (HE == WO). 
si avrebbe per a un valore troppo elevato (11 MeV-! per A = 100 e Z = 44). 
FERRARI e VILLI, tenendo conto del fatto che l’effetto delle interazioni nucleari 
é quello di diminuire la degenerazione del gas, trovano che a diminuisce al 
crescere della temperatura. Risultano cosi giustificati i valori che abbiamo 
preso per a, trovati empiricamente da WeıssKopr. Tali valori vanno anche 
d’accordo con quelli che GUGELOT (1%) trova sperimentalmente. 


4Ë (cattura da coppia P-N) 


AE (cattura da coppia P-P) 
2 


12 16 20 24 28 


Fig. 7. — Spettro delle «. 


Nella fig. 5a) e 5b) sono riportati gli spettri di neutroni e protoni per i 
due tipi di cattura; nella fig. 7 & rappresentato lo spettro per le particelle «. 


Conelusioni. 


I risultati su cui si deve discutere perché si prestano ad un confronto speri- 
mentale sono i seguenti: 


1) Distribuzione delle stelle in funzioni dei rami ionizzanti. 
2) Spettro dei protoni di alta energia (E > 30 MeV). 
3) Rapporto «/p. 


4) Numero medio dei neutroni per cattura. 


(2) F. FERRARI e C. Viczr: Nuovo Cimento, 10, 94 (1953). 
@°) Di €. GUGELOT: Phys. Rev., 81, 51" (1951); 
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Pertanto riassumiamo dapprima i risultati sperimentali esistenti sui punti 
citati (+4). 

1) Nella tabella I sono riportate varie serie di risultati che si riferiscono 

alla cattura in emulsione fotografica, senza discriminazione fra nuclei leggeri 
e pesanti. 


TABELLA I. 


| Numero rami 
ionizzanti | | j | 


MENON e Coll. 
(1800 stelle) 
= | | = 
| MENON +ADELMAN 
(3366 stelle) 


| 33,9 +1,4 31,3+1,3 | 21,8-+1,1 | 10,9+0,8 | 1,9+0,3 | 0,1+0,1 
| 


| | 


| 32,441,0 | 32,7+1,0 | 22,3+0,8 | 10,5+0,6 | 2,0+0,2 | 0,1+0,1 | 


| | 


| PERKINS | | | 


26 | 3% 
(112 stelle) ae 26,8 19,6 9,8 0,9 
| = | 
ENOMOTO a | a) | 
T ave | 40,0-45,8 | 28,0783 | 21,51% CPAS | 1,0 268 | — 


(200 stelle) = | 


| | | | 


‚ Nostri risultati | | 
1230|263123 | Me MER se a 
| (500 stelle) 44,2+3,0 | 26,3+2,3 | 19,54+2,0 | 8,4+1,3 | 1,6-L0,6 | 


À facile vedere che, ove si tenga conto del differente modo con cui viene 
definito dai vari autori un ramo, vi € un sostanziale accordo. (MENON e coll. 
e ADELMAN considerano come ramo qualunque gruppo di grani formante una 
traccia definita, che non sia di elettrone. Noi, insieme a PERKINS, e ENOMOTO 
e coll. consideriamo come rami solo le traccie > 5 u. | 

Pit delicato & il problema della separazione tra nuclei leggeri e pesanti, 
perchè in base ai vari criteri selettivi (4) & possibile discriminare con relativa 
sicurezza solo una frazione delle catture. Questa separazione è riportata in 
tabella II e non verrà discussa in questa sede. 

Una parte dei mesoni x che terminano in emulsione non da luogo a stella 
visibile. Per questa percentuale è stato trovato da PERKINS il valore 29, 
MENON e coll. 28 + 2, ADELMAN 27 + 2, noi troviamo 27 + 3. MENON e coll. 
attribuendo questo tipo di cattura esclusivamente a nuclei pesanti o a idro- 


(4) D. H. PERKINS: Phil. Mag., 40, 601 (1949): M. G. K. Menon, H. MUIRHEAD 
e O. Rocuat: Phil. Mag., 41, 583 (1950); F. L. ADELMAN: Phys. Rev., 85, 249 (1952); 
8. Enomoro, Y. Furmoro, $. Horie e Y. Tsuzuxki: Progr. Theor. Phys., 7, 353 (1952). 
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geno, trovano che nelle catture in elementi pesanti, nel 36% dei casi non si 
ha emissione di particelle ionizzanti. Questo numero puod essere considerato 
in realtà un limite superiore, poichè non si pud escludere che gli elementi leg- 
geri non diano luogo a stelle a zero rami, pur ammettendo che questo avverra 
molto pit frequentemente in elementi pesanti. 


TABELLA II. 


| 
| Numero rami ionizzanti 1 


ag: 3 4 5 | Totale 
| | | 
3 B Fate ea Ile Ik 
| fl 1 13 | 15 6 iC te 
| COOL RTE 1 - i 5 | j 
Se eects 184 4 1 54 
ea ea ine eo ee |- | [ices ee 
| | | 
| leggeri . . . | 8789 2 13536 1716,39, A Bd 46,2 
a oeeet Tisulatt a pesanti', . >. 1435405130 1032 © 2,000 090ln 558 
= er = = a 
| Percentuale à: MENON e coll. 50,0 | 33,3 Hees | ape! 1,9 100 
| per i pesanti | Nostri risultati | 65,8 | 24,2 5,9 3 0,4 100 


2) MENON e coll. hanno trovato che 11,9 + 1,4 per cento delle stelle 
(escluse le stelle a zero rami) hanno associata una particella di carica |e| con 
energia >30 MeV; ADELMAN trova 9,5 +1,3. Facendo la separazione tra 
elementi leggeri e pesanti MENON e coll. trovano che 21,4 + 2,5 per cento 
delle stelle o in elementi pesanti hanno associata una particella di carica le 
con energia >30 MeV. 


3) Il rapporto «/p relativo alle sole stelle in elementi pesanti & stato 
trovato da PERKINS e da MENON e coll. 0,3 + 0,04. 


4) V. T. Coccont e EDWARDS (5) hanno trovato che il numero medio 
di neutroni emessi in seguito a cattura di x- a riposo varia proporzionalmente 
al numero atomico dell’elemento assorbente. Dai dati forniti si ricava che per 
i nuclei pesanti delle emulsioni il numero medio dei neutroni emessi per cat- 
tura à — 5,8. 


Si tratta ora di confrontare i risultati del Montecarlo con i dati sperimentali. 

La prima colonna della tabella III si é ottenuta considerando che la cat- 
tura del mesone avvenga da parte di una coppia N-P che si ripartisce energia 
ed impulso e considerando che la superficie del nucleo abbia un potere riflet- 
tente quale é da attendersi quantisticamente per incidenza normale della par- 
ticella. Il confronto coi dati sperimentali è molto soddisfacente. 

Nella seconda colonna sono riportati i risultati ottenuti considerando la 


(5) V. T. Cocconr e D. A. Epwarps: Phys. Rev., 88, 145 (1952). 
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cattura da parte di una coppia P-P, con la stessa ipotesi sul coefficiente di rifles- 
sione della superficie. Come si vede la situazione é peggiorata, perché aumenta 
fortemente lo spettro dei protoni di alta energia, e pure fortemente diminuisce 
il numero di cattura senza rami ionizzanti. 


TABELLA III. 


| Cattura | Cattura | Cattura 
| da N-P(1) da P-P da N-P (II) 
| | | 
| rami 53,9 52,6 34 
Distribuzione per- Butte oe 4 ne a 
n ; » 2 35:D 37,8 314 
centuale in rami, 3 9.9 8.9 | 24 5 
escludendo le stelle ; À rags + | hie 


a 0 rami in | 0,7 | 9,6 


~ 


— 
| 


| Percentuale di cattura cui non sono | | 
| associati rami ionizzanti 24,0 | 6,25 | 16,5 


Particelle di carica |e| con energia 
> 30 MeV. % relativa alle sole 


stelle con rami ionizzati 18,4 45,5 | 182 | 
| Fe 4,59 | 4,24 | 4,535 
| 3 Pa 0,945 | 1,28 1,05 
| Particelle ss r cattur ‘ | 7 
| i emesse per cattura | _ | 0.255 | 0.194 | 0 70 
| | | 
| ap | 0,27 | 0,15 | 0,67 


Pud, perd, sorgere il dubbio che il trattare il potere riflettente della super- 
ficie del nucleo come abbiamo fatto noi non sia corretto e che si debba portare 
alle sue estreme conseguenze il modello ottico, considerando che il coefficiente 
di riflessione sia funzione dell’angolo di incidenza. (Siamo grati al dott. KIND 
per molte discussioni sull’argomento). Abbiamo considerato anche questo caso 
estremo, ottenendo un paesaggio che scarta notevolmente dai dati sperimentali, 
come si pud vedere dalla terza colonna della tabella III. D’altra parte si puö 
ragionevolmente pensare che il modello ottico debba cadere in difetto quando 
la lunghezza d’onda dei nucleoni diventa confrontabile con le distanze inter- 
nucleari; e pertanto se si deve tener conto dell’angolo di incidenza sulla super- 
ficie nucleare, questo ha senso solo fino che l’energia dei nucleoni non supera 
i 10 MeV. Ma in tal caso i risultati non si scostano sostanzialmente da quelli 
del caso semplice trattato all’inizio. 

Resta ancora un dubbio importante che potrebbe alterare sensibilmente le 
conclusioni; si tratta del problema delle densita dei livelli energetici dei nuclei 
che evaporano, caratterizzati dalla grandezza a. Il valore da noi impiegato 
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è quello elaborato da WEISSKOPF (8), perd negli ultimi tempi. vari autori (16), 
hanno cercato di misurare il detto valore con esperienze dirette sullo spettro 
dei neutroni emessi nelle reazioni nucleari, pervenendo a dei valori notevol- 
mente pit elevati di quelli di WEISSKOPF. D’altra parte in questi lavori non 
viene separato il contributo dei neutroni di knock-on. Sorge allora il dubbio 
che, se si tenesse conto accuratamente dello svolgimento dei vari processi, il 
risultato potrebbe essere spiegato con un valore più piccolo di a: inoltre resta 
da dire che se il valore di a dato da questi autori rende conto dello spettro 
dei neutroni, questo stesso valore non pud render conto dei risultati di ana- 
loghe esperienze in cui viene esaminato lo spettro delle particelle cariche anzichè 
quello dei neutroni; infatti l’aumento del valore di a conduce a un impedi- 
mento quasi totale dell’emissione di particelle cariche, specialmente alle basse 
energie ove la teoria evaporativa & soddisfacente. 

In conclusione possiamo tentare di affermare che il risultato da noi pro- 
spettato per primo è quello che è più coerente col quadro attuale delle reazioni 
nucleari di moderata energia; pur non potendo escludere l’intervento di mo- 
deste frazioni di processi di altro tipo (cattura da P-P, cattura da «). Questo 
risultato trova riscontro nel fatto gia sottolineato da vari autori (17), che per 
quanto riguarda la cattura nucleare di mesoni z in volo, Vipotesi più plausi- 
bile & che essa si realizzi per assorbimento da parte di una coppia N-P. Si 
dovrebbe cosi concludere che non vi è nessuna ragione contraria alla affer- 
mazione che la cattura dei x di qualunque energia in seno alla materia nu- 
cleare avvenga quasi esclusivamente da parte di una coppia N-P. 

Se si accetta questa proposizione bisogna poi chiedersi quali siano i fattori 
che determinano questa situazione privilegiata per la coppia P-N. Se in realta 
il modello a gas di Fermi descrivesse il comportamento dei nucleoni nel nucleo 
dovremmo concludere che cid dipende dal fatto che elemento di matrice per 
assorbimento da parte di una coppia P-N e molto maggiore di quello relativo 
all’assorbimento da parte di una coppia P-P; di questo fatto si potrebbe avere 
la conferma confrontando, in base al principio dell’equilibrio elementare, questi 
due processi di assorbimento con i corrispondenti processi a ritroso di produ- 
zione e cioè 


= +tP+N>N+N N+N>P+N+7 
= +P+P-P N RP N-P+-P+mw. 


Comunque, in realtä, nei nuclei le cose sono molto più complicate, poichè 
non & possibile escludere in essi l’esistenza di ben definite correlazioni tra i 
vari nucleoni. Se, ad esempio, come & stato suggerito (17), in seno alla materia 


(1) E. R. GRAVES e L. Rosen: Phys. Rev., 89, 343 (1953). 


(7) K. A. BRUECKNER, R. SERBER e K. Watson: Phys. Rev., 81, 575 (1951). 
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nucleare vi fosse una forte tendenza ad associazioni in strutture parziali a 
deutoni, la ragione della posizione privilegiata della coppia P-N nell’assorbi- 
mento diventerebbe banale. 

Ma non è facile attribuire una tale situazione privilegiata alla sola coppia 
P-N e sembra più ragionevole ammettere che esista la tendenza alla forma- 
zione di piccoli complessi che vanno dal deutone al raggruppamento x. In 
ciaseuno di questi raggruppamenti, se l’assorbimento avviene sostanzialmente 
da parte di una coppia di nucleoni, il gruppo P-N & statisticamente privilegiato, 
a parte la questione, precedentemente discussa, che l’elemento di matrice sia 
più grande. 

L’unica obbiezione a questa generalizzazione potrebbe consistere nell’osser- 
vare che, se la struttura assorbente è un poco più complessa del deutone, 
Venergia totale del mesone si pud ripartire su più di due nucleoni. Si potrebbe 
tentare di dare una risposta teorica su questo punto, in base alla importanza 
relativa delle forze a due corpi rispetto a quelle a più corpi tra nucleoni, ma 
© forse più immediato segnalare che una esperienza diretta (18) di cattura di rn 
in He tende a mostrare che il risultato piü frequente dell’assorbimento consiste 
proprio nella distribuzione della energia del mesone sostanzialmente su due 
neutroni. 

Siamo cosi confidenti che i risultati da noi ottenuti, pur essendo relativi 
ad un modello nucleare molto schematico, descrivano largamente la genera- 
lita dei fatti, perchè comunque si pensi ad una struttura € ad un ordine locale 
nella materia nucleare, il risultato generale dell’atto di assorbimento & sempre 
la partenza di due neutroni veloci dal punto di cattura. 


E per noi gradito dovere ringraziare il prof. VALLE, direttore dell’Istituto 
di Fisica di Bologna, per la cordiale ospitalità e per aver in una prima fase, 
messo a disposizione i mezzi per questo lavoro, poi svolto con il contributo 
dell’I.N.F.N. 

Ringraziamo ancora vivamente tutti coloro che ei hanno aiutato nello svol- 
gimento del lavoro, e in modo particolare il dott. TOMASINI per il prezioso 
aiuto nella parte relativa alla evaporazione, le dott. FERRETTI e LANZARINI, 
per l’aiuto prestato nella estrazione degli esemplari, il signor TESTONI per la 
costruzione del modello meccanico. Ringraziamo ancora infine sentitamente 
il prof. BERNARDINI e il dott. GOLDHABER per averci dato le lastre nucleari 
necessarie alla ricerca. 


(18) C. P. Leavetr: Phys. Rev., 86, 603 (A) (1952) e tesi M.1.T. 
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SUMMARY (*) 


This paper closes a work about the absorption of mesons ~~ bound on the K-orbit 
by middle-weight nuclei (A=100). The problem is treated by the Montecarlo 
method and results are compared with experience. The distribution of stars against 


the number of visible prongs, the percentage of fast protons (H > 30 MeV), the emiss-: 


ion ratio «-particles/protons, the percentage of captures without visible prongs agree 
with experimental data if we suppose that the primary act of absorption in matter 
is performed by a proton-neutron pair. Remembering the results of the capture by 
nuclei of z-mesons in flight one might conclude that the capture of x-mesons in nuclear 
matter at all energies (at least up to 100 MeV) is essentially due to a proton-neutron pair. 


(*) Editor’s translation. 
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Mass Quantization and Isotopic Spin in Non Local-Field Theory. 


J. RAYSKI 


. . . . w 
Nicholas Oopernicus University - Torun, Poland 


(ricevuto il 31 Ottobre 1953) 


Summary. (*) — Starting from general remarks on the theory of fields, 
the local formulation of which seems to be insufficient, the author presents 
a non-local generalization of the Dirac equation that leads, by means of 
convenient supplementary conditions, to a mass spectrum for the particles 
of spin 4. The analogies of this mass spectrum with those obtained 
from the Yukawa and Pais theories are discussed. 


(*) Editor’s care. 


1. — A general remark. 


Physical quantities may be divided into two classes: universal and specific. 
Examples of universal quantities are: space-time variables, i.e. operators of 
position #,, displacement operators d,, spin variables, and isotopic spin va- 
riables. Examples of specific quantities are: wave functions of electrons, 
(nucleons, etc.), momentum of the electron (nucleon, ete.), field and so on. 
Universal quantities provide general means for description of physical reality 
and constitute a geometry in a broad sense. (Owing to the appearance of 
displacement operators it is a dynamical geometry). The specific quantities 
are formulated in terms of those geometric, universal quantities. 


9. — Defective character of the local field theory. 


In the traditional field theory the field quantities y are functions of the 
spin variables « and of the space-time variables @,. But x, are operators 
non-commutable with the displacement operators d,, and from a broad geo- 
metrie viewpoint the lack of a dependence of y upon d, should be regarded 
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as a defect of the formalism. The displacement operators are involved only 
in the field equations, e.g. in the Dirac equation 


(1) VAL p(æ), d,]- + xy(x) = 0, 


where y(x) = p(#, «) = y,(#) is a spinor field quantity. The local fields are 
only a special case characterized by four supplementary conditions 


aa 
bo 
= 


Ly, æ,]- = 0, Dh = Mopson ti 


3. — Non-local fields. 


Assuming generally a dependence of field quantities on universal quantities 
we have 
(3) y = ylx, d) where [x,,id,]- = 70, , 
(moreover y is a function of spin variables, isotopic spin variables, and perhaps 
other yet unknown universal quantities). In consequence of (3) the field quan- 
tities are operators in a twofold sense: as functions of the operators #,, d, 
and as quantities subjected to field quantization. In the «-representation 
the geometric operators, x,, d, as well as the field quantities become matrices 


(4) CNC Ney = 21,00 w'\d,| 2 = — — (a x"), 
(24 | | 
Dir hala, a) | a> 


satisfying the matrix algebra. Alternatively, the field quantities may be re- 
garded as functions of YUKAWA (1) variables 


(6) p= (0) MEDS 0, = ha, a), it 


There exist two independent linear forms in y and d,, namely [y,d,l-. In 
the #-representation we have 


[ ante, 

= apt | all ° HN sr ¥ 

(7) ep, dlz | a (2 + aa] ee ai] 9 ya, r) 
| or, 


() H. Yuxawa: Phys. Rev., 77, 219 (1950). 
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There is a striking resemblance between the transition from classical to quantum 
theory and that from local to non-local theory. In the local limit <@'|y|x") 
becomes a diagonal matrix 


(8) <a’ |yla”> > p(a')d(a’— 2") or via, r) > p(x) 


so that (7) tends to 
0p(x) 
Lu, dle|o"> >] Wu 
0. 


(9) <a! 


Thus, the anticommutator between % and d, vanishes in the limit of local 
fields while in the general (non-local) case it is different from zero. 


4, — A generalized Dirac equation. 


From a non-local viewpoint the equation (1) is obviously defective since 
it contains only one of the two independent linear forms in d,. A natural 
generalization of (1) will be 


(10) CARD d,]- ap BLY, d+ a U 0. 


Corresponding to a twofold number of degrees of freedom in the (+, d)-space 
the number of spin variables and of Dirac matrices is also doubled. We 
assume 


(11) [&,; œ,|+ = CB, B,}+ = 20, ’ [&,; By ly =0. 
Written in the x-representation the equation (10) is 


1 0 N. 0 er mal IN 
(12) Fr. a than + oa) ya) = 0, 


where Pas y, denote another set of Dirac matrices satisfying also the relations (11) 
and connected with the former by 


; il 1 
(13) ae v3 (C7 = Pa) ’ Vie V2 (CA tz By) : 
In terms of Yukawa variables the generalized Dirac equation is 


a ° 
(14) (a, a 2B. Fr, ate «) p(x, r) = 0. 
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The operator 


0) 
(15) M = %— 2,5, 
1 


obviously plays the role of a mass. Its eigenvalues are supposed to represent. 
the rest masses of spinor particles. In this way the non-local field theory leads. 
in a natural way to mass quantization. | 


5. — Necessity of supplementary conditions. 


Examples of mass operators of a type similar to our (15) have been in- 
vestigated recently by YUKAWA (?). However, the eigenvalue problem for 
such four-dimensional operators meets with serious difficulties (connected with 
the non-definite character of the metric form). Indeed, dealing with the eigen- 
value problem in exactly the same way as one used to do in the older quantum 
mechanics with three-dimensional operators in an euclidean space Yukawa 
found expressions contradicting the Lorentz invariance (*). The reason of 
this discrepancy is that eigenvalues and states may be defined only with 
respect to a space-like surface. Time-like components cannot be used for 
defining discrete eigenvalues. Thus, we do not believe that for such fourdi- 
mensional operators the eigenvalue problem may be formulated at all. The 
non-local formalism is yet too broad and has to be limited so as to re- 
strict 7, effectively to a space-like fourvector. Thus, we need some supple- 
mentary conditions [being a substitute for (2)]. 


6. — Reciprocity. 


Most probably the character of such supplementary conditions is universal 
(not specific). In other words: the supplementary conditions should be the 
same regardless the (spinor or tensor, etc.) character of the field quantities. 
Thus, Dirac matrices are not expected to be involved into the supplementary 
conditions. Moreover, universal conditions (restricting the effective range of 
the variables r,) being a matter of geometry are supposed to possess covariance 
properties characteristic (intrinsic) for the geometry itself. Now, the (x, d)-space: 


(*?) H. Yuxawa: Phys. Rev., 91, 415 (1953). 
(*) E.g. in the formula (2) of his second Letter to the Editor appears an exponential 
factor with r? + 7? while in a relativistic theory it should be rather |r? — »2|. 
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is invariant against the reciprocal transformation (*) 
1 =1 N ; nr, 
(16) 0, —iM,, id, >— 1%, 
where À is a constant introduced for dimensional reasons (with € — # — 1 


the constant À is of dimension of a length). It seems plausible to require the 
supplementary conditions to be invariant against (16). Such conditions 


lia) [Cw, d,]-; CAE = [Ly, CAES d,)- U ? 
(17) 2[[y, d,\-; d,,\- = A I[y, æ,]-; æ,]- ? 


have been introduced earlier by the author (5%). With (17’) a universal length 
À is introduced into the formalism. With 


u 


(18) ar = = (Ae, F tAd,) , 


the supplementary conditions may be written also in an equivalent but more 
elegant form 


(19) pain ai]- = [[y, a,b, @,]-=9- 


in terms of Yukawa variables the supplementary conditions are 


fs ‘ 22 
(20) r =. pla, r) =0, (22 =i = A tri) (a, r) = 0. 
u 


m 


Introducing Fourier transforms y(p, r) the Dirac equation (14) and the con- 
ditions (20) become 


3 © 
(21) (iso — Dur FR «) w(p,r)= 0, 
u 


(22) DT wt) =O, Are) = 0 - 


(4) M. Born: Rev. Mod. Phys., 21, 463 (1949). - 

6) J. Rayskı: Proc. Phys. Soe., 64, 957 (1951); Acta Phys. Pol., 11, 109 (1951). 

(6) Instead of (17/) YUKAWA assumed a little different second condition which 
is not self-reciprocal : 


(17Y) Ip, &,1-, %J- =A’ - 
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For a particle at rest p=0, p—M we get 


(22') | x M— 28, > er “)y=0, 
or, 
(22") KM p= 0; (— AM? + Ary = 0. 


From (22') it is seen that (with p=0) y(p,r) must vanish unless 7,=0, 
r — 42|M| so that the range of variables of r, is effectively restricted to polar 
angles 9, 0 on a sphere with radiun r=/?|M|. With the alternative auxi- 
liary condition (17Y) the radius is simply A. 


7. — Mass spectrum. 


To determine the eigenvalues of M we have to solve the equation (21') 
with the constraints r,—0, r=/?|M|. To this end we go over to polar 
variables (7) 

(23) B; 2 = ps ie, AH je. Oy eviter 

er, or : ê 

where t denotes Pauli matrices in the space of the / indices, U is the orbital 
angular momentum in the space of r, eis a matrix commuting with + and 
anticommuting with B, with the property «?=1. e together with ß,, and + 
may be represented as sets of independent matrices with two rows and columns 
whose eigenvalues are +? for Be and + (j+) for vl! +1. The terms 
with 0/0r, and 0/0r disappear in consequence of the constraints and we are 
left with 


(24) (+ 


This yields two families of particles (and antiparticles) with masses M and 
MS satisfying the equations 


1 
M + 2 ( E= à ip rn «) Vir=almh = 0 with I = 


NI © 


1 
5? 


j + + W 
‘ Ge Ge, Ne 
(25) M; Een Du : MS x — Sr 2 


Each family splits into two subfamilies with absolute values of masses 
— 
= 7 Be. (2) | 2(; à alt 
| 2 
AMY) = a ee 2(j +)|; 
2)! 


() P. A. M. Dirac: Quantum Mechanics (Oxford, 1947), third ed., p. 266. 


(26) 


MASS QUANTIZATION AND ISOTOPIC SPIN IN NON-LOCAL FIELD THEORY 1735 
where + refers to i=1,2. The number of masses of the first family is 


infinite, while that of the second family is finite and restricted by the con- 
ditions 27 < (Ax/2}° — 1. 


8. — The formalism of Pais (°*). 


By using the auxiliary condition (17Y) instead of (17") we have to replace 
Æ|M| by A in (25) which yields a simpler, linear mass spectrum 


(26) AM =| Ax + 2(j + %) 


» EN 


which is identical with the mass spectrum in the formalism of Pais if 
our x 2/A is identified with the mass of the nucleons. (Pais’s A is smaller 
from Yukawa’s by a factor 14). Indeed, by interpreting the matrices 7 ap- 
pearing in (23) as isotopic spin matrices, the non-local formalism with the 
auxiliary condition (17Y) is completely equivalent with Pais’ formalism. It 
seems to us, however, that the condition (17”) leading to the more elaborate 
mass spectrum (27) is preferable. 
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(8) A. Pais: Proc. Lorentz and Kamerlingh Onnes Conf. Leiden (1953), June 22. 


PACA ISS) LUPINE 


Partendo da osservazioni generali sulla teoria dei campi, la cui formulazione locale 
appare insufficiente, viene proposta una generalizzazione di tipo non locale dell’equa- 
zione di Dirac, la quale conduce, con l'uso di opportune condizioni supplementari, a 
uno spettro di masse per le particelle a spin 14. Sono discusse le analogie di tale spettro 
di masse con quelli che si ottengono nelle teorie di Yukawa e di Pais. 


‘(*) A cura della Redazione. 
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(ricevuto il 3 Novembre 1953) 


Summary. — Further evidence is presented in favour of the existence 
of charged unstable particles of hyperprotonie mass. In two cases the 
decay product is a fast light partiele, in one case it is a slow proton. The 
decay schemes: QŸ +n + 7° + Q: QO” > p+r0+Q are proposed 
for the two types of events respectively. They are possibly alternative 
modes of decay of the same type of H-particle. The evidence available 
from photographic emulsion and cloud chamber experiments is discussed. 


The observation of the decay at rest of a particle (J,,) heavier than the 
proton has been reported in a preceding letter (1). Two other events which 
can be also interpreted as the decay of H-particles (**) are described here: of 
these, the first (J,,) is similar to J, , Since the secondary is a fast charged 
particle certainly lighter than a proton; the second (J,,) on the other hand shows 
a slow proton arising from the end of the primary track . 


(*) Communicated at the Cosmie Ray Conference of Bagnères de Bigorre, July 1953. 
(**) Nomenclature proposed at the Cosmic Ray Conference, Bagnères de Bigorre, 
July 1953. 


() A. Bonwrrr, R. Levi Ser, M. Panermı and G. Tommasini: Nuovo Cimento, 
10, 345 (1953). 
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Since these examples of H-particles have already been indicated with 
the symbol J in preceding communications, we find it convenient, to avoid 
confusion, to retain in this paper the old denomination. The subfixes / and 
p indicate that the charged secondary is a light meson or a proton. 


1. — Mass measurements on J-partieles. 


The mass of J-particles has been determined from the conventional «-R 
and G-R measurements. The proton tracks employed as standard for relative 
mass measurements by gap- 


counting were taken from the, 
vicinity of the J-tracks, at 
the same depth and with the 
same dip. oe \ 
. The mass of our K-me- 

sons and heavier particles has ca! eR Bes ace eye 
been also determined by C. t K 
DILWORTH, making use of 
the « constant sagitta me- 
thod» (24) (%). A very good 
agreement has been found 


between the new results and Fig. 1. — Mean angle of scattering v. range for 
particle J,. The first two points refer to inde- 
pendent sections of the track of about 4mm in 
the first plate; the third one corresponds to à 
value is taken of scattering section of the same length in the second plate. 
and gap determinations, the See (1). 


0 


0 2 4 6 Ü 10 12 4 16 Ru) 


those previously obtained. In 
ths paper, when a mean 


«-R mass value employed is 

that given by the constant sagitta method: this choice has been recom- 
mended in the standardisation meeting at the Varenna Summer School (Sum- 
mer 1953), 


a) Events J, and = Event J, has been already described in à pre- 
liminary note (1), where the results of the mass determination by «-R (angular) 
method was given. A test of the inner consistency of these measurements is 
shown in Fig. 1: the three points correspond to %o, determined on inde- 


(2) C. DILWORTH, S. GOLDSACK and L. Hirscupera: Cosmic Ray Conference, Ba- 
gnères de Bigorre, July 1953, and Nuovo Cimento, in press. 

(2) S. Biswas, E. C. GEORGE and B. Perers: Cosmic Ray Congress, Bagnères de 
Bigorre, July 1953. 

(4) T. HOLTEBEEK, N. ISACHSEN and 8. 0. SORENSEN: Phil. Mag., 44, 1037 (1953). 

(*) This method gives a better precision than the usual ones, particularly on short 
tracks: a calibration on proton and meson tracks gives a mean value of the proton 
mass of 1790 + 40 m, (?). 
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pendent sections of the track of about 4 mm. The plot strongly suggests 
that the direction of motion of the particle is towards the decay point at 
R= 0. 

Measurements by the G-R method consistently indicate a mass higher 
than that of the proton; however irregularities in the plot of the gap v. 
range, probably due to uneven development, made uncertain the attribution 
of a definite value of the mass. 

Due to shortness of the primary of event J,, the indications on its mass 
are less accurate; the direction of motion however was well established, showing 
that also this event represents a decay process. 


TABLE I. 
—— — —— 
| Primary y Dore 
| ee EEE =u eye rowers 
| Residual - | Mass, m, Mass, m, Mass, m, 

| Range, um | Constant Sagitta _G- R a-R 

PE a Zu ri Bu en a ee a a 
| | | 
SP: | 15 760 2210 +39 | >p 2 500 +8 
| | | 
| iy. 1250 2.340 2680 4 2 550 +530 2250 +300 


The information on both secondaries is meagre. They are too short (120 u. 
and 245 u respectively) for scattering measurements to be made; however 
no measurable deviation is observed along their path. Their grain density 
is near the minimum. It can be excluded that they are slow electrons arising 
from the radioactive decay of known particles, and may be either fast electrons 
or fast mesons, obviously not protons, since the energy available is insuffi- 
cient. The data concerning the primaries of the two events are collected in 
Table I. 


b) Event J,,.— Event J, is shown in the photomicrograph of Fig. 2. The part- 
icle J, comes to in end of its range at point 0 and there gives rise to a highly 
ionizing secondary which stops in the emulsion (p). The primary, travelling 
at an average angle of dip of 25° to the plane of the unprocessed plate, has a re- 
sidual range of 900 u, the secondary, travelling horizontally, has a range of 1670 + 
+18 u. The results of mass measurements on both primary and secondary are 
contained in Table II. The absence of a visible decay product at the end 
of the secondary particle, together with the results on mass determination 
(see Table II), strongly suggest that this particle is a proton. Under this 
assumption, we deduce from the range an energy of 18.7 + 0.18 MeV at 
the point of ejection 0. 
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Fig. 2. — Partial photomicrograph of event J,, Observed by R. Levi SETTI 


ON THE EXISTENCE OF UNSTABLE CHARGED PARTICLES ETC. 1739 


TABLE II. 


Primary Particle. 


Residual | Mass, m, Mass, m, | Mass, m, 
page | eR |e eee 
| | 
ip 900 233078 | 270078 | 1840 +88 


Secondary Particle 


Mass, m, 


Mass, m , 2 
Range, PE Enegy, MeV 
NE Constant Sagitta | G-R | 2 
ce x er SSS a 
Jo, 1670 +18 2030 1480 | 18404670 | 18,7 + 0,18 
| (Identified proton) | 


The assumption that event J, involves two distinct particles, one of 
which is considered to be the primary, the other the secondary, is supported 
by the results of the measurements which show that the event can hardly 
be interpreted as a scatter- 
ing or two pronged star. 

The graph of Fig. 3 
shows the variation of scat- 
-tering and of gap-density 
along the two tracks. The «0 = 
sense of variation of both \ 
measured quantities, before 
and after the well marked eo 
discontinuity observed at “à 
point 0, indicates that the 
two tracks are due to di- 


stinet particles: one of them, Fig. 3. — Gap and scattering measurements on the 
primary and secondary particles of event J,. 
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the primary, moves, slowing 
down, towards the point 0; 
the other one, the secondary, is ejected from the same point and then is 
stopped. The ionization of the primary near point 0 suggests also that the 
particle is either at the end of its range or very near to it. 

The possibility of the event being due to a chance juxtaposition of the 
end of a stopped particle with the origin of a knock-on proton has been con- 
sidered. Actually, from the density of these two types of events in our plates, 
one would expect one chance juxtaposition in about 107 cm? of emulsion. 
The probability of observing such a coincidence is therefore very small. 
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2. — Interpretation of the results. 


The two events J, and J,, belong to a group of particles heavier than the 
K-mesons, which decay into a single charged secondary. They are distin- 
euished from K-mesons on the basis of mass measurements on the primaries. 
: Six examples of K-mesons 
decay have been found up to 
now in our laboratories; their 
mass distribution is shown in 


os Nun e 


Gr Be 


[Gs Fig. 4 together with the va- 
cod nn "no" 10087200: me (me) lues obtained for J-events. The 
re four points indicated as « B» 
el; 8 correspond to the results on 


the K-mesons found by the 
Brussels Group (5) and kindly 
communicated to us. In these 
histograms the range of each 
particle and the type of measurements performed on it are indicated. 

The question arises whether or not the points corresponding to J-particles 
can be regarded as representing the tail of the distribution of the K-mesons. 
Let us assume that J, and J; belong to populations of particles of mass 
about 1000 m,, with dispersions equal to those of the two tracks respectively: 
one finds that the probability of getting the constant sagitta corresponding 
to the observed mass of J, is about 107; similarly the compound proba- 
bility of getting the observed values of constant sagitta and gap-density for 
J, is about 10-*. | 

Thus J :, and J;, should represent the decay of particles with a mass signi- 

ficantly different from that of K-mesons; we suggest that the two events are 
really due to hyperprotons of mass around 2300 m, . 

J,, requires a more complex analysis. We will discuss three possible inter- 
pretations of this event: 


Fig. 4. — Mass distribution of K- and J-particles. 


1) nuclear capture of a negative K-meson at rest or in flight; 
2) nuclear interaction of a deuteron, giving rise to the evaporation of 
a proton from the target nucleus; 


3) spontaneous disintegration of an unstable particle. 


If one considers the primary particle asa K-meson not at the end of its 


(5) C. Ditwortn, S. GOLDSACK and L. HIRSCHBERG: Cosmic Ray Conference, Ba- 
gnéres de Bigorre, July 1953 and Nuovo Cimento, in press. 
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range, one should expect: a) a mass determined from gap measurements lighter 
than that corresponding to a K; b) anon constant sagitta, the set of variable cells 
being that corresponding to a particle at the end of its range. A check has been 
made of the variation of the sagitta with residual range and, in the limits of 
the accuracy allowed by the short track, no significant departure from a 
constant level has been observed. On the other hand, the higher mass value 
determined by gap, already inconsistent with the mass of a K-meson at rest 
is still less compatible with a K-meson in flight. Let us then assume that the 
particle is a K-meson at rest: the probability of observing the high mass 
value obtained from both scattering and gap measurements is of the order 
ONE 

The second interpretation requires that the «primary » deuteron should 
arrive at the point 0 with a certain residual energy (*). Since the mass to be 
supposed in this case is higher than that measured, we find that the gap- 
density could be made to fit quite well that of a deuteron in flight. Let us 
assume a residual energy of about 10 MeV (energy comparable with the height 
of the Coulomb barrier for the heaviest nuclei of the emulsion): from scattering 
one should have found a non constant sagitta with an average value corre- 
sponding to the sagitta of a particle of mass around 10000 m,. The proba- 
bility that the observed value of sagitta is a fluctuation from the expected 
one is about 107+. 

It seems therefore highly probable that this event represents a decay pro- 
cess; the parent particle must then be more massive than the proton. 

Other examples of hyperprotons have been found both in photographic 
emulsion (*7) and in Wilson chamber experiments (*"’). Among the events 
observed by these authors there appear two types of decay. In some 
cases the light secondary was identified as a r-meson (5°). .Two decay 
schemes have been proposed by analogy with the decay of V} - particles. 
Adopting the nomenclature proposed at the Bagnères Conference these are: 


(a) OF >n+r"+Q, 
(b) Ve ae 


(*) See P. Fowrer’s intervention: Proceedings of the Cosmic Ray Conference, Ba- 
enéres de Bigorre, page 183, July 1953. 

(6) M. Crccarerıı and M. Merri: Cosmic Ray Conference, Bagnères de Bigorre, 
July 1953; Nwovo Cimento, 10, 1207 (1953). 

(7) D. Lar, Yasn PAL and B. PETERS: Cosmic Ray Conference, Bagnères de Bigorre, 
July 1953. 

(8) H. S. Brip@e and M. ANNIS: Phys. Rev., 82, 445 (1951). 

() B. Rossı and H. BrinGe: Cosmic Ray Conference, Bagnères de Bigorre, July 1953. 

.(%) A. M. York, R. B. LEIGHTON and E. K. BJORNERUD: Phys. Rev., 90, 167 
(1953) and Cosmic Ray Conference, Bagnères de Bigorre, July 1953. 

(1) W. B. FRETTER, B. G. GREGORY, R. Jounston, A. LAGARRIGUE, H. MEYER, 
PF. MüLter and C. Pryrou: Cosmic Ray Conference, Bagnères de Bigorre, July 1953. 
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The events of CECCARELLI and MERLIN (%) and of LAL et al. (7) are com- 
patible with the first decay scheme and give Q-values of 131 and 135 MeV 
respectively. BRIDGE and ANNIS (®°) found an event which could be inter- 
preted as an Q* with Q of 120 MeV. 

Our events J Ds and J, can be attributed to the class of Q= - particles. 


ri 


Their measured mass is in agreement with that calculated from the decay sche- 
me (a), assuming a Q-value of the order of that given by CECCARELLI and 
MERLIN and by Lat et al. The ionization of the secondary particles is consi- 
stent with that of a m-meson of about 100 MeV energy. 

The observed coincidence between the energies of the z-mesons ejected 
both in the disintegration of the Q+-particles and that of y-mesons (7) ma- 
kes it possible the separation between individual examples of Q=’s and 
K-mesons only on the basis of accurate mass measurements on the primaries . 

The fact that no example of slow K-particles decaying in flight in the 
emulsion has been reported, whereas among slow Qs there are two 
cases (7) of decay in flight, could be interpreted as an indication of a shor- 
ter mean life for these particles than for K-mesons. 


Tage III. — Charged Hyperons. 


| if | N | | Mass 
pe ass, m | 
JI x e | ©, MeV | from assumed 
of detector direct measur. | | 
| | decay scheme | 


| | 
+ 


| 
| Dy KO near | 
| | 

| 


Ja (ES) | Nuclear plate | 2210 +250 | ? | ? 

gs ; | 2400 + 580 ? | ? 

| Padua event (5) | » | 2100 + 400 131 + 24 | ~ 2370 
Bombay event (7) | 5 | 2420 + 330 - | 135 +35 | 2380-430 | 


M.I.T. event (%°) | Cloud chamber | — 120 ~ 2350 


| | 
b) OF >p+r+Q 


Te Nuclear plate 2500-5620 | 11543| 232043 
| | | 
Cal. Tech. | Cloud chamber = fe? | ~2180 
events (1°) | E40 | 
Paris event (1) | » | = | 26-38 | ~2170 


Among the events which can be interpreted as corresponding to the second 
decay scheme, two described by York et al. (1°) give a Q of about 40 MeV 


(2) M. G. K. Menon and C. O’CEALLAIGH: London Discussion on Heavy Mesons, 
January 1953. 
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A value of this order is reported by FRETTER et al. (1). If we attribute to 
event J, the decay scheme (b), we find a Q-value of 115 +3 MeV and à 
mass for the Q* of 2320 +3 m,. 

The data available on charged hyperons are collected in Table III. 

The Q-values obtained on the cloud chamber events of type Q* seem to 
be very close to that of the V}-particles, i.e. 37 MeV. On the contrary, the 
Q of J, is significantly higher and in agreement with that observed for events 
following decay a). Thus while the Cal. Tech. and Paris Group events could 
be regarded as the charged counterpart of V{-particles, J, as well as the Q= 
examples cannot be interpreted in this way. The agreement of the Q-values 
between J, and Q= suggests that the two types of events represent alternative 


TT 


modes of decay of the same particle. 
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RIASSUNTO 


Si presentano nuovi risultati in favore dell’esistenza di particelle cariche instabili 
di massa iperprotonica. In due casi il prodotto di disintegrazione è una particella leg- 
gera veloce, in un altro caso & un protone lento. Si propongono gli schemi 


OF sn +ri +0; OQ, +p +7 + 


per interpretare i due tipi di eventi rispettivamente. Essi rappresentano probabilmente 
processi di disintegrazione alternativi dello stesso tipo di particella H. Si discutono i risul- 
tati fino ad oggi disponibili riguardanti esperienze in camera di Wilson e in lastre 
nucleari. 


LETTERE ALLA REDAZIONE 


(La responsabilità scientifica degli scritti inseriti in questa rubrica à completamente lasciata 


dalla Direzione del periodico ai singoli autori) 


Seiami di almeno due elettroni prodotti in piombo da raggi cosmici. 


C. MILONE 


Istituto di Fisica dell Universita - 


Catania 


(ricevuto il 31 Ottobre 1953) 


Durante una ricerca sugli elettroni e 
fotoni negli sciami estesi sono state ese- 
guite delle misure preliminari sugli sciami 
di almeno due elettroni emergenti da uno 
strato di piombo quando su questo giunge 
un elettrone od. un fotone di origine co- 


4 
E(s) Eis) 
vat 
rt 3 
AR 
0.2 t/ $= E(s) 
$ = Fis) 
on |* E(h) 
a F(h) 
Oe L eee 
0 10 Ones: 
mm Pb 
(a) (b) 
Fig. 1 


smica (associato 0 no a sciami estesi). 
I risultati ottenuti possono essere con- 
frontati con i risultati di misure di altri 
Autori (1-4). 


() V. H. REGENER: Rice. Scient., 11; 66 
(1940). 
(*) N. NERESON: Phys. Rev., 61, 111 (1942). 


() J. Cuay e C. Levert; Physica, 12, 321 
(1946). 
(4) S. FRANcHETII: Nuovo Cimento, 10, 


551 (1953). 


Il dispositivo sperimentale adoperato 
è rappresentato nella fig. la; sono stati 
registrati contemporaneamente gli eventi 
A+B+Ce B+C al variare dello spes- 
sore s dello strato di piombo interposto 
fra A e.B. I tre contatori del gruppo 4 
erano riuniti in parallelo; un evento C. 
invece era costituito da una qualsiasi 
coincidenza del tipo 142 ovvero 2+3 
ovvero 1+3. Tutti i contatori avevano 
superficie utile di 20 x 20 em?. I risultati 
ottenuti sono riassunti nella tabella I. 

Coi valori dei rapporti 


ABO 
ae nn 
mın 


Bee rei) 
— = FG), 
min 

sono stati costruiti i grafici riportati nella 
fig. 1b. I valori indicati con H(h) ed 
F(h) sono relativi alle frequenze degli 
eventi registrati allorché sopra il grup- 
po A veniva posto uno strato di piombo 
di spessore h=150 mm. Questi eventi, 
attribuibili a particelle penetranti, costi- 
tuiscono una parte del fondo. Sottratto 
il fondo, le anticoineidenze F(s) possono 
essere attribuite per la maggior parte a 
fotoni producenti in s almeno due elet- 
troni; indichiamo tali eventi con f (n > 2). 
Se le coincidenze H(s), sottratto il fondo, 


SCIAMI DI 


ALMENO DUE ELETTRONI PRODOTTI IN PIOMBO ECC. 
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TABELLA 1. 


i | eel 
8 | | 
(Dee Seay min |A +B+0| B+C 
] e 
0 2 820 615 | 1089 
6 1 360 426 | 834 
| 9 1 370 466 | 951 
12 1345 445 957 
15 1405 430 958 
24 1500 | 424 948 
=== = ——— 
hmm Pb | : | : 
| sopra A | min |[4+B +0) B+C 
= a=} > 3 .- 
150 4065 845 824 


possono essere attribuite per la maggior 
parte a produzione di elettroni in s da 
parte di un elettrone, e indichiamo questi 
eventi con e (n > 2), riferendoci ai valori 
corrispondenti ai massimi delle due curve 


— er 
\(B+0)- -(A+B+C) | 


| A+B+G 
— = Eis) 


—— : ANS) 
min min 

0,218 + 0,009 0,618 + 0,008 
0,314 + 0,015 0,300:+ 0,015 | 
0,341 + 0,016 0,354 + 0,015 
0,332 + 0,016 0,380 + 0,017 

| 0,314 + 0,015 0,376 + 0,016 

| 0,283 — 0,014 0,348 + 0,015 

E(h) F(h) 


0,12 + 0,005 0,083 + 0,004 


di transizione si ricava 
n>2 
2 1,84 + 0,4. 
e(n >2) 


Questo valore non differisce eccessiva- 


Guo aN 


7 TS 
fin >2) | | Metodo | | 
| i t Ë 5 | D e | 
ein >2) | Autor | di misura Note À | 
= E ee a | 
| | 
| æ 1 | REGENER, 1940 | Contatori | In corrispondenza del mas- 
| | | | «imo delle curve di tran- 
| | sizione. 
| | | 
à NERESON, 1942 Contatori | In corrispondenza del mas- | 
| | simo delle curve di tran- | 
sizione | 
0/ ; | | 
Re | Cray e LEVERT, Contatori | Per s= 8 mm Pb | 
202722025 1946 | | 
2 +0,4 | FRANcHETTL 1953 | C.W. | Per s=20 mm Pb | 
| | | 
| | , | 
1,8 + 0,4 | Presente lavoro | Contatori | In corrispondenza del mas- | 


simo delle curve di tran- 
sizione. 
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mente da quelli ottenuti da altri Autori 
(tabella II) quando si tenga conto che, 
come hanno fatto notare DAUDIN e 
Tsar Cuv (+), le misure eseguite con soli 
contatori non possono fornire su questo 
argomento indicazioni esatte. Lo scarto 
0,4 aggiunto al valore 1,8 si riferisce 
quindi solo alle oscillazioni statistiche: 
eventi 4+(B+C) dovuti a sciami pro- 
dotti da un fotone associati con una par- 
ticella ionizzante che colpisce A, od 
eventi dovuti a più particelle indipen- 
denti possono avere influito sensibilmente 
sui risultati. Si nota perd (tabella IT) 
che le diverse misure eseguite con con- 


C. MILONE 


tatori forniscono in generale 


jin> 2) 
a LF 
> 2) 


e(n 


in accordo coi risultati ottenuti in ¢.W. 
da FRANCHETTI (4); cid mostrerebbe, in 
accordo con le conclusioni del predetto 
autore (2) l’esistenza di un eccesso di fo- 
toni rispetto agli elettroni per energie 
relativamente basse. 


() J. DAUDIN e TSAI-CHU: 
234, 1277 (1952). 


Compt, Rend., 


5 
Ve 
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Sull’equazione del moto dell’elettrone nell’elettrodinamica classica. 


P. CALDIROLA 


Istituto di Scienze Fisiche dell’ Universita - Milano 
Istituto Nazionale di Fisica Nucleare -.Sezione di Milano 


(ricevuto il 9 Novembre 1953) 


Sono note le difficoltà che si incontrano, gid nell’elettrodinamica classica (cioè 
non quantistica), quando si voglia determinare il moto di un elettrone sotto l’azione 
di forze esterne e di quelle derivanti dalla reazione dovuta alla radiazione emessa 
dall’elettrone stesso. 

Nella approssimazione non relativistica una equazione per tale moto è stata 
ricavata a suo tempo da Lorentz (!) nell’ipotesi che l’elettrone sia schematizza- 
bile come una piccola sfera carica indeformabile. Tale equazione & la seguente: 


(1) Mo = F+R, 


essendo la reazione R espressa da una serie di potenze erescenti del raggio À del- 
l’elettrone : 
2e? d’v 


= — 4 R + RE 
R ra )\R+( ) 1 


Si osservi esplicitamente che i coefficienti delle varie potenze di R contengono 
le derivate successive della velocita v e che inoltre la loro effettiva espressione cambia 
se si assumono modelli d’elettrone diversi da quello adottato da Lorentz di sfera 
rigida; in ogni caso pero i termini di R successivi al primo si annullano tutti quando 
si supponga, con un passaggio al limite, l’elettrone puntiforme (R — 9). 

Successivamente Drrac (?) ha sviluppato una teoria relativistica per l’elet- 
trone puntiforme arrivando all’equazione: 


du, 


2 MC 
(2) mo dé 


@) H. A. Lorentz: The theory of Electrons (Leipzig, 1916), 2* ediz. 
(2) P. A. M. Drrac: Proc. Roy. Soc., A 167, 148 (1938); Ann. Inst. Poincaré, 9, 13 (1938). 


sn =. Ss 
; 5 = 
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F e i 2 2e? [dus | dup ( 1) 
Dr == ES 5 = Uw UU, = — 
FE CAN ra 3e \ ds? | * P de? ene 


sono rispettivamente i quadrivettori della forza esterna e della reazione. La (2), 
nell’approssimazione non relativistica, coincide con l’equazione di Lorentz per il 
caso limite di un elettrone puntiforme, cioè alla: 


dv 2e? d?v 


La difficoltà maggiore che si incontra nell’adozione di equazioni del tipo 
della (1) o della (2) è legata al fatto che, comparendo in esse derivate temporali della 
velocità di ordine superiore al primo, i loro integrali generali contengono un numero 
di costanti arbitrarie superiore a due e, attribuendo a queste costanti valori generici, 
si hanno soluzioni deseriventi moti in netto contrasto con le previsioni dell’ordi- 
naria intuizione fisica. Per ovviare a questo inconveniente, diversi autori (*) hanno 
suggerito opportuni criteri tendenti a discriminare, tra le soluzioni analitiche delle 
equazioni, le soluzioni « fisiche » da quelle «non fisiche ». Tali criteri non hanno perd 
un carattere completamente generale. Per di più in alcuni problemi particolari, 
come quello, sia pure un po’ patologico, di un elettrone lanciato lungo una retta 
contro un protone « fisso », si hanno (4) come soluzioni delle (2) risultati ehe, indi- 
pendentemente dai valori assunti dalle costanti arbitrarie, non sembrano in aleun 
modo fisicamente accettabili. 

Le considerazioni precedenti ei hanno fatto ritenere che, se si vuole arrivare 
ad un’equazione soddisfacente per il moto dell’elettrone, & necessario riesaminare 
le ipotesi essenziali stesse su cui si fondano le equazioni (1), (2) e (8). 

Recentemente in diversi campi della fisica teorica (5) si & riscontrata la neces- 
sità di introdurre una lunghezza minima fondamentale (e connesso a questa anche 
un intervallo di tempo minimo fondamentale). Abbiamo voluto esaminare se, met- 
tendoci da questo punto di vista, fosse possibile ricavare qualche risultato interes- 
sante nello studio del moto dell’elettrone. In relazione a cid ci siamo proposti di 
stabilire un’equazione per il moto dell’elettrone prendendo come punto di partenza 
la seguente serie di postulati: 


1) Esiste una costante reale 7, dalle dimensioni di un tempo (intervallo di 
tempo fondamentale). 


(®) M. ScOHNBERG: Phys. Rev., 67, 211 (1945); C. I. ELIEZER e A. W. MAILVAGANAM: Proc. 
Cam. Phil. Soc., 41, 148 (1945); P. CALDIROLA: Nuovo Cimento, 5, 99 (1948); A. LoINGER: Nuovo 
Cimento, 6, 360 (1949). 

G@) ©. I. ELIEZER: Proc. Camb. Phil. Soc., 33, 173 (1943); G. Zin: Nuovo Cimento, 6, 1 (1949). 

(5) A. Marcon: Quantum mechanics of particles and wave fields (New York, 1951), p. 267 e segg.; 
A. LAND: Journ. Franklin Inst., 229, 767 (1940)); B. T. DARLING e P. R. ZILSEL: Phys. Rev., 91, 
1252 (1953). 

Della nozione di lunghezza elementare, e quindi di un intervallo di tempo elementare, si fa 
pure uso in tutte le odierne teorie dei campi non locali. 

Il punto di vista da noi adottato, che © sostanzialmente quello di DARLING e ZILSEL, non com- 
porta aleuna quantizzazione, nè alcuna concezione di struttura reticolare del continuo spazio-tem po 
einsteiniano, 
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2) Una forza, operante su di una particella all’istante 7 del tempo proprio di 
questa, determina un brusco passaggio dallo stato di moto relativo all’istante 7 —7+, 
allo stato di moto relativo all’istante 7. In altri termini le leggi della meccanica 
connettono grandezze dinamiche della particella relative a punti della linea oraria 
cronotopica distanti l’uno dall’altro di una lunghezza s, = er,. 


3) Per 7, 9 le nuove equazioni si devono ridurre a quelle della dinamica 
classica. 


4) Le leggi devono essere invarianti relativisticamente. 


Per semplieitä determiniamo dapprima le leggi del moto dell’elettrone quando 
si lasei cadere la condizione 4). Detta allora Fr, t) la forza che agisce sulla parti- 
cella all’istante f, essa produrrà una variazione sul moto di questa in modo da farle 
acquistare una quantita di moto myv(t). Tenendo conto dei postulati 1), 2), 3) assu- 
meremo allora a base della nostra teoria le seguenti equazioni alle differenze finite: 


(4) ol ae rt). 


To 


Sviluppiamo v(t — r,) in serie di potenze crescenti del parametro 7: 


; dw , Ta dv 
v(t — To) = v(t) —T ro at 


In approssimazione zero avremo 


che evidentemente & l’equazione del moto di una particella sottoposta alla forza 
esterna F e senza reazione. 
In prima approssimazione avremo 


dv(t) 7, d2v(t) 
dt 2 dt? 


= F(r;t), 


Mo 


€ 


che si riduce all’equazione di Dirac non relativistica (3), quando si assuma 


4 e? 
Ton = ate ery 
coe eon, cr 
Detta pertanto À la lunghezza fondamentale della teoria (connessa al tempo 7, 
dalla relazione 7, = 2R/c), avremo R = (2/3)(e?/mye?), per cui la massa dell’elet- 
trone (evidentemente di natura elettromagnetica) risultera data da m, = (2/3)(e?/Re?). 
Si osservi come, considerando tutti i termini dello sviluppo di v(¢— +), la (4) 
equivalga ad un’equazione differenziale in cui compariscono derivate di ordine 
qualsiasi (°). 


() Ad un’equazione alle differenze finite, equivalente alla (4), & arrivato anche C. I. ELIEZER 
(Report Intern. Conf. on Elementary particles, Bombay, 1950), partendo da un modello (non relati- 
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Procediamo ora alla generalizzazione dell’equazione (4) in modo da renderla 
relativisticamente invariante. Tenendo conto che F, = (ejc)F,gug, avremo che il 
primo membro della (4) dovra essere sostituito da una nuova espressione X,, che 
dovra soddisfare alle seguenti condizioni: 


a) per v2/e +0 le prime tre componenti di X, si dovranno ridurre a 
(mo/T)[v(t) — v(t — To)]: 


b) per „— 0 si dovrà avere X, > (2e2/3c)[(d2u,/ds?) + wyug(d?eg/ds*) |; 
c) essendo P,gwgu, = 0, dovra essere X „u, = 0. 

Queste condizioni determinano la X, che risulta: 
Mo 


X, = — {u,(t) — g(t — To) + ulT)ug(?) [ wg(t) — uglr To)]} = 


To 


re = Gs — To) + Ug(t)Ug(t)Ug(t — 7o)] : 
0 


essendo 7 il tempo proprio s/c. 
L’equazione relativistica alle differenze finite per il moto dell’elettrone si seri- 
ve quindi 


. Mo Cha 
(5) 5 [walt — To) + Ug(t)Ug(t)U_(t — To)] ok F pug =F, - 
0 


Ammesso poi di poter scrivere l’equazione del moto sotto la forma consueta 


du, 
me — = Ff, te: 
ds 


confrontando con la precedente si ha la seguente espress one della forza di reazione 


€ = an a aa 
ds Le 


du,  Ux(T — To) + UaT)Up(T)UplT — To) 
R, = m, - —— — = 
L'integrazione dell’equazione (5), 0 anche di quella corrispondente non relativi- 
stica (4), è in generale molto complicata. Si puö pero osservare che la legge del moto 
che se ne ricava & completamente determinata assegnando il valore di due costanti 


vistico) di elettrone caratterizzato dalle equazioni: 


dv 2 e? dy À GE (2R/c) ae 
Mo — — } >= DK : 


dt SCIE 3 OR 4 (n + 1)! ad” 


che formalmente possono scriversi: 


Moe (1 Je d | ee 
ex wt) =F. 
eR 4 “at 


Quest’ultime, per le note proprietä dell’operatore exp [— »a/dt)] ‚e cioè exp [—A(d/dt)] JE) =f(t — A), 
risultano equivalenti alle (4). 
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a priori arbitrarie e cid corregge il difetto più grave della vecchia teoria (7). Finora 
siamo riusciti ad integrare l’equazione non relativistica (4) in aleuni casi particolari. 
Eccone i risultati: 


1) Particella libera (F = 0): 
v(t) = v(t —t)) = costante. 


Il moto & rettilineo e uniforme. 


2) Particella sottoposta a forza costante (F = K) 


v(t) —v(t—t) K 
= = costante . 
ta My 


Il moto & uniformemente accelerato. 


3) Particella sottoposta a una forza dipendente dal tempo [F = F(+)]: 


To 


(NT) = Do + —- >; Fin). 
OEL 


In particolare: 
Se F(t) = K exp [— at]: 


To _ 1 —exp [— ant] 
vlnr) = V — K ————_ — . 
(ax) ot My exp [ar] — 1 


Se F(t) = K sin ot: 
To 77 COS WT)— COS ((2n+1)/2) WT 
Mo 2 sin (w7y/2) : 


Dr) =o + 


4) Particella in un campo magnetico costante (F = ev A H,): 


2e2H2\—"/2 tT eH, 

| sce) = Vo 1+ 70° 70 cos In arctg | Q à] 
ms Mo 

= œi + yj con 

2e2F72\ — 2/2 teH 

| Gio) ON De sin |n aretg =) 3 

| mr Mo 

5) Particella sottoposta a wna forza elastica di richiamo (F = — ka). 


La soluzione della (4) risulta: 


Al) = 


© (1) Questa affermazione & vera in quanto il postulato 2) ci autorizza a escludere le soluzioni 
periodiche dell’equazione It) — tr) = 0. Cid è in relazione anche al fatto che, assunto come 
origine dei tempi ¢= 0, le funzioni sono definite solo per i valori di f= NT. 
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essendo « una costante complessa determinata dall’equazione trascendente: 


Tok 
a( 1 — exp [— ar) = ——3 
Mo 


ne consegue che il moto dell’elettrone & oscillatorio smorzato. 


6) Blettrone lanciato lungo una retta contro un protone fisso 


e? f HT RC al) 
F = — — e(a) con (a) — > 
a? = NT ENT, 


Pur non essendo riusciti a risolvere esplicitamente l’equazione (4) corrispon- 
dente, crediamo tuttavia di poter affermare, in base a una discussione qualitativa 
di questa equazione, che il moto risulta ancora oscillatorio smorzato e che conse- 
guentemente non si hanno le soluzioni non fisiche che caratterizzano le corri- 
spondenti equazioni (2) e (3). 


Aggiunta in bozza. 


Il Prof. WATAGHIN ci ha cortesemente segnalato un lavoro di F. MôüGricH 
e R. Romer in Zeits. f. Phys., 118, 740 (1939), nel quale si studia il problema del 
moto dell’elettrone mediante un’equazione (non relativistica) alle differenze finite. 
Tale equazione, assai più complessa della (4), è basata su un’ipotesi di «struttura. 
granulare » del tempo essenzialmente diversa del postulato 2) da noi assunto. 
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Interaction of z+ Mesons in Photographie Plates (I). 


A. Mincuzz, G. Puppr and A. Ranzı 


Istituto di Fisica dell’ Universita - Bologna 


Istituto Nazionale di Fisica Nucleare - 


Sezione di Padova 


(rieevuto l’11 Novembre 1953) 


In a systematic survey of photo- 
graphic plates exposed to the 62 MeV 
a+ mesons beam of the Nevis cyclotron, 
besides stars, inelastic and elastic scat- 
tering, which will be discussed in forth- 
coming papers, some peculiar events 
have been observed which will be re- 
ported here. 


a) Scattering of m+ on H. — 22 of 
such events have been selected both 


«scanning along the track» and by 


«area scanning ». The scatterings on H 


5 5 
An 
ae DD. 
ee MAL 
0° 30° 60° Fo 120° 150° 180° 
Ls, 


Fig. 1. — Number of (rt, H) scattering versus 
angle in lab. system 


were identified owing to the characte- 
ristics of the two body collisions (1). 
We have plotted in fig. 1 the angular 
distribution as a block diagram showing 
the number of scatterings per 10° of inter- 
val in the lab. system and in fig. 2 the 


() G. GOLDHABER: Phys. Rev., 89, 1187 


(1953). 


relative differential cross-section in c.m. 
system; i.e. the quantity (47/0) -(do/dQ)- 
as a function of the angle of scattering 
in the c.m. system. 


3 
4 
alan 
2| 
1 
90° 0° 
Dem 
Fig. 2. — (nt, H) differential scattering cross- 


section in c.m, system. 


It is easy to see, in both represent- 
ations, that the angular distribution is 
strongly peaked backward. This result 
is in agreement with similar findings at 
75 MeV by GOLDHABER (?) and at 47 MeV 


(2) Private communication. We are grateful 
to dr. GOLDHABER for many information on his 
work. 
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by Lorp and WEARER (*) with photo- 
graphic plates, but is at variance with 
the results of the Brookhaven group at 
53 MeV in diffusion cloud chamber (4). 
On the other hand, recent experiments 
with counter at 58 MeV by the Columbia 
group (5) and at 40 MeV by the Rochester 
group (5) show a similar strong forward- 
backward asymmetry. 

STEINBERGER and coworkers (7) have 
discussed the phase shift analysis of the 
scattering data at 58 MeV. Since the 
statistical accuracy of their data is much 
greater and our angular distribution 
agrees with their's, we do not attempt 
here to repeat the calculations. 

The distribution as shown in fig. 2 
is independent of the knowledge of the 
absolute cross-section, but we can try to 
estimate this later . Out of the 22 scat- 
terings, 5 were found «scanning along 
the track» giving directly a mean free 
path in G.5 emulsion of 1039 + 460 cm. 
The remaining 17 events were found 
«by area scanning» as a by-product of 
a systematic scanning for elastic scat- 
tering, with projected angle of more 
than 10° in which 541 stars -+inelastic 
scatterings were also found. Since the 
mean free path for all incoherent pro- 
cesses was found, «scanning along the 
track», to- be of 30.8 + 2.2 cm, we 
can calculate, on the base only of the 
relative frequency of the events an 
«equivalent » mean free path in the 
scanning by area, for r++H scatterings 
of 1054-2255 cm. - 


&) We thank Dr. Lorp and Dr. WEARER 
for sending us their unpublished results. 
"We wish to thank also Prof. Fermi for infor- 
mation on the progress of the Chicago work. 

(*) E. ©. FowLER, W. B. FOWLER, R. P. 
SHUTT, A. M. THORNDIKE and W. L. WHITTE- 
MORE: Phys. Rev., 91, 135 (1953). 

) D. BODANSKY, A. SACHS and J. STEIN- 
BERGER: Phys. Rev., 90, 996 (1953). 

() J. PERRY and C. E. ANGELL: Phys. 
Rev., 91, 1289 (1953). 

(7) A. BODANSKY, A. SACHS and J. STEIN- 
BERGER: Phys. Rev., 90, 997 (1953). 


A. MINGUZZI, A. PUPPI and A. RANZI 


The two data combined give a total 
mean free path in our G.5 emulsion of 
1050 + 220 em, which corresponds to an 
uncorrected value for the cross-section 
of 29 mb. This result must be corrected 
for the moisture content that looks to 
be very high on the base of the non 
uniform distribution of the events of 
interest inside the plates. 

Taking a 50% of moisture content as 
a prudential figure, we obtain a cor- 
rected value of 22+5 mb for the cross- 
section at a mean energy of 65 MeV. 
This value is in agreement with the most 
recent determinations. 


b) Stops of m+. — In a total of 712 
incoherent processes as stars and in- 
elastic scatterings observed, 6 stops in 
flight were identified, two of which by 


_scanning along the track. BERNARDINI 


and cow. (8) found none on 115 stars + 
inelastic scatterings. Owing to the 
differences in the statistics the two results 
can agree. 

Supposing that these events repre- 
sent charge exchange phenomena inside 
the nuclear matter, with the exit of a 7° 
and without visible excitation we can 
estimate the elementary cross-section per 
neutron. 

Assuming a mean free path in nuclear 
matter for m+ as well as for 7° equal 
to 2.5075; the probability that a 7° 
escapes absorption before leaving the 
nucleus is 0.32 for light elements and 0.16 
for the heavy ones (°). Taking into ac- 
count that the relative number of inter- 
actions in light and heavy elements is 
0.25 and 0.75, it is possible to calculate 
a mean free path for charge exchange 
phenomena in nuclear matter of A30r, 
that corresponds to an elementary cross- 
section for the process nt +N>P+r° 
of roughly 5 mb. 


() G. BERNARDINI and F. Levy: Phys. Rev., 
84, 610 (1952). 

() E. CLEMENTEL and G. Pupri: Nuovo 
Cimento, 10, 197 (1953) 


INTERACTION OF T+ MESONS IN PHOTOGRAPHIC PLATES - I 


It is interesting to note that this 
magnitude for the cross-section inside the 
nuclear matter is not in disagreement 
with the value on free nucleon, as can 
be deduced, by charge symmetry argu- 
ments from the measured cross-section (1°) 
for the process 7 +P>N-+7°, also 
taking into account on one hand the inef- 
ficiency in the scanning by area for this 
kind of phenomena and on the other hand 
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the possibility that some of the stops ob- 
served (at least 1 with 4 low energy 
electrons) is due to starless capture. 


We wish to thank warmly Pro- 
fessor BERNARDINI, who provided us 
with the plates. 


@°) A. ROBERTS and J. 
Rev., 90, 951 (1953). 


TINLOT: Phys. 


UB Rie RO BV ated 


K. MAYRHOFER — Inhalt und Maß, 
Mit 17 Figuren. VIII. 269 Seiten. 
Springer. Wien, 1952. 


La iettura di questo libro permette di 
valutare i grandi progressi che la teoria 
assiomatica della misura degli insiemi ha 
compiuto a partire dall’epoca, ormai lon- 
tana (1918), in cui C. CARATHÉODORY ne 
tracciava le linee direttive. Ma, mentre 
il CARATHÉODORY limitava la sua teoria 
agli insiemi di punti dello spazio euclideo 
n-dimensionale, l’Autore invece, ad imi- 
tazione di altri, procede a un’ulteriore 
generalizzazione, in quanto gli insiemi da 
Lui considerati sono insiemi di elementi 
appartenenti ad uno spazio astratto S 
arbitrariamente prefissato. 

L’Autore studia, con grande preci- 
sione di dettagli, dapprima (Cap. I) le 
proprieta d’una funzione d’insieme i(A), 
non negativa (e nulla sull’insieme vuoto), 
chiamata « Inhaltsfunktion » 0 semplice- 
mente « Inhalt » (= contenuto), definita al 
variare di A in un corpo k d’insiemi, cioè 
in una famigha k d’insiemi tale che, per 
ogni coppia d’insiemi A,, A, apparte- 
nenti a k, anche la somma A,+4,, il 
prodotto A,A, e la differenza A, — A, 
appartengono a k. Per definizione, i(A) 
è additiva su k cioè tale che 


HA, Ay) = SA) 


(qualunque siano A, e A,, disgiunti, 
in hk). 


E RECENSIONI 


Condizioni particolari, in tale studio, 
vengono successivamente introdotte, in 
quanto da un lato si precisa k in un 
c-corpo, cioè in un corpo tale che, se 
{A,} (n=1,2,...) è una qualunque suc- 


cessione d’insiemi di k, anche la somma 
1.0 


> 4, appartiene a k, d’altro lato i(A) 

n 

è completamente additiva in hk, cioè 
1,0 1,0 A 

tale che (D 4») = dil Ay) (supposti gli 


A, due a due disgiunti). In queste ipo- 
tesi particolari, i(4) è chiamata una 
« Maßfunktion » o semplicemente « Maß » 
(= misura) e indicata col simbolo m(A). 
Per un qualunque insieme B d’elementi 
di S, con la sola condizione che esista 
un insieme A (del corpo k) che lo con- 
tenga, vien definito il contenuto esterno 
i(B) = estremo inf. di i(A), al variare 
di A contenente B, e il contenuto interno 
iB) = estremo sup. di 7(A), al variare 
di A contenuto in B. Se 7(B) = i(B), 
pud porsi i(B) = i(B) = i(B), ottenen- 
dosi cosi il prolungamento del contenuto 
i(B) su una certa famiglia d’insiemi B, e 
tale famiglia è anch’essa un corpo. Ana- 
loghe definizioni e analogo prolungamento 
vengono studiati per una misura m(B), al 
variare di Bin una certa famiglia di sot- 
toinsiemi degli insiemi A (del o-corpo k). 

L’applicazione della teoria generale 
del Cap. I, al caso che S sia uno spazio 
euclideo n-dimensionale, conduce alla 
teoria della misura (contenuto), meglio 
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detta «estensione » dagli autori italiani 
(M. PICONE), secondo PEANO-JORDAN 
(Cap. II), a quella secondo BOREL, infine 
a quella secondo LEBESGUE (Cap. REN). 
Queste teorie ormai classiche, troppo 
note perchè valga qui la pena di ricor- 
darle, sono il fondamento della teoria 
moderna dell’integrazione: esse vengono 
qui studiate con grande profondità, con 
riechezza di dettagli e non trascurando 
di dare, nei punti piü importanti, esempi 
illustrativi accuratamente scelti. Il Capi- 
tolo IV Studia le proprietà che presen- 
tano il contenuto (l’interno e l’esterno) 
di PEANO-JORDAN e, più ancora, la mi- 
sura (linterna e l’esterna) di LEBESGUE, 
per effetto d’una qualunque trasforma- 
zione lineare. Il Cap. V & dedicato ad un 
confronto eritico fra la teoria dell’Autore 
e quella esposta in altri classici trattati. 
Nel Cap. VI sono dati gli elementi del- 
l’algebra di BOOLE, nei suoi sviluppi più 
moderni. Un’appendice, dedicata ad una 
ampia trattazione degli insiemi di Borel, 
chiude questo attraente libro, veramente 
ammirevole per chiarezza d’esposizione, 
per tono scientifico elevato, per accura- 
tezza nei dettagli e per aggiornamento sui 
progressi più recenti. > 

M VIOLA 


J. B. BIRKS — Seintillation counters, 
pp. vur-149, Pergamon Press, 
London, 1953. 


Il contatore a scintillazione, uno dei 
primi strumenti di ricerca nel cam po della 
fisica nucleare, dopo una lunga parentesi 
durata una quindicina di anpi durante 1 
quali era stato sopravvanzato, nella sua 
primitiva versione visuale, da altri pit 
efficienti rivelatori di particelle, à tornato 
oggi nei laborateri. 

Questo ritorno è dovuto in gran parte 
alla odierna possibilità di costruire i foto- 
moltiplicatori, rivelatori di debolissimi 
impulsi di luce, aventi elevate caratteri- 
stiche teeniche; oltre alla scoperta delle 


1757 


proprietà seintillanti di cristalli e di so- 
luzioni. 

Alle possibilità di questa nuova te- 
cnica, ai risultati raggiunti, allo studio 
del meccanismo di emissione di luce da 
parte degli seintillatori & dedicato il libro 
di BIRKS. 

I vari capitoli, a cominciare dalla in- 
troduzione storica, tendono a dare in 
suceinto tutte le possibili informazioni a 
carattere applicativo, sia sui fotomolti- 
plicatori e gli scintillatori presi a se, sia 
in relazione ai possibili usi. 

Nella larga parte del libro dedicata 
al meccanismo con cui avvengono le scin- 
tillazioni, à da notare l’esposizione del 
notevole contributo portato dall’Autore 
alle ricerche in questo campo. 

Gli altri pregi del libro, stampato a 
tempo di primato, vanno dalla estesa 
bibliografia al buon numero di informa- 
zioni, dati sperimentali e valutazioni in 
breve dei limiti attuali della tecnica, co- 
modissimi per chi deve impratichirsi con 
questo strumento di ricerca. 

Un’ultima osservazione riguarda la 
simpatica dedica del libro che è bene 
trascrivere direttamente, perché qualsiasi 
traduzione la guasterebbe: «This book 
is dedicated to “et al.”, that great team 
of scientists, who in the 
brevity, remain anonymous (except in 
the bibliography) D. 


interests of 


G. CORTELLESSA 


O. KLEMPERER — Electron Optics. 
2nd edition, xut-+471 pp., 167 fig. 
Cambridge Univ. Press, Camb- 
ridge, 1953. 


Questo trattato, che fa parte della 
Collezione Cambridge Monographs on 
Physics, è l’estensione € praticamente il 
completo rifacimento di un primo volu- 
metto scritto fin dal 1938 quando, come 
dice l’Autore nella prefazione, «Vottica 
elettronica era al primo stadio del suo 
sviluppo ». 


Sempre secondo l’Autore il libro do- 
yrebbe servire «al duplice scopo di intro- 
durre lo studioso al particolare argomento 
in oggetto e di essere guida al ricercatore 
o al progettista in questo campo, fornen- 
dogli i principi fondamentali, le princi- 
paliinformazioni quantitative e una guida 
alla letteratura specializzata dove pos- 
sono essere trovati tutti i particolari che 
lo interessano ». 

In realta solo il secondo scopo pare, 
e brillantemente, conseguito. Il libro ha 
un’impostazione essenzialmente pratica; 
mentre le dimostrazioni sono spesso solo 
accennate, le formule finali delle analisi 
quantitative sono ampiamente discusse, 
di esse viene quasi sempre indicato lo 
specifico campo di utilità nelle pratiche 
applicazioni, e fornito un accenno dei 
metodi sperimentali mediante i quali la 
loro validità pud essere controllata. Inol- 
tre con numerosi esempi concreti si inse- 
enano a sfruttare praticamente gli stru- 
menti di lavoro forniti dai principi gene- 
rali. Questi ultimi sono invece ricavati 
in modo piuttosto sommario ed in questo 
senso non pare che il libro possa essere 
utilmente diretto verso chi per la prima 
volta abbordi l’argomento e voglia criti- 
camente assimilarlo. 

Il 1icercatore troverà invece nel libro 
un ausilio prezioso per la mole di infor- 
mazioni, e per le ricchissime e precise 
citazioni bibliografiche, estese a 500 voci. 

E inoltre una gradita novità quello 
di veder unifieati in un solo trattato i 
campi dell’ottica elettronica vera e pro- 


pria, della dinamica delle particelie 
quando la carica spaziale nor possa venir 
trascurata (argomenti del noto libro 


Theory and design of electron beams del 
PIERCE), della spettrografia di massa, dei 
tubi elettronici, etc. 

Il libro & diviso in dodici capitoli. 

Il primo introduce brevemente i fon- 
damenti teorici dell’ottica elettronica. 

Il secondo tratta degli elementi cardi- 
nali di una lente nell’approssimazione 
gaussiana e i metodi sperimentali per la 
loro individuazione. 


LIBRI RICEVUTI E RECENSIONI 


Il terzo riassume i numerosi metodi 
per il tracciamento dei campi e per la 
costruzione per via grafica o numerica 
della traiettoria dell’elettrone in essi- 
Chiude il capitolo un brevissimo cenno 
sui metodi di «sintesi» degli elettrodi 
capaci di impartire all’elettrone una 
traiettoria assegnata. 

Nel quarto capitolo vengono illustrate 
le proprieta delle più comuni lenti elettro- 
statiche. 

Nel quinto, dopo l’introduzione del 
potenziale di Stormer, si passano in ras- 
segna le proprieta delle più comuni lenti 
magnetiche. 

I capitoli sesto e settimo sono dedi- 
cati alle aberrazioni, trattate da un punto 
di vista essenzialmente sperimentale. 

L’ottavo e nono capitolo sono dedi- 
cati alla carica spaziale nell’ottica elettro- 
nica e particolarmente sviluppato & il 
progetto di intense sorgenti di elettroni. 

Il decimo capitolo è dedicato all’ot- 
tica dei sistemi bidimensionali caratteriz- 
zati da linee anziché da punti focali. Ven- 
gono discussi anche in questo caso gli 
effetti dovuti a carica spaziale e il pro- 
getto di sorgenti lineari di elettroni. Le 
considerazioni svolte sono applicate a 
tubi a raggi X, a tubi elettronici a fascio 
e a fotomoltiplicatori. Chiude il capitolo 
un cenno sul metodo dell’« analogia in- 
versa» che permette mediante un mo- 
dello in vasca elettrolitica la determina- 
zione di campi magnetici bidimensionali 
a vorticita non nulla. 

L’undicesimo capitolo, tratta dei cam- 
pi dettrici e magnetici di deflessione e 
delle loro proprieta focalizzanti. Le de- 
duzioni sono applicate alla spettrografia 
di massa, alle deflessioni .del fascetto 
elettronico nei tubi oscillografici, e alla 
stabilita delle orbite negli acceleratori di 
particelle. 

Il capitolo dodicesimo da un quadro 
descrittivo finale delle più importanti 
applicazioni dell’ottica elettronica nella 
ricerca e nell’industria. 

L’edizione a cura della Cambridge 
University Press & molto curata e due 
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indici oltre l’imponente bibliografia ren- 
dono molto agevole la consultazione del- 
l’opera. 

EMILIO GATTI 


F. O. Rick e E. TELLER-La struttura 
della materia. Edizioni Einaudi, 
Biblioteca di cultura scientifica, 
vol. XX XVII, Torino, 1953; pagg. 
543. 


Nel corso dei due ultimi decenni, pa- 
rallelamente al progressivo sviluppo della 
meceanica quantistica, è via via divenuto 
possibile rendersi conto di una parte 
sempre maggiore dei fenomeni chimico- 
fisici connessi con la struttura atomica 
della materia. Anche nei casi — invero 
tutt’altro che sporadici — in cui il pro- 
blema fisico si traduce in uno matema- 
tico troppo difficile per essere risolto in 
modo esatto, & possibile di regola gettar 
luce sugli aspetti essenziali e più inte- 
ressanti del problema, attraverso l’im- 
piego di opportuni metodi di approssima- 
zione. Perd, causa l’enorme difficolta del- 
l’algoritmo matematico richiesto per la 


risoluzione sia pure approssimata della. 


maggior parte dei problemi concernenti 
la struttura della materia, studi di questo 
tipo sono finora rimasti relegati, esclu- 
sivamente o quasi, al mondo degli spe- 
cialisti. 

Con questo libro « sulla struttura della 
materia» gli AA. si propongono eviden- 
temente di rendere accessibili ad una 
classe più vasta di lettori, le conoscenze 
fino ad oggi acquisite in questo campo 
della fisica. Il libro infatti & seritto quasi 
senza fare ricorso all’algoritmo matema- 
tico, riportando solo occasionalmente delle 
formule atte a mostrare l’andamento dei 
fenomeni discussi, indicando soltanto le 
linee essenziali dei ragionamenti che con- 
ducono ai risultati finali e facendo fre- 
quentemente appello aliintuizione, che 
qua e là vien sorretta da diagrammi € da 
rappresentazioni grafiche. 
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Oltre ad un breve capitolo int rodut- 
tivo (in verità troppo breve per «intro- 
durre » il lettore inesperto ai concetti 
della fisica quantistica...) opera com- 
prende ancora dodici capitoli. Il primo 
di questi contiene una rapida trattazione 
dell’atomo di idrogeno dal punto di vista 
della meccanica quantistica. I concetti 
di « autofunzione di stato », di «stati de- 
generi », di «ortogonalita » vengono qui 
introdotti (forse di nuovo in modo un 
po’ troppo succinto per il lettore che con 
essi non abbia alcuna familiarita) per 
essere poi utilizzati ampiamente nel corso 
del libro. Nel capitolo successivo, dedi- 
cato al sistema periodico degli elementi, 
si introduce il principio di Pauli (dap- 
prima per il caso semplice in cui si tra- 
scura lo spin dell’elettrone, e poi nella 
sua formulazione più generale) metten- 
done chiaramente in luce aleuni degli 
aspetti pit significativi, quali Videntita, 
che da esso consegue, tra tutti gli elet- 
troni. Seguono aleuni capitoli nei quali 
si presenta in forma sobria e piuttosto 
rapida quanto di essenziale si conosce 
sul moto dei nuelei nelle molecole, sul 
comportamento degli atomi e delle mole- 
cole nei campi elettrici, sulle proprieta 
magnetiche della materia e sulle forze 
attrattive di van der Waals. Con molta 
maggior diffusione sono invece trattati 
gli argomenti concernenti il legame chi- 
mico (dove vengono illustrati e discussi 
i metodi degli orbitali molecolari e degli 
orbitali di valenza), le forze agenti nei 
solidi, le vibrazioni molecolari (spettri 
infrarossi, spettri rotazionali e vibrorota- 
zionali, effetto Raman, ecc.) e gli spettri 
elettronici. I] penultimo capitolo del 
libro, piuttosto ampio ma proporzionato 
alla vastità dell’argomento in esso trat- 
tato, & dedicato alla chimica nucleare. 

L'opera termina con una rapida ma 
assai suggestiva descrizione dello stato 
della materia nelle stelle, e con un cenno 
alla teoria dell’universo in espansione e 
alle possibili spiegazioni sull’origine degli 
elementi. 

La sensazione che si ha leggendo il 
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libro — e che non puö essere molto smi- 
nuita dalla natura degli argomenti trat- 
tati nell’ultimo capitolo — è che gli AA. 
abbiano voluto con esso perseguire anche 
un fine pratico: quello cioè di avvicinare 
al mondo della fisica quantistica i ricer- 
catori ed i tecnici che in certe industrie 
(per esempio, chimiche e metallurgiche) 
hanno frequentemente occasione di im- 
battersi in problemi connessi con le pro- 
prieta della materia. Come se gli AA. 
memori degli immensi benefici che da 
mezzo secolo a questa parte industria 
ha tratto dalla corretta applicazione delle 
legoi dell’elettromagnetismo, si fossero 
proposti di preparare il terreno sul quale 
i «nuovi ingegneri» dovranno operare 
applicando le leggi della nuova mecca- 
nica. 

Si deve pero osservare che a dispetto 
del tono sobrio e discorsivo con cui à 
scritta quasi ogni parte del libro, la let- 
tura di aleuni eapitoli è tutt’altro che 
facile per chi non sia in una qualche 
misura iniziato agli argomenti in essi 
trattati. 

Di queste difficolta di lettura sono 
forse responsabili, principalmente, le dif- 
ficoltà intrinseche dei numerosi concetti 
che inevitabilmente occorre introdurre 
nella stesura di un’opera di questo ge- 
nere e la vastità degli argomenti discussi. 
Ma qualunque sia la loro causa, resta il 
dubbio che lo studio dell’opera pssa es- 
sere per il lettore non fisico utile nella 
misura che sarebbe desiderabile ai fini 
del superamento di quella barriera tuttora 
esistente tra il mondo delle applicazioni 
ed il mondo della fisica quantistica. 

Il libro è invece indubbiamente assai 
istruttivo per il giovane fisico e per il 
chimico-fisico che, già in possesso delle 
idee fondamentali della meccanica quan- 
tistica, troverà in esso un quadro vasto 
e completo delle attuali conoscenze sulla 
proprieta della materia. 


M. CONVERSI 


w. GLASER- Grundlagen der Elektro- 
nenoptik. Un vol. in-8° di x-699 pa- 
gine, ed. Springer, Wien 1952. 


I trattati di ottica elettronica apparsi 
nel ventennio di vita di questa giovane 
scienza si dividono (assai più nettamente 
che quelli di ottica della luce) in due cate- 
gorie abbastanza distinte, secondo che il 
loro autore ebbe in vista principalmente 
l’aspetto tecnico o quello fisico-matema- 
tico della materia. Il volume del GLASER 
è uno dei migliori rappresentanti, e certo 
il più completo, di quest’ultimo indirizzo. 
Esso si propone di esporre in modo siste- 
matico e rigorosamente deduttivo i me- 
todi di calcolo delle traiettorie e delle 
onde elettroniche, cosi da mettere in 
erado lo studioso di leggere criticamente 
la letteratura specializzata e da fornire 
al progettista tutto cid che occorre per 
il calcolo dettagliato di uno strumento 
ottico-elettronico. Le descrizioni di appa- 
recchi sono ridotte al minimo e confinate 
a un capitolo introduttivo. 

Nonostante questo carattere altamen- 
te fisico-matematico e l’esigenza di rigore 
continuamente tenuta presente dall’Au- 
tore, il libro & di lettura relativamente 
facile. Il GLASER ha infatti avuto cura di 
esporre i principii con grande chiarezza, 
senza presupporre nel lettore speciali co- 
enizioni di fisica o di matematica che va- 
dano oltre (per intenderei) quelle for- 
nite ai nostri studenti nel primo biennio 
di fisica; inoltre ha svolto le deduzioni 
matematiche senza economia di spiega- 
zioni dei passaggi analitici e senza evitare 
ripetizioni ove possano facilitare la com- 
prensione e risparmiare fastidiosi rinvii. 

Lo spazio risparmiato sacrificando la 
parte descrittiva e strumentale, nonchè la 
mole considerevole del libro, consentono 
al GLASER di svolgere la teoria dell’ottica 
elettronica con un’ampiezza assai mag- 
giore del consueto. Si trovano percid in 


questo volume molti argomenti che non 


avevano ancora avuto sistemazione trat- 
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tatistica, tra cui non pochi contributi 
originali dell’Autore. 

Il volume & diviso in tre parti. La 
prima, dopo due capitoli sui principii 
generali dell’ottica elettronica, tratta 
esaurientemente, nell’approssimazione di 
Gauss, l’ottica dei campi elettrici e ma- 
gnetici con un asse di simmetria. A ca- 
ratterizzare l’indole del trattato, si puo 
rilevare, ad esempio, che non solo sono 
stabilite e discusse le varie forme delle 
equazioni delle traiettorie parassiali, i 
loro metodi di integrazione e le loro ap- 
plieazioni a diversi tipi dilenti, ma inoltre 
i concetti di «fuoco » e « distanza focale » 
di una lente elettronica sono discussi cri- 
ticamente, e sono trattate a fondo varie 
questioni ad essi relative che non hanno 
Vanalogo nell’ottica ordinaria @ che di 
solito vengono sorvolate nei trattati di 
ottica elettronica. Particolarmente svi- 
luppata, come & naturale, & la tratta- 
zione di quella speciale schematizzazione 
delle lenti magnetiche che è generalmente 
nota col nome dell’Autore e che si è rive- 
lata di tanta pratica utilità. Un capitolo 
speciale è dedicato alle questioni di 
«apertura » nel microscopio elettronico ; 
un altro alla aberrazione cromatica. 

La seconda parte del libro riguarda 
la teoria delle aberrazioni geometriche 
(del terzo ordine) che sono trattate siste- 
maticamente col metodo hamiltoniano. 
In questa parte si trovano anche i soli 
due capitoli del libro dedicati a campi 
agenti essenzialmente in direzione tra- 
sversale al moto: uno di essi tratta la 
teoria della (debole) deflessione elettrica 
o magnetica di un fascetto elettronico 
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quale si verifica, per esempio, negli oseil- 
lografi catodici, e le aberrazioni che essa 
introduce nel focheggiamento del fa- 
scetto ; l’altro riguarda invece il focheg- 
giamento di fasci elettronici ad asse eir- 
colare. Peccato che l’Autore tratti questi 
due problemi solo nella approssimazione 
non relativistica, e tralasei inoltre la 
teoria generale dei campi deflettori che 
si applica, per esempio, nella spettro- 
metria di massa. 

La terza parte del libro, dal titolo 
«Ottica Elettronica e Meccanica Ondu- 
latoria », fornisce al lettore un’ampia e 
in gran parte originale trattazione dei 
fondamenti della meccanica ondulatoria 
e della loro applicazione alla formazione 
delle immagini elettroniche. In partico- 
lare, & trattata in dettaglio la teoria del 
potere risolutivo del mieroscopio elettro- 
nieo in relazione sia alla diffrazione che 
alla aberrazione sferica. I capitoli intro- 
duttivi di questa terza parte possono 
interessare il fisico teorico anche indi- 
pendentemente dalla loro applicazione 
all’ottica elettronica. 

Il manoseritto di quest’opera, come 
racconta l’Autore nella prefazione, era 
stato seritto in una prima redazione 
avanti la guerra. Avendone varie circo- 
stanze ritardato a più riprese la pubbli- 
cazione fino ad oggi, l’Autore Vha per 
due volte rimaneggiato e completato. Il 
libro perd non serba quasi traceia di cosi 
complicate vicende, e si presenta ben 
ordinato e generalmente aggiornato. La 
veste tipografica & quella ben nota, niti- 
dissima, della casa Springer. 
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